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Introduction and Summary 


This is the final report on grant number NAG5-1 96 which started 
on July 1, 1981. 

Theses and Dissertation 

The grant has provided partial financial support for three 
MS students and one PhD student. The names of the students and the 
titles and abstracts of their theses and dissertation are listed 
below. 

D. D. Dyche, MS thesis entitled "Experimental Determination of 
Atmospheric Scattering Effects on Scanner Edge Response." 

The Influence of background radiance from surface 
areas outside the ground-projected IFOV of an aircraft 
sensor was experimentally Investigated on a 600 x 550 m 
field, evenly divided between mature cotton and bare 
soil. An aircraft carrying a multlspectral scanner made 
passes at 1,000, 8,000, and 16,000 ft A6L obtaining 
radiance data In the 420-450 nm range. Ground 
reflectances were approximately 0.025 (cotton) and 0.105 
(bare soil). At 440 nm the atmospheric extinction 
coefficient was measured to be 0.297. The data indicate 
that a boundary effect may exist, where flux reflected 
from one side of the field is scattered by the 
atmosphere so as to appear to the sensor to come from 
the adjacent side. Experimental and theoretical data 
agreed closely for the cases of total detected radiance 
changes and border radiance values but not for radiance 
changes between altitudes. Thus, a boundary effect may 
be present, but its existence is not clearly proved or 
disproved. 

A. L. Phillips, MS thesis entitled "Absolute Calibration of and 
Atmospheric Measurements with a Multiband Field Radiometer." 

Three different techniques were employed to 
calibrate in an absolute sense a multiband field 
radiometer. The three techniques were a source-based 



one using a calibrated standard lamp as a reference, a 
detector-based one using a collimated monochromatic 
source, and a detector-based one using an extended 
uniform source. Because of the optical system of the 
radiometer, the detector-based technique using the 
monochromator system was Inappropriate for absolute 
calibration. Based on the uncertainties associated with 
the other two techniques, the detector-based one using 
the extended source was preferred. However, the 
monochromator system was also used to get the shape of 
the responslvlty as a function of wavelength. This 
technique 3 rlelded an absolute calibration with an 
uncertainty of ±4.5 Z. 

Using the same radiometer, an experiment was 
conducted to teat the influence of “surround spectral 
reflectance” upon the measured ratio of the diffuse 
component of solar radiation to the direct component. A 
comparison was made of the measured ratio to the ratio 
determined by an atmospheric radiative transfer program. 
The comparison showed similar trends in 3 visible bands, 
but showed disagreement In a near Infra-red band. 


S. L. Wltman, MS thesis entitled “Hadlometrlc Calibration of the 
Thematic Mapper 48-lnch Diameter Spherical Integrating Source 
(48-SlS) Using two Different Calibration Methods.’* 

Two methods are described for calibrating the 
48-lnch Integrating sphere used to radlometrlcally 
calibrate the Landsat Thematic Mapper. The output of 
the sphere was first measured monochromatlcally by 
engineers at the Santa Barbara Research Center. The 
author of this thesis reduced these monochromatic data 
such that they coincided with the spectral bandpasses 
of a calibrated radiometer. The second calibration 
method used this radiometer to measure the output 
radiance of the sphere. The reduced monochromatic and 
radiometric data sets were compared. The results show 
that Che radiance of the Integrating sphere measured 
with a monochromator and a radiometer agreed within the 
limits of Che uncertainties of the measurements. The 
radiometer was calibrated against the same reference 
standard as used for calibrating Che monochromator. In 
addition, Che monochromator versus radiometric sphere 
calibration results show good agreement with other 
independent source-based and detector-based calibrations 
of the radiometer. 


C. J. Kastner, PhD dlssertatloa entitled ”In-*Fllght Radiometric 
Calibration of the Landsat Thematic Mapper." 

The In-fllght absolute radiometric calibration of 
the Thematic Mapper (TM) Is being conducted using the 
results of field measurements at White Sands, New 
Mexico. These measurements are made to characterize the 
ground and atmosphere at the time the TM Is acquiring an 
Image of White Sands. The data are used as Input to a 
radiative transfer code that computes the radiance at 
the entrance pupil of the TM. The calibration Is 
obtained by comparing the digital counts associated with 
the TM Image of the measured ground site with the 
radiative transfer code result. The calibrations 
discussed here are for the first four visible and near- 
infrared bands of the TM. 

In this dissertation the data reduction for the 
first calibration attempts on January 3, 1983, and 
July 8, 1984, Is discussed. Included are a review of 
radiative transfer theory and a discussion of model 
atmospheric parameters as defined for the White Sands 
area. These model parameters are used to assess the 
errors associated with the calibration procedure. Each 
Input parameter to the radiative transfer code Is varied 
from Its model value In proportion to the uncertainty 
with which It can be determined. The effects of these 
uncertainties .on the predicted radiances are determined. 

It Is thought that the optical depth components t^y, 

TMle> '^oz> ^820 measured to within 10%, 2%, 

10 % , and 30 % , respectively. For the white gypsum sand, 
surface reflectance uniformity Is on the order of 1.5%, 
and Che overall uncertainty In measured reflectance Is 
about 2 % . This Is due to an uncertainty In the 
reflectance factor of Che calibration plates. The 
greatest uncertainty In calibration Is attributed to our 
uncertainty In the aerosol parameters. In particular the 
Imaginary component of refractive Index. The cumulative 
effect of these uncertainties Is thought to produce an 
uncertainty In computed radiance, of about 5 % . 


Publications: 

Ten publications were also directly supported by the grant. The 
authors, titles and journals are listed below and the papers 
themselves have been duplicated to comprise Che balance of this 


final report 
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P. N. Slater (1982), "Absolute radiometric calibration of advanced 
remote sensing systems," International Geoscience and Remote 
Sensing Symposium, Munich, Germany. 

P. -N. Slater and R. D. Jackson (1982), "Atmospheric effects on 
radiation reflected from soil and vegetation as measured by 
orbiting sensors using various scanning directions," Appl. Opt. 
21:3923-3931. 

C. J. Kastner and P. N. Slater (1982), "In-flight radiometric 
calibration of advanced remote sensing systems," Proc. SPIE 356. 

P. N. Slater, (1984), "Importance and attainment of accurate absolute 
radiometric calibration," Proc Vol 475. SPIE Critical Review of 
Remote Sensing, ed. P. N. Slater pp. 34-40. 

N-Z Che, R. 0. Jackson, A. L. Phillips, and P. N. Slater (1984), "Field 
radiometer methods for reflectance and atmospheric 
measurements," Proc. SPIE 499, pp. 24-33. 

P. N. Slater (1984), "Radiometric considerations in remote sensing," 
Invited paper In the Proceedings of the IEEE Vol. 73, No. 6, pp. 
997-1011. 

S. F. Biggar, C. J. Bruegge, B. A. Capron, K. R. Castle, M. C. 

Olngulrard, R. G. Holm, R. 0. Jackson, L. J. Llngg, Y. Mao, M. S. 
Moran, J. M. Palmer, A. L. Phillips, R. K. Savage, P. N. Slater, S. 
L. Wltman and B. Yuan (1985), "Absolute calibration of remote 
sensing Instruments," Proc. Third International Colloqlum on 
Spectral Signatures of Objects In Remote Sensing, ESA SP-247, 
pp.309-314. 

P. N. Slater (1985), "Variations In in-f light absolute radiometric 
calibration," ISLSCF Conference In Rome, Italy. 

Ray D. Jackson and Philip N. Slater (1986), "Absolute calibration of 
field reflectance radiometers," Photogrammetrlc Engineering and 
Remote Sensing, Vol 52, pp 189-196. 

The grant has been of Immense value in providing support for 
some basic and applied research aimed at improving our ability to 
determine the in-fllght absolute radiometric calibration of 
advanced land observing systems. In particular it has provided 
research funding in partial support of our very successful 


satellite sensor calibration program at vnilte Sands. We are 
pleased to say that this grant activity Is being continued as part 
of the Earth Observing System program and wish to express our 
thanks to NASA for their continued support of the satellite 
calibration program at the University of Arizona. 
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9. a. Slaear 


CsmecM on itaaoea Sanaing 
Optical Seianeaa Castar 
Onivacaity of Ariaoaa 


ABStSJIc; 

Oio diatinctian bacMaaa eha uaaa of ralatly* 
and ahaoluca apa cec o r a rtl i T a rTt c callbcatioa of ra- 
flnta aaeaind lyacaaa ia diacnaaad. Tha advaataaaa 
of datacoDg- h a a ad aHaolota calihracioa ara daaerUaad. 
aad cba cataqaciaa of ralativa aad abaolnta ay a e a a 
ealibrationa ara LLacad. Ttia Uaitaciaaa aad 
pcodlaBa aaaoeiacad with thraa eoaaoa eathada uaad 
for eha ahaoluta calibration of raaeta aanaiitg 
aystaaa aca diacnaaad. 

TWO oathoda aca pmpnaaii far tha ia-flioht 
abaolnta calibration of advaaead Boltiapaetcal 
Uaaar array (MtAl ayacaaa. Ona oakaa uaa of a sna> 
illiiari natad panal ia f ro nt of tha aanaor. tha 
radlaaea of which ia aonitorad by a spaetrally 
flat pycoalacBie raril naarar. tha other uaaa a 
lacda. uaifoxB. hioh-cadiaaea cafaranca prodad anr> 
faea. tha ^ o on d aad ataoaphacic iriaaaiir»aiira 
raqoirad aa inpnt ta a cadiati«a traaafar pr opr M 
CO pradiet cha cadiaaoa laval ac cha antraaea pnpil 
of tha ottital aanaor ara diacnaaad. aad cha 
qrouad laatmaancatioa ia daaerlbad. 

Kay wocdai Radlaaacryi Calibcaclaai Raaota Sanalac. 

I. SPeCTRORJUUOPBTlUC CKLIBmTrOIl— «80 NZSCS m 

Pet pnrpoaaa of Chla diacnaalon. cha caaaBaity 
of renece Sana Inn data uaara can ba dlwtdad Into 
two ocouca ; 

1. riaara recnirino calativa. hue not 
neeaaaarily abaoluca. spaetcocadlaratrle sanaot 
calihcaclon. Thaaa Ineluda workers in uuuuucer* 
aided scene elasaifleacion. cartoerashers. Inaee 
processors, phocoinesrprsesrs. and oeools conesmsd 
with cooposlno Laroe mosaics. 

2. Osers cecuitlnq absoluts speetso> 
cadlcneeric calibraeion. Thsse Ineluds ptnrsical 
scienrisca ceneamsd with relaelne preand-neasared 
pacaascera and/or aesnsohsrie eharactsrlaclcs to 
the sseetrai radiance at tha entrance pupil of tha 
scacs sanaor. 

The dlatlnetion can ba draen Uiat tha feraer 
are concerned oriaarilv with takinp an I n v en tor y or 
coaoillnp a nap of oronnd feacurea. The lattar are 
concerned priaarily with undaratandlno and 
charactarlzinp the physical interacclcns taklnq 
place, uaiiallv throuqh nodela and often with a view 
to ootiairino (for eaaaple an Irrication sehedulel 
or to oredictine a chanoe in a oarticular process 


(for eaaapla an aorlcultnral- yield) . Expreasad 
differently, these coneamad with relative 
ealibcaeion use data deseribad in dlqltal eounca. 
while thoaa naadlno absoiuca calibraeion use data 
refarancad in taras of radiance unita at the 
aneranea oopil of the sensor (in this case, tha 
radloBatrle calibration of tha sensor la Invokad to 
convert the dlqltal value to spectral radiance) . 
l*a mat o^haaixe that the radiance value ia of no 
aers worth to the averaqo user than tha dtqital 
value, as both rafar to an ineldant radiance laeol 
at tha sanaor over broad, unequal spaetral 
paaabanda. over which the spaetral reflaetanea of 
the feature bainq obsatvad ean ehaiiqa by a 
aiqalfieant but unknown aaeunt. Bowevar. there ara 
two reasons for eonvertinq the dlqltal valua to 
radianeei first, in onltitacperal sansinq, to 
aecoont for anv d oe u aentad ehanqas of radlocatrie 
calibration with tioot saeond. to tast or utlUae 
physical nodela la which thn qrouad raflaeeanea and 
aeaoapharie effaets are naanorad and/or calculated 
over identical spaetral paaabanda as a s qilrur a il by 
the space sensor. 

Perfaet relative fadloaatrie reaponsa occur s 
whan the outputs free all detaeters In a hand are 
equal or ean be adjusted durlnq prepcoce a alnq to ba 
equal, whan the Ineident spectral radiance la 
eenscant aecosa the sanaor 's fltld of view, (tioea 
that the nuaher of detaetera la a band can be aa 
Caw aa sis -for tha Hultlspaetral Scanner Syatasi 
(MSS) on Landsae and aa many as 18,500 on future 
MLh syataas.) This eundltion nust he net 
Indapandancly of the spectral content of tha seane. 
Whan this etnditian is net net, tha laaqa appears 
striped. If It is not s can #' daoe nd ant, stripinq 
often can be cc^lacely cenoved by the hlatooraa 
equalization nathodi thus, relative radlcoatric 
precision can ba hlqh even thouoh the accuracy 
involved nay ba low. In this equallxation 
procedure the hlatooraa of each datactor oueput ia 
ppapared with tiiat of every other detector, after a 
laroe nuaCer of data saaplea (^.2 f 10^) have been 
recorded. It is assuaad that. If tha scene is 
spatially and siiaetrally randoa, tha hlatooraos for 
a laroe nuaber of saaplas will ba identical. If 
the histoqraas ara net identical, adjoataanca are 
nade durlnq the preproeassinq step to naka thea so. 
This proeedore ean be repeated for seenea of 
different averaoe radiance, and tha relative 
responsee ean then be aqoalizad over the dynaaic 
ranqe of the dateetora. This procedure does not 
work if the stripinq is sc en e ■ d e p e n dent (Ref. 1) . 
?ortunately the design and the spectral hand 
location for the Theaatic Mapper and future MLA 
sensors will not qlve rise to a sce n e- dependent 



strioino nroblen o£ ebe .'Jiaonieue'# of shat in the 
MSS. 

Inar'eouatelv corrected relative- deteetot-to- 
detector response causes unsightlv stripina in the 
iaiaqerv. Strioino can cause inaccuracies in 
lueomaced scene classification. In addition, even 
t.houqh the huaan ohotointernreter is usuallv mote 
foroivino than a conpucsr, if the striping is 
severe it can cause errors in ohotointerpretation. 
Clncorrected strioing can be aogravated hv the 
apolication of some imaoe orocessing algorithms. 
Thus, it is imoortant to reduce striping as much as 
possible for image processing and classification 
purooses. The relative response of the individual 
detectors is also important for ratioing purposes. 
However , once the in-band tolerance has been met, 
the concern is with the relative stabilitv of 
response because menv appilcetiana utilise the 
ccmoerlson of multi temporal band ratios. 

The utilization or verification of phyeical 
models usually reouices the use of data calibrated 
in an absolute sense. 7ntil recently the highest 
in-orbit abseiuts radlonetcie accuracy has been 
iittle better than 10%. This low aeenraey has been 
due tot (al the fact thee the callbretion in orbit 
has oftan been far the focal plane onlv, not for 
the camolete syatemi ibi the loss in accuracy 
aceoaoanvino the transfer of calibration from the 
standard source at the netionsl labofaterv to the 
factory or laboratory calibration site? (e) the use 
of soure^bssed calibration proesdnrea. 

As dsserlbad latsr, the use of deteetor-besed 
calibration prosiaea to reduce the appteuiastely 
10% error to 1%. To whet extent we can afford to 
relax or strive to ex c eed the 1% potential can beet 
he decemined from a set of well-coordinated 
meesuresient and modeling exercises conducted as 
part of an experimental MIA niasicn. The results 
of wort to data on this subject are in 
dlseoreenent. On one hand we Vcnow that netursl or 
scene variabilitv is seldom lees thm i% and we 
also know froai atmospharlc modeling (not yet 

experimentally verified) that we can — 
considerable atmospherlc-seatter-indnced spectral 
radiance croaacalk between naighbering ground 
saoBles owino to the so-called adj a c en c y or 
boundary effect (Refs. 2-3) . On the other hand, 
the sensor simulations meda to data (Rsf. 3) 
indicate there is an improve m e n t in classification 
aceuracv in going from S-bit to 8-bit gusntlzatlon. 
Sight-bit ousntlzstion is available on meat future 
sensors. A purist mloht centsnd that, if 8-btt 
radlomatric resolution is really justified, it 
should, for maximum utility, be aasoeiatad with 
eommensnrate (<0.4%) error in abeolute radiometric 
accuracy ( 

In sunssarv, we know that soactroradtomecrie 
calibration is imoortant and that the need for 
celative or abeolute calibration depends on the 
application for the data. How e v e r, we do not know 
how accurate the calibrations should be, 
particularly the abaeluca eallbracion. The most 
eneeuraoing aspect of the situation la that we new 
have the potential of utilizing 3-bit quantization 
resolution with cemmenaurata absolute radiomstrie 
accuracy in an orbital imager at laaat over the 400 
to 900 nm range. This potential should be 
exploited as fullv as possible as an important 
component of anv system and applications research 
associated with an experlmentiu. land-obeerving 
system having an 3- to 10-bit quantization 
caoability. 


2. 0ETECT0H-3A.3SD HAOIOMBTSIC CALIBRATION 

The recent work at the OS National Bureau of 
Standards on self-calibrated photodiodes is 
described only briefly here. For more details the 
interested reader is referred to references 7-12. 

The csLlibration of the ohotodiode is 
accompiished by either of two differenciv applied 
biasing orocedures, deoendlno on the wavelength 
region of interest. At short wavelenochs, a 
negative bias is aoplied to remove the 
re c o mb ination centers at the Sl-SiO interface at 
the front of the detector. To do this, a contact 
is made with the front surface using an eleetzode 
iamersed in a conducting liquid, or the surface can 
he exposed to a corona dlseharoe. For long 
wavelengths, a back-bias is aoplied to extend the 
deelstion region to a depth beyond which incident 
flux penetrates. The experimental procedure- is to 
irradiate the detector with a conatant 
moneehreoatie flux level and to increase the bias 
TOltags until further increase no longer gives rise 
to an increase in output signal. For both the 
short and the long wavelength ranges, the internal 
quantum efficiency sstnratsa at a value extresely 
close to unity, as shown In Fig. 1. Thus the 
maximum Inereaaa in signal output obtained as a 
result of biasing can be used to detsmina the 
intamal quantum effleisney of the detaetor without 
biasing, as it will be used in practice. 



!• Typical pbotodlode Intamal onantum 
affleianey without biasing (lower curve) and with 
biaalng [upper curve) , reference 9. 


The only significant loaa in the external 
quantum efficiency of the photodiode Is caused by 
reflection. This can be red u ced to insiqnifleanee 
by making use of three photodiodes according to the 
geometry sketched in Fig. 2. 



?lg. 2. A three-diode arranaanent to 
minimize specular reflection losses. 


The output signals from the three diodes are 
sunsaed to provide the signal corresponding to a 
black detector of overa.ll guanttsn efficiency that 
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-ran be assuaoc! io he uni tv. The second diode 
reflects t.he specuiac reflection from diode 1 to 
diode 3. which acts as a tettoceflectot. The 
incident flux thereby underooes five reflections, 
and at a 10% value the final specular reflectance 
is dovm to 10“' of the initial incident radiance. It 
is claimed (Ref. 12) t.hat the diffuse reflectance 
losses for clean detectors are tyoically much less 
than 1%. The photodiodes used in similar 
calibration facilities at !fflS and at the aniversity 
of Arizona, Optical Sciences Center, are EGG OV 
•;44B. 

The discussion at the end of this paper refers 
to the use of self-callbratinp MBS detectors for 
the spectroradioaetric calibration of an MIA 
system. Rowever, as shoim in ?iq. 1. the unbiased 
ouantuffl efficiency is wavelencch dependent. 

Because of oossible changes in the pasaband 
oosition of soeetxal filters during long duration 
space flights, it may be advisable to use the MBS 
detectors to calibrate a spectrally flat 
oyroelactrlc detector, at the 0.1-0. 2% level, and 
for that to be used in the in-orbit calibration. 


3. BROAD (XASSZ7ICAT10N OP CALXB8ATX0R paoem O Hg S 

Soto of the most comaonly used proeeducea for 
the calibration of resteta sensing systems are 
referred to in Pig. 3. The two major divisions in 
the fioare are between relative and absolute 
calibration and between the static aaer^ image 
resoonae and ths dynamic micro- iaacs response of 
the system. Some aspects of relative calibration 
were dlsenased in the first pact of this paper and 
together with dynamic micro- image response will not 
be discussed further, beyond remartcing that the 
dynamic mlcro-imMe cesponae Is of vital intareat 
in anv pixal-by-pimel analvsis of remotaly 
eolXeet^ Imsgecy. Ths rest of this paper deals 
with the absoluts calibration of remote sensing 
systems, and only ths static macro- image response 
will be considered in this ccntext, as Is usually 
ths case. 

The ccceedurea for the absolute calibration of 
a remote sensing system fall into ths three 
cacecerles shewn in the bottom right of Pig. 3t 

(1) The absolute csllbratlon cf the system is 
made onlv before launch. In flight the calibration 
is cheeked by irradiating the focal clane with a 


radiomeerlcallv calibrated source and ootical 
system. The drawbacks to this orocedure are that 
any chance in the transmission of the imace-foraiinc 
ooties of the sensor system, due to the 
condensation of outoassed contaminants, will be 
undetected and the on-board calibration svstem is 
also assumed to be stable through launch and 
unaffected by t.be vacuum, high enerov particle 
irradiation, and zero-c environment at orbital 
altitudes. The Thematic Mapoer and the 
Multispectral Scanner Svstem on C.andsat-C are 
exasDles of remote sensing systems calibrated in 
this manner. 

(2) The sun or an on-board calibrated source 
can be used to irradiate the focal plane through 
the image-forming optics. Ths drawbacks to this 
approach are the uncertaintv in the knowledge of 
(a) the irradianee of the sun above the aemosphece 
and (b) the output of ths calibrated source system, 
for ths reasona mentioned earlier. Furthermore, in 
examples of ths use of this proeadure (MSSa 1, 2 
and 3 and SPOT) , the calibration beam passes 
through only a small nertlon of. the aeertura of the 
syatam, thus not simulating the actual operation of 
the system, when imaging the ground, the system 
entrance aperture is irradiated over its entire 
area by flux incident over a roughly three 
steradlan solid ancle. In the imaging mode there 
is much mere stray light presant in the system and 
incident on the fecal plane. If this additional 
flux lavel la unknown. It may Introduce a 
substantial uncertainty into the absolute 
calibration of the system. 

(c) Reference can be mads in-flight to a 
ground area of known spectral radiance. If at the 
tiro the sensor syeta is imaging the knoMi area, 
measurements are made of the atmospheric eendlticns, 
these data can be used with an atmospheric 
radiative transfer program to predict ths spectral 
radiance at the entrance pupil of the sensor. The 
main uncertainty in this approach is that of 
datanslning the atmospheric aerosol concent well 
enough. The aporoach is also 11ml Cad to scenes 
having large uniform areas of high radiance. For 
exasBle, although many water bodies are of 
sufficient size and unlfotnlty, they ace net 
appropriate for ealibratien purposes because their 
radiance Is too low to provide a calibcacion of 
sufficient accuracy or to cover much of ths dynamic 
ranee of ths sensor. Foctunacsly, some suitable 
areas do exist, particularly in the acid regions of 
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Fig. 3. Claasificatlon of radiamacrle ra Libra tion proeeducea. 





'.:iv ‘njij.i-, -uL <3Aauiwi.e ac nixce aanna> aew .nexxcu, 
in Che Cnite<i States. 

TTie rest st this paoer is ^evotetf to 
discussion oif the factor” and in-fliaht calibration 
of an MLA system us i no the self-calibrated 
ohotodiode approach and the use of a ground 
reference area for calibration ourposes. 


4. CAI.IHRAT10H IS THE FACTORY AND IS ORBIT 

The concent oroposed for the factory 
calibration is similar to the oroposed orbital 
procedure, the main differenee being that an 
artificial source is used in the factory and the 
sun Is used in orbit—siaply a matter of 
convenience in the former cane and of con ve n ience 
and reliability in the latter cane. In the 
factory, redundanev is not at a presUua and our 
reoui resents for a source are siamiy power, 
soeetral flatness, and stability. He do not need a 
standard source although an arrav of standard KBS 
FSL tungsten halogen lasipa could be used, if their 
polariration characterlstiea can be tolerated (Ref. 
Ill . A xenon arc selected for miniasns are wander 
and with a highly stable oower supply and a 
fsedbaet loop would suffice. 

The source would be used to irradiate a 
near-Laisbertian, near-uni ty^refleetanee, 
white-surfaced panel perhaes 1 m x O.S o la sire in 
front of the system. <An Integrating sehere could 
be used, but it would have to be very large, and 
uniformity cheeks can sometlaee themselves 
introduce ncn-unifotaltlea. ) A self-calibrating 
KBS-styie radiometer, with the incident beam 
perhaps defined by two or three apertures, and 
using spectral bandpass filters matching those used 
in the MIA, would be used to determine the radiance 
of the panel' in eseti band. The MLA would iaege the 
panel out of focus, but being an estatided object, 
its inege would have exactly the same Irxadlanee la 
or out of foens. The arranuement is sketched In 
Pic. 4. 



Flo. 4. The factory procedure for absolute 
soeetrotadlometxie calibration. 

To avoid oroblems due to the nonuniform 
irradiation of the oanel, the MLA should be rotated 
to sequentially ifsadiata the focal plane with the 
image of the same small area that is sampled by the 
radiometer. The reason for a large panel is to 
simulate the viewing conditions from space In 
which, depending on the baffle design, significant 
out-of-the-fleld-of-vlew stray lloht could be 
incident on the image olane to modify the 
calibration. For this reason it would be 
worthwhile to conduct at least one calibration 
using a white panel several meters in diameter. To 


c.necx ror Linearity or response, several aitterent 
irradlance levels on the panel should be used. 

This irradlance level can be changed conveniently 
by insertino heat-resistant neutral densitv filters 
in front of the stabilized xenon are source. 

The profiles of the spectral filters should be 
measured in a speetroohetometer usino the same 
F/Ho. beam as the MIA. If thev ate integral with 
the array, they should be measured using a double 
aonochromatoc , aoain with the MIA F/No. , before 
installation in the focal plane. Care should be 
taken to cover the whole wavelength sensitivity 
eanaa of the detectors, the off-band suppression 
being particularly important for detectors with the 
wide spectral response of silicon. 

The procedure proposed hare for in-flight 
calibration Is similar to the oanel method sketched 
in Flo. 4, but it uses the sun as the source. He 
believe that the irradlance over the panel can then 
be considered to be uniform and known soectrally to 
batter than 1% absolute. (Several solar 
measurement programs are currsntlv being conducted 
with this accuracy as a goal. Sowever. if the 
uncertainty la thought to be greater than It, a 
pyroelectric detector could be used to measure the 
direct solar flux in orbit, over the wavelength 
intervals of intaraat and at the sanw tima that the 
system la baing calibrated.) The calibration would 
ba carried out in the faw ninntas while the 
spacecraft is sun-iUuadnatad but before It images 
the sun-illuminated earth. The absolute C 2 idlafflatar 
containing the pyroelaetrie dataetor would now be 
naadad only to ebaek for any datarioration in the 
rsflaetanea of the panel owing to expesnra to tha 
space envirensMnt and short expoaurea to 
unattanuatad tJ7 and other high energy radiation 
from tha sun and from space. In this last respect, 
tha panel would usually ba s to wed in a 
wall-shialdad cempartmant and exposed only during 
actual calibration chaeks. Also, whan daployed, it 
would net Intsrfare with normal operation of the 
system, as it wold ba visaed by the stsreo mirror 
in one of its extreme pesieions. In this respect 
it is fail-safe. 


S. H8ZT8 BANDS AS A CAUBRATION tU gg RBNCB 

To keep the deeription of the theory of the 
method brief, spectral depandaneies are not 
included in the following dlsensalon. In all cases 
tha spectral daoandanea is Implied, tha spectral 
value having tha same wavelength deoendenee as the 
spectral response of each band of tha sensor being 
calibrated. The theoretical basis for the 
calibration method is straiohtforward. 

Tha rsdlanea at tha sensor, L,, is determined 
with respect to the ground radiance, the 
reduction of this radiance by the upward path 
through the atmosphere, '^'^xt ^ (where '^'ext 

the aBBO s p h erie extinction and 0 is the nadir scan 
angle), and the atmospheric path radiance 

(where ^ is the aslsmth angle of the observation) . 
Thus, 

ts ■ i'oekBCf'axt aeoB) +• 

The path radlanee term, in the equation 

is unkROwn and can be determine onlv by the use of 
a radiative tranafer calculation. Such 
calculations regnire a knowledge of ground 

reflectance, the irradlance Incident at the top of 
the atmosphere, the aerosol optical deeth and phase 
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function, tfia oocicai deoth due to ozone and water 
'aoor, tile solar zenith angle, the nadir anola o£ 
the sensor, and its azismthal anole with reaoeet to 
the sun's direction. These radiative tranafar 
calculations assune an infinite, flat, Lafithertian 
reflecting surface with an atstesshece comoosad of 
plane hoategeneeus layers. Althouoh these 
assumntions are never exactly aet, thav ate closely 
aooroximated hv the neariv Lasmertian flat surface 
of the gvosum at White Sands and the snail field of 
view of the ssace sensors. 

The auantity r' can be detemined from the 
solar radiometer dau' tn the following way. The 
solar radiometer records the Irradlanee at the 
around as a function of the solar zenith anale 
during the morning of the overpass. If the 
atnoseherle conditions ars stsady durlno that 
period, the Lanslev plot of In against seed is 
a straight 11ns of slope r* The extrapolation 

of the line to a sees^ value of zero gives the 
value of the irradlanee above the atatosphete in the 
measureaimt pasaband. The value of the around 
radiance Lg is detscmined at the time of the 
overflight by the solar radiometer looking at the 
ground at the sasM S and 4 anglas as that of the 
space sensor. The ground radiance value can be 
converted into the raflaetanee eeguiced as input to 
the radiative tranafar program in one of two ways. 
First, the global irradlanee can be measurad by the 
salat cadiomatar with either a horizontal cosine 
diffuser or an integrating sphere, with horizontal 
entrance port, placed over the entrance aperture of- 
the solar radiometer. Or second, the solar 
radiometer can be pointed down at a horizontal, 
diffusing, wbits surfaea such as aaSOi, or Salon. 

The cosine diffuser aeptoach is preferred because 
its small surface can be c o n ve n iently stored to 
minimize contamination. Intereemparisona will be 
made, during the six monthly ceealibrations of the 
field eguipment, between the cosine diffuser global 
irradlanee and the (.ambertian whits rafarenee 
surface meesuzemenea of the ground reflaetanee. 

The nadir angle from the satellite to the ground 
observation point, the solar zenith angle, and the 
azimuth of the satellite nadir angle to the sun's 
direction are all readily detacaiinabXe. Such 
goantities as atmospheric oressure, humidity, 
surface and air temperature, etc. , are routinely 
.monitored by the Atmospheric Sciences Laboratory at 
Whits Sands. These data will be formatted with the 
radiation data and tslsnatared to Wallops Island. 
Virginia, anitsd States. 

The difficult guantieiss to detatmlne for 
radiative transfer caleuXatlons are the 
characteristics of the aes oa p h erie aerosols. 
Fortunately, we can make reasonable assumptions 
regarding tbeir values based on prior results. For 
example. 79 aerosol profiles have been measured at 
tfhite Sands (Kef. 14) to establish an atmospheric 
model; the effects of aerosol complex refractive 
Index and size dlstrlhutlon on extinction and 
absorption for ths wavaianoth range O.SS to 10.6 um 
for desert atmospheres are reported in rafersnce 
15; and reference 16 reports that tha imaginary 
refractive index is 0.007 at 0.6 urn and shofa 
little dependence on wavelength over the range 0.3 
to 1.1 Urn. 

We intend to investioata some new teehniguea 
for determining aerosol chaetsristies that ars 
under developaient and that make use of sky . 
oolarization measurements. We also Intend to 
explore methods of monitoring the aerosol size 
distribution from inversion of the spectral optical 


deoth measurements (Ref. 17) and aureole 
measur^ents (Ref. 13) . Some of these 
measurements and analyses, particularly ths 
poiacimetric, fall into the research category. 

They reereaent attaupts to improve t.he accuracy of 
t.he calibration procedure and. because of their 
untested nature, regui're validation. However, 
without them we anticipate an overall accuracy of 
about =34 absolute. This reocesents an ImptoveBent 
over the i5% value claimed in reference 19, and is 
due to our anticipated ti% absolute accuracy of 
calibration of the around instrumentation. 

Ths proposed basic design of the 
instrumantatien is similar to that being designed 
and built at tha aniversity of Arizona for the 
Investioation of solar irradlanee variations over a 
23-year period. It ffiOBloya a oreeision alt-azimuth 
tracking stand, with stapoer motors to drive the 
two axes, so that It can be oointed in almeat any 
direction or be held in alignment with the sun. A 
mlcrootoeessor baaed c o m p u tet system will be used 
to control the motors as well as ths data 
aeauisLeion and proesaaing systes. A conceptual 
view of the instrument Is shown in Fla. S, further 
details of its design will be preamted at the 
symposium. 



Fig. S. Conceptual vlw of speetZQCadlo»> 
star for ground and aumospherlc maasuramenta . 
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Atmospheric effects on radiation reflected from soil and 
vegetation as measured by orbital sensors using 
various scanning directions 

Philip N. Slater and Ray 0. Jackson 


Ground-measured spectral reflectance data for Avondale loam and drought-stressed and unstressed wheat 
were converted into digital counts for spectral bands 5 and 7 of the Landsat Multispectral Scanner System 
(MSS). For dry loank, the differences between ratios of MSS bands 7-5 as determined from space and born 
ground level measurementa ware 2.3% for clear and 5.6% for turbid atmospheric conditions. By contrast, 
for wet loam the differences were 10.4 and 29.5%. We found that atmoepheric conditions may cause a delay 
of from 3 to 7 days in the discrimination between drought-stressed and unstressed wheat. For oblique angle 
obaervadona the atmospheric modification of ground-measured reflectances increased with angle at a great- 
er rate in the 0/180" azimuth than in the 90/270" azimuth. Implications of this result are discussed for 
oblique angle Systime Probatoire d’Observation de la Terre (SPOT), Mapsat, future multispectral linear 
array system imagery, and wide-angle imagery collected from scanners in high-altitude aircraft 


I. Introduction 

During the first years that Landsat satellites were in 
orbit there was a paucity of ground-based information 
on which to base interpretations of satellite imagery. In 
the past few years a considerable amount of ground- 
based data has accumulated that describes the spectral 
reflectance characteristics of earth surface features; a 
recent colloquium in fact was devoted largely to 
ground-based measurements of the spectral signatures 
of soil and vegetation.^ The reason for using ground- 
based data is that they can describe the experimental 
results of different agronomic treatments on numerous 
small plots. Such experiments can be conveniently 
conducted repeatedly and rapidly to yield relationships 
between agronomic variables and remotely sensed pa- 
rameters. 

To permit comparison of ground-based information 
and satellite data, atmospheric scattering and absorp- 
tion must be accounted for. Atmospheric scattering has 
been considered by Chance,* Dave,*’"* Otterman and 
Fraser,® Potter,® Richardson et al. ,* Rochon et al.? and 
Turner et al. ® The effect on the Landsat Multispectral 
Scanner System (MSS) band 7 of atmospheric water 
vapor absorption at 0.96 fim has been investigated by 
Pitts et al. (Absorption in the thermal IR bands has 
been studied by sevei^ investigators but is not relevant 
here.) 


Philip Slater is with University of Arizona, Optical Sciences Center, 
Tucson, Arizona 85721; Ray Jackson is with U.S. Department of Ag- 
riculture, Agriculture Research Service, Water Conservation Labo- 
ratory, Phoenix. Arizona 85040. 
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A recent report by Schnetzler^^ concerning the effect 
of the atmosphere on the radiance at the sensor for 
oblique viewing angles is of particular interest He 
determined the conations that allow off-nadir across- 
track angles to be used so that the radiance difference 
between nadir and off-nadir is <5%. His variables in- 
cluded three atmospheric conditions, two wavelengths, 
a range of latitudes &om 60” north to 60” south, solstice 
and equinox dates, two equatorial crossing times, and 
bidirectional reflectance data for a grass target.^***® 
Our results complement those of Schnetzler, although 
fewer variables are considered. Our purpose is to show, 
through use of ground-based data and an atmospheric 
radiative transfer model, how satellite data vary with 
atmospheric conditions. Our data cover a wide range 
of Lambertian reflectances for bare soil and for 
drought-stressed and well-watered wheat. Employing 
the commonly used vegetation index MSS7/MSS5, we 
illustrate the magnitude of atmospheric effects on the 
soil line (a fundamental parameter for the physically 
based model of Kauth and Thomas^'*) and on stressed 
and nonstressed wheat. The data were calculated in 
terms of Landsat digital counts and include an example 
of signal quantization. Two optical depths were con- 
sidered, and the ranges of nadir and azimuthal look 
angles were chosen to include the Landsat MSS, Sys- 
tdme Probatoire d’Observation de la Terre (SPOT), 
Mapsat, and future multispectral linear array systems 
as well as high-altitude aircraft scanners. The specific 
reflectance data used were the soil line for Avondale 
loam described by Jackson et al. and the discrimina- 
tion of stressed from unstressed vegetation described 
by Jackson and Pinter.'® 
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II. Geometric and Atmospheric Considerations 

Figure 1 is a plot of solar azimuth angle (measured 
eastward from north) as a function of the variation of 
solar declination angle throughout the year. The lati- 
tude values for the four cxirves are 55, 30, 15° N, and 0. 
The longitude values corresponding to the position of 
a satellite in a 99° inclination orbit are, respectively, 94, 
104, 108, and 111° W. Actually this particular orbit is 
Landsat path 33 in the Worldwide ^ference System 
index map of the U.S A compiled by the U.S. Geological 
Survey and included in the Landsat Data Users 
Handbook. The orbital path crosses the U.SA. from 
eastern North Dakota to western Texas with an equa- 
torial crossing time of ~0930 h. Because of the nature 
of the Landsat orbit, this particular orbit is typical of 
all others covering this latitude range at the times of 
year shown in Fig. 1. 

The nadir-viewing mode of operation of the proposed 
Mapsat system, Colvocoresses,^^ and Itek,^^ is geo- 
metricallyldenticai to that of the Landsat MSS. 

A. Solar Zenith Angle 

We used the single solar zenith angle of 45° to mini- 
mize the number of costly radiative transfer calculations 
required. This resulted in introduction of a variation 
in the time of year with latitude. Thus by reference to 
Fig. 1, we see that at 55° N the corresponding times of 
year are late April and mid-August; similarly, 30° N 
corresponds to mid-March and late September and 15° 
N to mid-February and late October. 

B. Solar Azimuth Angle 

The solar azimuth angles, measured eastward from 
north, corresponding to a solar zenith angle of 45° at 
latitudes of 55, 30, and 15° N, are roughly 150, 132, and 
128°. Taking into account the roughly 99° orbit^ in- 
clination of the Landsat MSS, the azimuthal angles 
between the orbital path and the sim’s direction become 
141, 123, and 119°; for convenience we chose a single 


value of 120° (see Fig. 2). The MSS scans across a 
swath perpendicular to the orbital plane. Thus, for the 
case considered here, the scan dir^ion is at a 30° azi- 
muthal angle to that direction. 

C. Nadir Viewing Angie 

A 5° angle to nadir was selected. For reference, the 
mavimiiTn angles to nadir for the Landsat MSS and 
Thematic Mapper are 5,8 and 7.5°, respectively. The 
MSS scan direction is not perpendicular to the sun- 
target-sensor plane, so we can expect to obtain different 
anpilar distributions of the radiance levels at the sensor 
on either side of nadir. This is because of nadir angle 
variations in atmospheric path radiance, depending on 
whether the nadir angle is 5° to the east or the west of 
the nadir point. We shall discriminate between these 
cases by referring to 5° E or 5° W, recalling that these 
scan positions are at 30° in azimuth to the sun-target- 
sensor plane. 




Fig. 2. Irradiance and sensing geometry 
selected. Thesolarzenithangleis45'’. The 
scan directions are at azimuth angles with 
respect to the sun-target-sensor plane of (6 
“ 150 and 330®. For convenience, these di- 
rections are designated as westerly (W) and 
easterly (E), respectively. Thus d/v = 5® E 
refers to an angle of view from nadir of 5® in 
a direction 330® in azimuth (measured in a 
clockwise sense) from the sun’s direction. 


Fig. 1. Plot of solar azimuth angle vs solar decli- 
nation angle for a typical Landsat orbit. The day 
and month of the year and the solar zenith angle are 
marked for 5® increments of declination angle. The 
latitudes considered are 0, 15, 30. and 55® N. 


3924 APPLIED OPTICS / Vol. 21, No. 21 / 1 November 1982 




Fig. 3. Role of the atmoephere in attenuating and scattering incident 
and reflected radiant flux.^ 


0. Atmospheric Considerations 

The various ways in which the atmosphere influences 
remotely sensed data collected from high altitudes can 
be appreciated by reference to Fig. 3. Solar radiant flux 
is partially absorbed or scatter^ as it passes through 
the atmosphere to the earth’s surface. The surface at 
B is, therefore, irradiated by direct solar flux and by 
scattered flux from the hemisphere of sky above B. An 
example of this latter contribution is the scattered flux 
from G. The ground scene reflects part of the flux in- 
cident on it in the direction of the remote sensing sys- 
tem. This reflected flux, passing through the atmo- 
sphere, is again absorbed and scattered as shown at F, 
but to it is added scattered flux from the atmosphere 
shown at D that has not been reflected from the ground 
scene. The radiance contribution at the entrance pupil 
of the system owing to scattering at points such as D 
provides a major contribution to what is referred to as 
path radiance. Scattering at point E in Fig. 3 also 
contributes to path radiance and is accounted for in the 
radiative transfer calculations. However, we do not 
consider the case of A having a different reflectance 
than B. 

The primary quantity that determines the influence 
of the atmosphere on the total radiance L, at the en- 
trance pupil of an orbital sensor is the atmospheric 
spectral extinction coefficient or optical depth Ti,,t(X), 
which is deHned in terms of the spectral transmittance 
r(X) along a slant path at a zenith angle d through the 
entire atmosphere by 

r(X) = exp(-ri„(A) secfl). (1) 

A change in Te„(X) causes a number of interrelated 
changes in the ground irradiance and in L,. Thus with 
larger r'e„(\), we get smaller direct solar irradiance at 
the ground but a larger component of diffuse downward 
scattered fliix, greater attenuation of the flux reflected 
from the ground, and greater path radiance. The 
greater path radiance is important in quantitative 
studies of imagery collected from high-altitude aircraft 
or space sensors, because, for low-ground reflectances, 
short wavelengths, or turbid atmospheres, it is the 
predominant component of radiance at the entrance 
pupil of the sensor. 


The radiative transfer calcidation technique devel- 
oped by Herman and Browning, used to obtain the 
results for a simulated atmosphere, and described in the 
following, assumes a flat ear^ with an infinite uniform 
Lambertian surface and a horizontally homogeneous 
atmosphere. All orders of atmospheric scattering are 
included. (Polarization effects could have been cal- 
culated, but they were not for this analysis.) The at- 
mospheric parameters used as input to the model as- 
sumed a Junge size distribution of nonabsorbing par- 
ticles of refractive index 1.54, which correspond^ to a 
Mie extinction coefficient that varied in inverse pro- 
portion to the wavelength. The output from the radi- 
ative transfer calculations gave the total radiance 
emerging from the top of the atmosphere at 5, 15, 25, 35, 
and 45° to nadir as a function of five ground reflectances 
(0, 0.1, 0.25, 0.5, and 0.75) for two wavelengths roughly 
at the centers of MSS5 and MSS7 (0;£5 and 0.9 fim), for 
four azimuthal planes (0/180, 30/210, 60/240, and 90/ 
270°), and for a clear and a turbid atmosphere. 

Figure 4 is a plot of spectral extinction coefficient 
against wavelength for the clear and turbid atmospheres 
examined and, for reference, for a Rayleigh atmosphere. 
The extinction coefficients as a function of wavelength 
are listed in Table I for the Rayleigh and aerosol scat- 
tering components. The meteorological range for the 
clear atmosphere is ~100 km, which is typical of a clear 



Fig. 4. Atmospheric extinction coefficient as a function of wave- 
length for various atmospheric conditions. 


TaMe L RayMgh, NDe, and Total Extinction Coeffidants for the OUtarant 
Atm oa p h araa Studied 


Atmosphere 

Wavelength 

(pm) 

Rayleigh 

Mie 

Total 

Clear 

0.65 

0.048 

0.026 

0.074 

Turbid 

0.65 

0.048 

0.397 

0.445 

Clear 

0.90 

0.013 

0.020 

0.033 

Turbid 

0.90 

0.013 

0.300 

0.313 


1 November 1982 / Vol. 21. No. 21 / APPUED OPTICS 


3925 


day in the desert Southwest of the U.SJL and many 
other arid regions of the world. The meteorological 
range for the turbid atmosphere is 10 km. Flowers et 
al. ^ have found that, for selected urban sites across the 
U.S A., on 55% of the days in the year the meteorological 
range is <10 km; for suburban and rural areas the per- 
centages are 18 and 1%, respectively. 

Several assumptions and approximations have been 
made to simplify the analysis: 

(1) Lambertian reflectance has been assumed 
throughout. 

(2) A flat earth has been assumed which means that 
results for large angles to nadir predict too small an 
atmospheric effect. The actual slant range (h in me- 
ters) can be calculated approximately from 

A ■ (/? costf — (S* co8*9 — 2HR sin*S)*^/8iii*tf, (2) 

where R - earth radius in meters, 

H a sensor altitude in meters, and 
9 s slant angle to nadir. 

For high-altitude aircraft the difference between the 
curved and flat earth slant ranges is insignificant, even 
when d a 450 , However, for a nadir angle of 25” at an 
altitude of 700 km, used in a later example, the effective 
atmospheric path for a curved earth corresponds to a 
nadir angle of 23.7°. 

(3) Square passbands between 0.6 and 0.7 and 0.8 
and LI pim, having the atmospheric characteristics of 
monochromatic wavelengths at 0.65 and 0.9 Atm, have 
been assumed to model the Landsat MSS sensor. In 
view of the accuracy of the radiative transfer calculation 
method used, this approximation may seem unneces- 
sarily rough. However, at present, we cannot charac- 
terize the atmosphere well enough from measurements 
on a global scale to merit using more accurate simula- 
tions. For example, even for one particular remote 
sensing measurement, it is difficult to determine 
whether absorbing aerosols are present and, if they are, 
what complex indices of refraction to ascribe to ^em. 
Furthermore, the adjacency effect described by Pearce^ 
has not been experimentally validated and is difficult 
to account for even for simple scenes. 

(4) The effect of atmospheric water vapor was not 
considered, although it must be in all practical appli- 
cations of MSS7 data. The reasons for neglecting water 
vapor absorption were twofold: first, it can be largely 
accounted for from the results of Pitts et who 
showed that the water absorption band centered at 0.94 
Atm can reduce the radiance at the sensor by a factor of 
0.77; second, later systems such as the Thematic Map- 
per and SPOT use near-IR passbands that avoid strong 
water absorption bands. 

The following calculations are designed, then, to in- 
dicate representative trends rather than predict accu- 
rately the influence of specific atmospheric character- 
istics on the radiance at the sensor. We start by de- 
scribing the procedure for determining the relationship 
between the digital counts for MSS5 and MSS7 for the 
two atmospheric conditions defined and for Avondale 
loam and wheat (Triticum durum Desf. veir. Pro- 
dura). 



Fig. 5. Radiance at the sensor as a function of Lambertian ground 
reflectance for dear and turbid atmoapheric conditions in bands 5 and 
7. Curves are for unity irradiance, normal to the sun’s direction at 
the top of the earth’s atmosphere, and a solar zenith angle (9,) =• 

45*. 


III. Calculation Procedure for Avondale Loam 

(1) Figure 5 was plotted fiom the results of the 
radiative transfer mt^el to relate the radiance at the 
sensor L, to the ground reflectance p. Note that the 
radiance values correspond to an incident irradiance of 
unity at the top of the atmosphere. Figure 5 is not at 
a sc^e to show differences between the 5° E and 5° W 
curves, and hence only the four curves for the 5° E case 
are shown; however, the least-squares fits for all eight 
cases are listed in Table II. 

(2) Using the best fit equation in Fig. 6 for Avondale 
loam, values for pj were determined corresponding to 
P5 values of 0, 0.1, 0.25, 0.5, and 0.75. 

(3) The radiances L5 and L7 corresponding to the sets 
of ps and p^ values were determined by substituting the 
appropriate reflectance values in Eqs. (3) and (7) listed 
in Table IL (Later, in plotting the relationship between 
Ls and L^ for low p values, all eight equations listed in 
Table II were used.) 

(4) A least-squares equation was determined that 
related Ls and L^. 

(5) Steps (3) and (4) were repeated using Eqs. (5) and 
(9) to yield a least-squares fit for the relationship be- 
tween Ls and L7 for the more turbid atmosphere. The 
least-squares fits for the two relationships are 

L^ = 2.6058 X 10-3 + i.2363Ls - 0.17032L1 (11) 

for the clear atmosphere, 

Lt = -2.4725 X 10-3 + i.2815Z,5 - 0.17685L1 (12) 

for the turbid atmosphere. 

(6) Four sets of corresponding Ls and L^ values from 
the equations in Table II were converted to digital 
counts over the low reflectance range 0 < p < 0.1. 
These were used to determine the change in radiance 
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TaM a. LMit-SquarM Fits (or tiM RacHanea at tha Samer bi Tama at Qraond Raflaetaneas 


Scan azimuth 
wj.t. sun* 


X (^m) 

T«t(^) 

(deg) 

Eq. 

a 

b 

c 

0.65 

0.074 

30(E) 

(3) 

4.545 X 10-® 

0.2106 

1.194 X 10"* 



30 (W) 

(4) 

4.972 X 10-® 

0.2108 

1.155 X 10-® 


0.445 

30(E) 

(5) 

1.202 X 10-* 

0.185 

3.343 X 10-2 



30 (W) 

(6) 

1.270 X 10-* 

0.185 

3.343 X 10-2 

0.90 

0.033 

30(E) 

(7) 

1.393 X 10-® 

0.2203 

4.641 X 10-2 



30 (W) 

(8) 

1.538 X 10-® 

0.2202 

4.645 X 10-2 


0.313 

30(E) 

(9) 

6.659 X 10-® 

0.1996 

2.434 X 10-2 



30 (W) 

(10) 

8389 X 10-® 

0.1999 

2.395 X 10-2 


° a, b, and c are coefficients in the equation L, >* a + bp +• epH nadir angle is 5*. 

* (E) and (W) refer to easterly and westerly directiona; that is, the 5* nadir angle is pointing into the direction of the sun and at 30* to the 
-:un-target-sensor plane in the (E) case and away from the direction of the sun at 30* to the sun-target-sensor plane in the (W) case. 



5 Count* 


Fig. 7. Plots of digital counts in band 7 vs band 5 for Avondale loam under clear and turbid atmospheric conditions. 



<60 (OAHO SI SEFLECTAHCe 

i* ig- 6. Band 7 vs band 5 reflectance for Avondale loam, Phoenix, 
Ariz. The coefficient of determination (R-SQR) was 0.983. Data 
■vere taken between 1330 tuid 1400 h from Nov. 1978 to May 1979.'® 


at the sensor owing to a change in the viewing angle from 
5° E to 5° W and to examine the effect of quantization. 
Results are discussed later and presented in Fig. 3. 

The procedure for converting the radiances in Eqs. 
(11) and (12) to digital counts from the Landsat MSS5 
and MSS7 was as follows. First, the actual radiances 
at the entrance pupil of the MSS in bands 5 and 7 were 
found by multiplying the normalized radiance values 
by 15.1 and 25.1. These are the solar irradiances in mW 
cm~^ in bands 5 and 7 as determined hrom solar spectral 
irradiance data tabulated by Thekaekara et al. 2“* Then 
a range of pairs of MSS5 and MSS7 radiance values was 
converted to digital counts (DC) by use of the following 
relationships^: 

Band 5 DCs = (Is - 0.07)/0.0U7; (13) 

Band7DCs = (L7-0.14)/0.0637. (14) 
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Fig. 3. Digital counts in band 7 vs band 5 for dear and turbid at* 
mospheres for scan positions of ±5” to nadir. The staircases show 
the elfect.of Landsat MSS 6-bit quantization. 


IV. Results in Terms of Digital Counts 

Figure 7 is a plot of the digital counts in band 7 vs the 
corresponding counts in band 5 for the two atmospheric 
conditions. Note that the substitution of 127 and 63 
counts in Eqs. (13) and (14) yields the levels of entrance 
pupil radiances (1.56 and 4.15 mW cm~^) that corre- 
spond to the onset of detector saturation for the 
Landsat 2 MSS during the 22 Jan. to 15 July 1975 pe- 
riod. 

The digital counts are the consequence, of course, of 
the quantization of the output signals from the MSS. 
However, Fig. 7 does not reflect this condition in two 
respects. First, the straight lines should have a stair- 
case pattern, each step corresponding to one or more 
counts. Second, band 5, like bands 4 and 6, is initially 
quantized into 64 steps and later during ground pro- 
cessing expanded to 128 steps. If the initial quantiza- 
tion is linear, not logarithmic, each of the 128 steps has 
a precision of one part in 64, not one part in 128. If the 
initial quantization is logarithmic, the precision will 
vary throughout the dynamic range. The noise Equiv- 
alent radiance difference detectable by the MSS is ap- 
proximately half of a digital count in each channel, so 
with a probability of 95% the radiance error for a single 
pixel is ±1 count or ±1.6%. 

Figvire 8 shows the plots of digital counts in bands 5 
and 7 for clear and turbid atmospheric conditions and 
for a nadir look angle of 5” in both the easterly and 
westerly azimuthal directions. Because of the coarse 
quantization of the MSS, shown by the uneven staircase 
patterns, we caimot differentiate in most cases between 
the 5® E and 5® W results. However, the quantization 
resolution is quite adequate when related to the stan- 
dard deviation of data, presented in Fig. 6, for the re- 
flectances of Avondale loam. In this regard it is inter- 
esting to note that, when imaging identical ground 
features, the future s}^tems with their greater quanti- 
zation resolution of 256 steps will output significantly 
different digital counts between easterly and westerly 
scan positions. This is particularly true for systems 
that operate in a non-nadir viewing mode, as we shall 
see later. 



Tta* (Ju)i«n 0«v») 

Fig. 9. Effect of the atmosphere on the MSS digital counts from 
stressed O and unstressed -f vegetation, (a) Results of ground- 
measured reflectance ratios (umnodifled by the atmoephere); (b), (c) 
results using simulated Landsat counts through clear and turbid at- 
mospheric conditions, respectively. 



HMir '/itwing Angit. 7^ (degrc««) 

Fig. 10. Normalized values of L«,,/Lo* vs nadir viewing angle for \ 
= 0.9 pm, p = 0.0, and 4> = 0/180, 30/210, 60/240, 90/270“ for clear and 
turbid atmospheres. 
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V. Results for Vegetation 

Jackson and Pinter^^ have shown that the ratio of 
band 7 to band 5 reflectances can be used to discrimi- 
nate water-stressed vegetation. Results of radiative 
transfer calculations, shown in Fig. 5 and listed in Table 
II, were applied to the measured values of band 7 and 
band 5 values reported in Ref. 16 following the proce- 
dure outlined earlier for Avondale loam, llie resulting 
curves for the band 7 to band 5 ratios, for high quanti- 
zation resolution digital counts, are shown in Figs. 9(b) 
and (c) for clear and turbid conditiona. For compari- 
son, Fig. 9(a) is a plot of the ground reflectance ratios 
for stressed and unstressed vegetation. 

Calculations show that the 6-hit quantization of the 
MSS is sufficient to differentiate between the ratio 
values in Fig. 9(a) at Julian day 94 at the 95% proba- 
bility level, but it is insufficient for the same day under 
clear and turbid atmospheric conditions. The 8-bit 
quantization of future sensors will differentiate between 
the ratios in Fig. 9(a) over the entire period of interest 
from day 77 to day 118. It will also provide differen- 
tiation of the ratios at about day 88 under clear condi- 
tions [Fig. 9(b)] and at about day 94 under turbid at- 
mospheric conditions [Fig. 9(c)]. 



Fi«. 11. Ratio of digital counts in band 7 to band 5 for the Landsat 
MSS as a function of the reflectance in band 5 for Avondale loam. 
Results are for three atmospheric conditions: no atmosphere and a 
clear and a turbid atmosphere. The nadir view angles are 25” E, 25” 
W. 45” E, and 45” W. 


VL Oblique Angle Imagery 

The effect of the atmosphere can be illustrated most 
dramatically by considering the radiance at the entrance 
pupil of a space sensor that is imaging a region of zero 
reflectance. This radiance is then entirely path radi- 
ance. Figure 10 is a plot of the radiance at the entrance 
pupil of the high altitude or space sensor as a function 
of the angle between the scan direction and the nadir 
(which we call the nadir angle) to that for zero nadir 
angle. Several nadir angles 9 and azimuths, 0 were 
considered. (The 0 = 0/180° azimuth describes the 
vertical plane containing the sun and the target with 0 
= 0 representing the dire^on into the sun.) The curves 
are for a wavelength of 0.9 nm and for the clear and 
turbid atmospheric conditions described earlier. 

The following conclusions can be drawn from Fig. 10. 
First, the largest variation in the radiance at the en- 
trance pupil of the sensor occurs for the 0 = 0/180° az- 
imuth, and the least occurs for the symmetrical case of 
the 0 = 90/270° azimuth. Second, the variation for the 
asymmetrical cases is greatest for the direction away 
fiom the sun. Third, the variation is greater for a given 
azimuth the more turbid the atmosphere. 

The shapes of the families of corresponding curves 
for X = 0.65 fim and/or for p > 0 are similar to those in 
Fig. 10; however, the variations are not as pronounced. 
For example, the uppermost curve in Fig. 10 has maxi- 
mum values for Lgti/Loo of 2.1 in the direction into the 
sun and 3.2 in the direction away from the sun. The 
corresponding ratios for X = 0.65 pm and p = 0 are 1.8 
and 2.7, and those for X = 0.9 pm and p = 0.1 are 1.2 and 
1.5, respectively. These general conclusions are ap- 
parent in the following specific examples. 
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Fig. 12. Normalized values of vs iiadir viewing angle for 

Pxaes “ 0-0 0.5, ij> = 150/330, and 60/240” for (a) clear 

and (b) turbid atmospheres. 
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Pv,, = 0.5, Ptoi* “ 0 

pvu“0.5,pw -0.1 

Atmospheric 

Azimuth ansias 

Azimuth anides (dex) 

condition 

330/150 60/240 

330/150 60/240 

Clear 

Turbid 

+13 -25 +7 -16 

-8 -26 -6 -15 

+2 -5 +1 -3 

-3 -10 -2 -5 

T^lelV. Pot 

remaBS milei ew.ee DeteieOT at Nedh and at 

48“ 


Pv.. ■ 0.5, Px«« ■ 0 

Pv^, - 0.5. px-« - 0.1 

Atmospheric 

Azimuth amties (des) 

Azimuth anzies (deg) 

condition 

0/180 90/270 

0/180 90/S70 

Clear 

Turbid 

-7 -51 -17 -17 
-42 -59 -28 -28 

0 -15 -3 -3 

-19 -33 -11 -11 


Figure 11 is a plot of the ratio of the digital counts in 
band 7 to band 5 as a function of the reflectance in band 
5 for Avondale loam. The results are for no atmosphere 
and for a clear and turbid atmosphere. Only one azi- 
muth (Fig. 2) is considered; the angles to nadir are 25 
and 45°. The 25** case corresponds to an operating 
mode of Centre National d*Etudes Spatiales’ SPOT or 
NASA’s future multispectral linear array (MLA) sys- 
tem, in which imagery can be acquired by pointing the 
scanner to either aide of nadir, orthogonally to the flight 
path. (In some cases a 35^ angle is being considered to 
increase the frequency of coverage.) The 45* case 
corresponds to a high-altitude aircraft system, such as 
the Daedalus scanner, whose extreme scan angle is 
nearly 45*. 

The cvirves in Fig. 11 show how the atmosphere in- 
creasingly modifles the ratio of MSS7 to MSS5 digital 
counts at the sensor as (1) the ground reflectance de- 
creases, (2) the angle to nadir increases, and (3) the at- 
mosphere becomes more trubid. For wet loam with a 
band 5 reflectance of 0.1, the band 7 to band 5 ratio, 
unmodified by the atmosphere, is 0.47. The ratio drops 
to 0.425 at 25® E and 0.39 at 25® W for a clear atmo- 
sphere and to 0.365 and 0.34 for a turbid atmosphere. 
For 45® E and 45® W, the ratios become 0.415 and 0.37 
for the clear and 0.33 and 0.315 for the turbid atmo- 
sphere. It is interesting to note that, as the values of the 
ratio become further removed from the unmodified 
value, the differences between the values across the scan 
line become smaller. With reference to Fig. 10, we see 
that for a clear atmosphere and with respect to a zero 
nadir angle a change in scan angle of 25® to the east 
produces a 3.4% change in radiance at the sensor, and 
one of 25® to the west produces a 5.1% change. These 
changes are significant when compared with the ~0.5% 
noise equivalent radiance differences of future sensor 
systems. 

Figure 12 relates for clear and turbid atmospheric 
conditions, respectively, the normalized ratios of radi- 
ances at X = 0.9 fiva to those at \ = 0.65 nm as a function 
of nadir angle for two azimuths <i> - 150/330® and <t> - 


60/240* with respect to the sun. (Normalized radiance 
values were chosen to facilitate the extension of the 
calculations to the bandpasses at or near 0.65 and 0.9 
Atm to be used by some fiiture systems.) The 0 =• 60/ 
240* azimuth is representative of the plane of the fore 
and aft viewing directions for in-track stereo scanners 
such as the proposed Mapsat system. The 0 ^ 150/ 
330* azimuth is representative of sensors that achieve 
stereo and/or increased temporal coverage by pointing 
off-nadir in the across-track direction. Exmnples are 
SPOT and MLA. 

In Fig. 12 the flatter lines are for ground reflectances 
of 0.1 in band 5 and 0.5 in band 7; the more curved lines 
are for ground reflectances of zero in band 5 and again 
0.5 in band 7. The region between curves of the same 
azimuthal angle encompasses many representative 
vegetation ratios for radiwces at 0.9 and 0.65 nm. The 
values plotted in Fig. 12 are summarized in Table m. 
From these results it is clear that, from a rSdiometric 
standpoint, it is preferable to obtain oblique imagery 
from a Landsat-type orbit by pointing in the forward 
in-track direction. Results from Table IV show how 
desirable it is to collect wide-angle high-altitude aircraft 
imagery when the flight path is in the sun’s direction 
rather than when the path is perpendicular to that di- 
rection. 

VII. Conciusiora 

The results presented here indicate the rauige over 
which the atmosphere can be expected to modify band 
ratio results (MSS7 to MSS5 or Lj to Ls); the atmo- 
spheres considered had meteorological ranges of 100 and 
10 km. In all cases considered, the modification was 
significant primarily because, for the examples chosen, 
the reflectance at 0.65 fiia was lower than at 0.9 )um, and 
the effect of the atmosphere was, therefore, more pro- 
nounced at 0.65 Axm. 

In the case of Avondale loam imaged by the MSS, the 
influence of soil moisture is critical. For dry loam, the 
difference between the band 7 to beind 5 ratios as de- 
termined from space and as determined from ground- 
level measurements is 2.3% for clear and 5.6% for turbid 
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atmospheric conditions. In contrast, for wet loam the 
differences are 10.4 and 29.5%. Because of the low 
quantization resolution of the MSS and its small 
(<5.6°) angle to nadir, the difference between the ratios 
in the east and west scan directions is hardly detectable 
even for very low band 5 reflectances, which is indicative 
of damp loam. This is not the case for a near-nadir 
high-quantization resolution sensor such as the Ther- 
matic Mapper, and the effects become much more 
pronounced for oblique viewing systems. For example, 
the band 7 to band 5 ratio, for Avondale loam of re- 
flectance 0.1, determined from space measurements at 
a 25° nadir angle, is less than the vertical ground-based 
measurement by 10% in the east direction and 17% in 
the west direction for a clear atmosphere and by 23% in 
the east direction and 28% in the west direidion for a 
ttirbid atmosphere. Compared to the values of the ra- 
tios determined from space measurements at a nadir 
viewing, angle, the increases are —4.4, -(-3.9, 0, and 
-1-6.1%, respectively. The effect of the atmosphere on 
band ratios for vegetation indicates that atmospheric 
conditions may delay the discrimination of stressed 
from unstressed vegetation by S-7 days. 

In general the results indicate that the atmosphere 
modifies band ratios significantly for low-ground re- 
flectance values. The degree of the modification in- 
creases for oblique angle observations, and the path 
radiance and the change in path radiance are greater in 
the 0/180“ azimuth than in the 90/270“ azimuth. The 
results also indicate that if in-track stereo data aure to 
be used for radiometric purposes, they should be col- 
lected in the fore and nadir directions rather than the 
aft and nadir directions. The more similar radiometric 
responses between the fore and nadir directions may 
facilitate the roughly 1-D digital correlation of the two 
images for stereo purposes. 

Results emphasize the need for development of 
techniques for measurement of atmospheric parameters 
from space to facilitate the correction of spectral sig- 
natures for atmospheric effects. This is particularly 
true if in the future comparisons are to be made among 
near-nadir, oblique, and ground-based measure- 
ments. 
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Abstract 


Several methods for the in-flight calibration of advanced reaoce sensing systems are reviewed. Pre-flight 
and in-flight detector-based calibration are discussed. Calibration by reference to Instrumented ground 
targets is suggested for current systems. Preliminary data are presented on the accuracy and senaicivicy 
required of measurements of ground cacgees and the atmosphere between Che targets and the sensor to achieve 
an uncertainty in absolute terms of less chan z3 Z in the calibration of Landsat-4 systems. 

Introduction 


Many remote sensing applications and studies, particularly chose chat depend on Che modeling of the 
interaction of radiant flujc with the earth's surface and/or atmosphere, require an accurate conversion of 
the sensor's output (in digital counts) to the radiance level at the entrance pupil of the sensor. The 
accuracy of this conversion depends on the suicabilicy and accuracy of the absolute speccroradlomecrlc 
calibration of the system. Unfortunately, none of the in-flight absolute calibration procedures used for 
remote sensing systems have been entirely appropriate. They have involved directing a small beam of solar 
radiation through the optics, while leaving the rest of the entrance pupil dark. When imaging the earth, the 

entrance pupil is irradiated with a uniform radiance level over its entire area from a field of roughly 3 

steradiana. this inevitably introduces a level of stray radiation across the image surface that is not 

accounted for in the calibration procedures used thus far. (The stray radiation is due to scattering off 
optical surfaces, supporting structures, etc., and can increase dramatically with small amounts of 

contamination chat may accumulate on them in the space environment.) 

Although chrough-che-optlcs radiometric calibration techniques have been available on the early Landsat 
Hultlspectral Scanner Systems (HSSs), they have not been successful and have not been used to any great 
extent. For the HSS on Landsat-4 the calibration wedge approach^ will be used to determine "absoluca" and 
relative band-to-band response. Scene content hiscogramming will be used to monitor detector-to-detaetor 
gain and bias to achieve the zi quantum level specification. These techniques have worked well for earlier 
HSSs, but we should remember t^t hiscogramming is a relative correction technique, suitable only for 
removing scene-independent scrlping,^'^*^ and the calibration wedge provides only a focal plane calibration. 
It provides no information on Che cransmiccance of the image-forming system (including scanning mirror). 
Because the condensation of outgasaed contaminants and the presence of high-energy irradiation in the space 
environment can rapidly change the reflectances of mirrors and Che cransmltcances of refracting elements, a 
chrough-che-optlcs procedure must be used for the absolute calibration of the space remote sensing system to 
be meaningful, there is no provision for a through-the-opcics calibration of the Thematic Happer (TH) on 
Landsac-4. A focal plane calibration only will be used, in conjunction with a pre-flight absolute 
calibration. 

The Systhme Frobacoire d'Observacion de la Terre, Haute Sasolutlon Visible (SPOT H&V) utilizes a fiber 
optics arrangement that transfers solar flux from an external surface to a relay system. The main Imaging 
system together with the relay system then focuses an image of the internal ends of the fibers onto the CCS 
detectors." The absolute radiometric calibration of SPOT H&V is quoted^ as 10 Z. Specifications for advanced 
multlspectral linear array systems call for uncartaincles in absolute calibration of zSZ,^*^ while one 
detailed design scudy^^ has estimated that the uncertainty can be as low as tl.4Z. 

In this paper we present a brief review of relative and absolute calibration methods for advanced remote 
sensing systems, pointing out anticipated problems and limitations. We briefly describe the technique for the 
self -calibration of a photodiode in absolute terms using several laser wavelengths and how this has been 
extended to calibration in terms of a tungsten source filtered with narrow, well-characterized interference 
filters. We suggest this detector-based method for the pre-flight and in-flight calibration of future 
systems. For current systems, where there is no on-board provision for absolute calibration, we discuss 
absolute calibration with reference to large uniform ground areas Instrumenced to provide ground spectral 
reflectance and atmospheric data. 

Broad classification of calibration procedures 


Some of the most commonly used procedures for the calibration of remote sensing systems are referred to 
In Figure 1. The two major divisions in the figure are between relative and absolute calibration and between 
the static macro-image response and the dynamic micro-image response of the system. Relative calibration 
and dynamic micro-image response will cot be discussed further, beyond remarking chat the dynamic micro- 
image response is of vital interest in any pixel-by-pixsl analysis of remotely collected imagery and will be 
studied by several investigators as part of the Landsac-4 Image data quality analysis program. The rest of 
this paper deals with the absolute calibration of remote sensing systems, and only the static macro-image 
response will be considered in this context, as is usually the case. 
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Flgtira U Classificadoo of radiomacilc calibracloa procedures. 


The procedures for Che absoluce calibracloa of a rsBOCe seasing syscen fall Inco Che chree cacegories 
shown In che boccom rlghe of Figure 1: 

(O The absoluce calibracloa of che sysceo la nade only before launch. In fUghc, che callbracion la 
checked by Irradlaclng che focal plane wlch a radlooecrically callbraced source and opclcal syscea. The 
drawbacks eo ~ chls procedure are chac any change in che cransnlssion of che loage-fomlng opcics of che 
sensor syscem, due co che condensacloo of oucgassed coocaalnancs, will be undececced and che on-board 
calibracloa syscea Is also aaauoed co be scable chrough launch and unaffecced by che vacuum, high energy 
parclcle Irradlaclon, and zero-g envlronoenc ac orblcal alcicudes. The Tbeaaclc Mapper and che Mulcispeccral 
Scanner Syscea on Landsac-4 are examples of reaoce sensing sysceas callbraced in chls oanner. 

(2) The sun or an on-board callbraced source can be used co Irradlacs che focal plane chrough che laage- 
forming opcics. The drawbacks co chls approach are che uncercaincy in che knowledge of (a) che Irradlance of 
che sun above che acaosphere and (b) che oucpuc of che callbraced source syscea, for che reasons oencloned 
earlier. Furchermore, in examples of che use of chls procedure (MSSs 1, 2 and 3 and SPOT), che calibracloa 
beam passes chrough only a saoU porclon of Che apercure of che syscea, chus noc slaulaclng che accual 
operaclon of che sysceo. Uhen iaaglng che ground, che syscea encrance apercure Is Irradlaeed over Ics entire 
area by flux incident over a roughly chree sceradlan solid angle. In che iaaglng node there is ouch more 
scray Ughc present in che sysceo and incident on che focal plane. If chls additional flux level is unknown, 
ic introduces a substantial uncertainty into che absoluce callbracion of che system. 

( 31 Reference can be made in-f Ughc co a ground area of known spectral radiance. If oeasuremencs are 
made of che acaospheilc conditions ac che time che sensor syscea is iaaglng che known area, these data can 
be used with an acmospherlc radlaclve transfer program co predict che spectral radiance ac che entrance 
pupil of the sensor. The main uncercaincy in chls approach is chac of determining che atmospheric aerosol 
concent well enough. The approach is also Halted Co scenes having large uniform areas of high radiance. 
For exaaple, although many wacer bodies are of sufficient size and uniformity, they are noc appropriate for 
callbracion purposes because chelr radiance is coo low co provide a calibration of sufficient accuracy or co 
cover much of che dynamic range of che sensor. Fortunately, some suitable areas do erisc, particularly in 
che arid regions of che world, for example ac White Sands, New HeUco. 

The rest of this paper is devoted co a discussion of che factory and in-fUghc calibration of an HLA 
system using che self-callbraced photodiode approach and che use of a ground reference area for callbracion 
purposes. 


Dececcor-based radiometric callbracion 


The recenc work ac che US National Bureau of Standards on self-callbraced photodiodes is described only 
briefly here. For more details che interested reader is referred co references 12-16. 

The calibration of che phocodlode is accomplished by either of two differently applied biasing procedures, 
depending on che wavelength region of interest. Ac short wavelengths, a negative bias is applied co remove 


















Che recooblnaclon cencers ac che SI-SIO 2 lacerface ac the froac of che dececcor. To do chls, a concacc la 
made »ich che front surface using an eleccrode Immersed In a conducting liquid, or che surface can be exposed 
CO a corona discharge. For long wavelengths, a back-bias is applied co extend che depletion region co a 
depth beyond which incident flux penetrates. The experimental procedure is co irradiate che dececcor with a 
constant monochromatic flux level and co Increase che bias voltage until further increase no longer gives 
rise CO an increase in output signal. For both che short and che long wavelength ranges, che internal 
quantum efficiency saturates ac a value extremely close co unity, as shown in Figure 2. Thus che m ax imum 
increase in signal output obtained as a result of biasing can be used co determine che internal quantum 
efficiency of che dececcor wlchouc biasing, as it will be used in practice. 

The only significant loss in che external quantum efficiency of che photodiode is caused by reflection. 
This can be reduced co insignificance by making use of three photodiodes according co che geometry sketched 
in Figure 3. 
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Figure 2. Typical photodiode internal quantum Figure 3. A chree-diode arrangement co minimize 

efficiency without biasing (lower curve) and specular reflection losses, 

with biasing (upper curve), reference 14. 


The output signals from the three diodes are summed co provide che signal corresponding co a black 
detector of overall quantum efficiency chat can be assumed co be unity. The second diode reflects che 
specular reflection from diode 1 co diode 3, which acts as a recrorefleccor. The incident flux thereby 
undergoes five reflections, and ac a 10 Z value Che final specular reflectance is down to 10*^ of che inicial 
incident radiance. Ic is claimed^ ^ chat che diffuse reflectance losses for clean detectors are typically 

much less chan IZ. The photodiodes used in similar calibration facilities at NBS and ac che Onlversicy of 

Airlzona, Optical Sciences Center, are EGG OV 444g. 

The discussion ac che end of this paper refers co che use of self-callbraclng NBS detectors for che 
speccroradlomecric calibration of an MLA system. However, as shown in Figure 2, che unbiased quantum 

efficiency is wavelength dependent. Because of possible changes in che passband position of spectral filters 
during long duration space flights, ic may be advisable to use che NBS detectors co calibrate a spectrally 
flat pyroelectric detector, ac che 0.1-0.2Z level, and for chat co be used in che in-orblc calibration. 

A recent paper^^ shows chat photodiode self-calibration using filtered tungsten radiation, in place of 
laser radiation, can provide an absolute spectral response scale from 0.4 co 0.S um with an uncertainty of 

less chan cl Z. The filters used in chls work were roughly of 10 am half-width, which is about a factor of 

10 less chan chose of the MSS and TH. The possibility of extending this filter approach or a monochromator 
approach^^ co provide a means to calibrate field and flight speccroradlomecric instrumentation with an 
uncertainty in absolute accuracy of the order of cl Z appears then co be close co reality. 


Calibration in che factory and in orbit 


The concept proposed for che factory calibration is similar co che proposed orbital procedure, che main 
difference being chat an artificial source is used in che factory and che sun is used in orbit — simply a 
matter of convenience in che former case and of convenience and reliability in che latter case. In che 
factory, redundancy is not at a premium and our requirements for a source are simply power, spectral 
flatness, and stability. Ue do not need a standard source although an array of standard NBS FEL tungsten 
halogen lamps could be used. If cheir polarization characteristics can be colerated^*^. A xenon arc selected 
for mini mum arc wander and wich a highly stable power supply and a feedback loop would suffice. 


The source would be used co irradiate a aear-Lambertian, near-unity-refleccance, white-surfaced panel 
perhaps I m x 0.3 m in size in front of che system. (An integrating sphere could be used, but ic would have 
CO be very large, and uniformity checks can sometimes themselves introduce aon-uniformlties.) A self- 
callbraclng NBS-scyle radiometer, with che incident beam perhaps defined by two or three apertures, and using 
spectral bandpass fllcers matching chose used in che MLA, would be used co determine che radiance of che 
panel in each band. The MI.A would image che panel out of focus, but being an extended object, its image 
would have exactly che same irradiance in or out of focus. The arrangement is sketched in Figure 4. 
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Fig. 4. The factory procedure for absolute speccroradloaetrlc calibration. 


To avoid problems due co the aonunlform Irradiation of the panel, the KLA should be rotated to 
sequentially Irradiate the focal plane with the Image of the same small area chat la sampled by the 
radiometer. The reason for a large panel Is to simulate the viewing conditions from space In which, 
depending on the baffle design, significant ouc-of-che-field-of-view stray light could be Incident on the 
Image plane to modify the calibration. For this reason It would be worthwhile to conduct at least one 

calibration using a white panel several metres In diameter. To check for Unearity of response, several 
different Irradlance levels on the panel should be used. This Irradlance level can be changed conveniently 
by Inserting heat-resistant neutral density filters In front of the stabilised xenon arc source. 

The profiles of the spectral filters should be measured In a spectrophocomater using the same F-No. beam 
as that for the MIA. If they are Integral with the array, they should be measured using a double 
monochromator, again with the MIA F-Ho., before Inatallatloa in the focal plane. Care should be taken to 
cover the whole wavelength sensitivity range of the detectors, the off-baad suppression being particularly 
Important for detectors with the wide spectral response of silicon. 

The procedure proposed here for In-fUghc calibration Is similar to the panel method sketched In Figure 
4, but It uses the sun as the source. We believe that the Irradlance over the panel can then be considered 
to be uniform and known spectrally to better than 1 X absolute. (Several solar measurement programs are 
currently being conducted with this accuracy as a goal. However, If the uncertainty Is thought to be greater 
than IZ, a pyroelectric detector could be used co measure the direct solar flux In orbit, over the 
wavelength Intervals of Interest and at the same time chat the system la being calibrated.) The calibration 
would be carried out In Che few minutes while the spacecraft is sun-Ulumlnatsd but before It images the 
sun-illuminated earth. The absolute radiometer containing the pyroelectric detector would now be needed 
only CO check for any deterioration In the reflectance of the panel owing co exposure co the space 
environment and short exposures co unaccenuaced OV and ocher high energy radiation from the sun and from 
space. In this last respect, the panel would usually be stowed In a well-shielded compartment and exposed 
only during actual calibration checks. Also, when deployed. It would not Interfere with normal operation of 
Che system, as It would be viewed by the stereo mirror la one of its extreme positions. In this respect It Is 
fall -safe. 


White S ands as a calibration target 

White Sands, New Me:dco, has been proposed as a calibration sice for the Landsat-4 sensors. This ares 
provides a large uniform target, the reflectance of the gypsum sands Is close co Lambertian for the solar 
angles of Interest, and the concentration of acmospheric aerosols under calm wind conditions Is low. These 
properties simplify radiative transfer computations. To compute the radiance at the sensor, a sec of ground- 
based measurements is combined with a proven radiative transfer model. 

The sensor radiance, Lj, Is given by 

Lg - E • exp(-f„c S«c8s^ ’ 7 S 

(spectral quantities assumed). Here 9g Is the sensor down look angle (as measured from nadir), o the ground 
reflectance, and E the solar spectral Irradlance. A Langley plot of In E versus sec bg (where is the 
solar zenith angle) can be drawn from data taken throughout the morning of the Landsac-4 overflight. The 
slope of this plot is a measure of the total atmospheric extinction (t'gxc)> ^nd the exo-acmospheric spectral 
solar Irradlance (denoted Eq) Is found by extrapolating co the value of E where secb^ • 0. Path radiance, Lp, 
is computed from data inputs such as Eg and p and from assumptions about the composition of the atmosphere 
at Che time of the Landsat overflight. 








Seoslcivlcy Analysis 


Previous efforcs co esclmace radiance levels ac doun-looliLlag orbical sensors have been, ac best, ac che 
rS2 uncercaincy level, alchough some comparaclve studies have shovn char much larger uncercaincles can 
exisc. For example, che raclo of oucpucs of che S192 co che S191 enerlmencs on Skylab was 0.79 - 0.14 when 
Che syscems were slmulcaneously viewing che same ground scene.^^ Uich improved inscrumencacion and 
calibracion accuracies (co rlZ), Ic now seems feasible co use ground-based reference cargecs for sensor 
calibration. A senslcivlcy analysis is required, however, co predict Just what this new accuracy might be. 
This Involves a study of che parameters affecting che radiative transfer codes, and chelr effect on che 
predicted radiance. Such work has been initiated, and preliminary results are presented here. Currently the 
radiative transfer codes agree among themselves co within 1 Z for a given set of input parameters.^^ 
Problems arise when specific values of these inputs have co be defined so as co represent a given atmosphere. 
Typically, approximations and seasonal averages are used for parameters relating co aerosols (number density, 
refractive index, radius distribution, and vertical discribucion), ozone, and water vapor concentrations, and 
molecular gases. Often, however, these parameters vary considerably from chair seasonal average values. 
This la especially true of aerosol and ozone and water vapor concentration. Fortunately, large variations in 
some of these parameters (10 Z) often lead co tolerable errors (<IZ) in calculations of radiance ac che 
sensor. Table 1 lists chose variables chat affect che transfer calculations. Comments are included co 
summarize cheir expected variability and our preliminary study into the' effect each has on overall accuracy, 
tc is seen chat ocher parameters such as Che ground reflectance, dme of day, and solar spectral irradiance 
should be known co 0.1 Z. This is not an unreasonable requirement, given che recent advances in calibracion 
techniques. Our simulation was made by using an average desert aerosol number density, with complex index of 
n~l.S-0.01i. Published daca^^ were used for che ozone optical depth and absorption coefficients. Dp co now 
water vapor variability, which affects che fourth MSS band, has not been analyzed. 


Table 1. Preliminary Examinacion of Input Parameters co Radiative 
Transfer Codes 


Aerosols 

Total number Density and Vertical Distribution. Soch varlable,^5 aerosol optical 

thickness varying, by as much as a factor of 2. 

Particle Size Discribucion . Normalized lunge and Deirmenjisn distributions compared. 
Number of particles ac a given radiiu found to vary by an order of magnitude. 
Change in output was <0.1 Z. 

Refractive Index. Extinction may vary by an order of magnitude for che various 
values used co approidzisce desert aerosols,^^ extinction being more affected by che 
real part of che refractive index. Our calculations show chat a change in index from 
n-l.S CO 1.5-O.Oli varies output by <1Z. 

Other. Paraaecers 

Ozone. Can expect up co 10 Z daily variation from seasonal mean,^^ this changes 
output by 1 Z . 

Water vapor . With humidity range of 30Z co 90 Z, and cemperacure of 4°C 'to Zl'C, 
water concentration varies from 0.13 to 3 gm cm~^ km'^. 

Ground reflectance, exo-acmoapheric solar spectral Irradiance. For negligible oucpuc 
error, need co know chese parameters co within 0.1 Z. With recent detector 
calibracion techniques, this is no longer unreasonable. 

Field of view . A study has been done using an Infinitesimal target having che 
reflectance of White Sands (p~0.7) and a background of p~0.33 and 0.7. The resultant 
output varied by Z Z . 

Pressure. A 10 Z change in pressure changes output by only 0.2 Z. Negligible error 
anticipated as pressure can be measured on a continuous basis co high accuracy. 

Time of Day . Need co know wichin 1 second of measurement for Langley ploc.^^<^^ 

Polarization. Scudies^^’-^^ need co be extended co determine che influence of scene 
and atmosphere Induced polarization on che calibration of an orbical sensor. 
Polarization introduced by che sensor itself needs co be known from pre-flighc 
measurements. 
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Conclusions 

The use of self-callbraced dececcors promises Che achievement of pre-flight and on-board in-flight 
absolute calibration of future systems to an uncertainty of less chan 

Uichouc appropriate pre-flight and on-board calibration, it now appears possible chat present earth 
resources satellite systems can be radiometrlcally calibrated in absolute terms with an uncertainty of less 
chan C3Z by reference to grotind targets. This conclasion is based on recent advances in the establishment 
of self-calibrated detectors as accurate absolute radiometric standards together with the development of 
radiative transfer codes chat account for multiple scattering and include all relevant atmospheric 
parameters. The degree of uncertainty In the calibration procedure depends heavily on our ability to 
accurately characterize the atmosphere above the target during the sensor overflight. We plan to investigate 
this matter both theoretically and experimentally during the image data quality analysis phase of the 
Landsac-4 program. 
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Abstract 


The importance of accurate absolute radlooetric calibration is discussed by reference to the needs of those 
'^riahing to validate or use models describing the interaction of electromagnetic radiation with the atmosphere 
and earth surface features. The In-flight calibration methods used for the Landsat Thematic Mapper (TM) and 
the Systeme Probatoire d'Observatlon de la Terre, Haute Resolution Visible (SPOT/HHV) systems are • described and 
their limitations discussed. The questionable stability of in-flight absolute calibration methods suggests the 
use of a radiative transfer program to predict the apparent radiance, at the entrance pupil of the sensor, of a 
ground site of measured reflectance Imaged through a well charactered atmos^ere. The uncertainties of such 
a method are discussed. 


Introduction 


The use of satellite-acquired image data in scene models is inextricably dependent on atmospheric correction 
and sensor absolute radiometric calibration. Unlike automated scene classification, which depends on a statis- 
tical analysis of the digital counts in a scene image, usually without correction for the intervening atmos- 
phere, modeling is concerned with determining the radiance of the scene. For this purpose we need to know, 
first, the output digital counts Oram the sensor when it is imaging the feature of interest, second, the abso- 
lute radiometric calibration of the sensor in ordo* to convert the digital counts to radiance at the entrance 
pupil of the sensor, and third, the radiance modificatlai introduced by the intervening atmosphere (the atmos- 
pheric correction) in order to relate the entrance pupil radiance to the radiance of the ground feature. 

The problem divides into three parts: (1) determining to what accuracy we have to know the feature radiance 
in ordo- to produce satisfactory data fhom the various models available— not a well understood relationship in 
many cases, (2) provl<Ujig a means to make atmospheric corrections to a certainty commensurate with (1). above, 
preferably using the imaging sensor or a co-located system— even simple ground-based atmospheric measurements 
are rarely attempted today, and (3) providing in-flight absolute radiometric callbratlcn to a certainty commen- 
surate with (1) and (2)— the estimated &10S uncertainties of current systems are unlikely to meet these needs. 

This paper discusses the need for accurate absolute radiometric calibration in the use of scene radiation 
models. It then describes calibration methods presently in use and the main limitations and sources of error 
associated with them. Finally, .a method that utilizes simultaneous ground and atmospheric measuremoits is sug- 
gested for the in-flight calibration of the newer sensors. 

Use of absolute calibration in scene models 


The main reason for developing models is to obtain quantitative and/or diagnostic information regarding spe- 
cific areas or features or to study global 'phenomena. Examples of the. former are to diagnose the cause of a 
loss of crop vigor or to determine and map chlorophyll concontration in prospective fisheries. Examples of the 
latter are global earth-atmosphere-ocean studies such as those envisioned by COSPAR ' 3 International Satellite 
Land Surface Climatology Project and those proposed as part of NASA's Global Habitability program. 

A good source of information on the various models developed over the past decade can be found in a recent 
review by Smith. ^ Most of the models predict an upward radiance, just above the feature, on the basis of cer- 
tain input values. For example, in the case of a crop cOTopy, the input values would include the leaf spectral 
reflectance and canopy geometry (leaf area index, LAI, and leaf orientation), the soil spectral reflectance, and 
the geometry of illumination and viewing. For many purposes the inverse form of this model is more valuable: 
Given the radiance and illumination and viewing geometry, what is the constitution' (LAI, vigor, etc.) of the 
canopy? Coel and Strebel^ have described such a model; however, relatively little other work on inversion 
models has been reported. 

The sensitivity of both direct and inversion models to measurement accuracies and assumptions needs much 
furtha* exploration. In this respect it is interesting to note how various scientists working with models re- 
spond informally to the question '"^/hat sensor absolute radiometric calibration uncertainty can you tolerate in 
using your models?" The answers should represent the model's sensitivity. Some responses are as follows: 

0 . 5 % D. G. Goodenough (Canada Centre for Remote Sensing) to provide the correction for change in 
sensor response with time, necessary for multi-temporal studies. 

1 5 J. L. Barker (Coddard Space Flight Center) in support of his layered concept. For example, 
the removal of atmospheric effects to obtain BBDFs followed by the use of the BHDF data tp 
interpret subtle texture changes and natural variations. 
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IS J. Cower (Univa’sity of British Columbia) for ocean color studies, In particular chlorophyll 
concentration determination. 

3S J. J. Bunnik (NLB, the Netherlands) for vegetation and ocean model studies. 

3S J. Dozlo' (UC Santa Barbara) Car snowfield model studies. 

3S J. A. Smith (Colorado State Univs'sity) for general model studies and plant canopy models in 

particular. 

5S B. 0. Jackson (Agricultural Research Service, Phoenix) for evapotranspiratlon models. 

These responses deserve more detailed study. In particular it is important to determine, fa* each response: 
(1) whether absolute radiometric calibration meant a calibration in physical units (radiance) or a stability of 
relative calibration (in digital counts) with time, (2) whether it was assumed that atmospheric effects had been 
perfectly corrected, and (3) whether these values were well corroborated and a study had been made to deter- 
mine informatioD loss as a function of calibration uncertainty. . 

The literature describing the use of models for the aialysls of satellite-acquired digital image data is 
meager. The works of Aranuvachapun and LeBlond,^ Ooak et al.,'* Gordon,^ Kowalik et al.,^ Otterman and Fraser,^ 
and Roblnove et al.^ represent perhaps the moat significant contributions. The lack of reported results can be 
related to (1) the difficulty of making accurate atmospheric corrections, (2) the large uncertainty in sensor 
absolute radiometric calibration, sui (3) the fact that not many models have been developed and their experi- 
mental testing and exploltatian involves considerable effort. 

In-flight calibration 

We describe in-flight radiometric calibration by reference to the methods used on the SPOT/HRVs and the Th. 
The limitations of both methods are ^so discussed . 

Calibration of the HRVs 

On the HRVs, two calibration procedures®*^® are to be used. One pro- 
vides absolute calibration for a few detectors in each of the four 
Image plane arrays. The other provides relative callbratian for all 
the detectors in the arrays. 


Both the absolute and relative methods provide calibration of the 
complete optical system. The radiation sources for both calibration 
methods, the internal ends of optical fibers and a lamp, are located in 
a calibration unit positioned so as to relay radiant flux through the 
telescope optics to the detector arrays. 

The absolute calibration, to be conducted once each month in orbit, 
is made by reference to the solar Irradlance incident on bundles of 
fibers 3 m long whose external ends are perpendicular or nearly ?«•- 
pendioular to the sun's rays and whose internal ends are located in the 
calibration unit. To cover the range of solar incidence angles Cram 
winter to summer, three seta of fibers are used with gradient refrac- 
tive index tips manufactured by Selfoc. One fiber in each set emerges 
from a different port on an external surface at a different aigle to 
the sun. The Selfoc tips provide an output constant to 1% for a solar 
angle variation of The fibers are drawn from doped silica, and 

tests show they do not discolor in a simulated high-energy-particle 
space environment. 

Because of the possibility of a relative displacement occurring 
between the calibration unit and the arrays, particularly during launch, 
and an uncertainty in the angular position of the pointing mirror, the 
images of the fibers are arranged to follow a staggered line along the 
arrays. The arrangement of the fibers in the focal plane of the cali- 
bration unit is shown in Figure 1. The calibration unit and HRV produce 

a magnification of about 4.3X, so the diameter of a 250-uffl-core fiber, 

when aligned with the array, covers more than 50 detectors in their 
13 -tim along- line dimension. The image of the staggered fiber array 
covers a lateral displacement of about 2.'* mm in the HRV focal plane. 

Figure 2 is a schematic of the calibration unit. It shows the inter- 
nal ends of the fibers imaged, afta* reflection by a beamsplitter cube, 

onto the focal plane of the HRV as described above. A Philipa Osram 

20-W, 12-V helical filament tungsten halogwi lamp is used to uniformly 
irradiate the aperture of lens C. By use of Kohler illumination, this 


Fibers are about 3SO urn in diameter 
with cores about 2SO um in diameter. 



Muitispectral Panchromatic 

band band 

Figure 1. Staggered images of 
the riba’s in the focal 
plane of the SPOT/HRVs. 
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Figure 2. S^ematic of the SPOT/HRV calibration unit. 


aperture is Imaged onto the focal plane 
of the HRV, Thua the entire focal plane 
L3 uniformly irradiated with white 

light. This arrangement can then be » - — J 

used to calibrate the system in a rela- . 

tive sense and thereby transfer the yj\ ^ ^ 7 

absolute calibration of the detectors — — — 

calibrated by use of the fibers to ail / j ' ' '"^V- — ~ 

the detectors in each array. L^p 

The accuracy of the absolute call- '\ / V ” 

bratlon method is limited by (1) the ^ W SEAMSPUTreR 

uncertainty In our knowledge of the \j 

exo-atmospherlc solar spectral irradi- I 

ance. (2) the accuracy to which the I 

spectral transmittance of the 3-m-long n 

fibers and the optics in the calibration ' n\ 

unit can be measured, (3) the constancy //L^® FIBERS EACH 

of the transmission of the fiber optics /I \ 

when exposed to the space environment, / I \ 

particularly including any contaminatian ^ 

of the external fiber surfaces, aid 
(4) the amount of stray light present in 

the system in its normal operating Figure 2. Schematic of the SPOT/HRV calibration unit, 

mode. We shall return to this last 
point later. 

Calibration of the TM 

The calibration of all the detectors of the TM 'is carried out by an opto-mechanical system referred to, 
because of its action, as a "flag." The flag oscillates back and forth at the same ihequency as the scanning 
mirTor. tXirlng the 10.7 ms of scan mirror reversal time, the flag blocks the scene radiant flux from the de- 
tectors and provides various levels of irradiance for the calibration. 

Figure 3 shows the internal calibra- 
tor flag. It consists basically of 

three light sources whose radiant . _ _ . 

output is transmitted by optical fibers 
to a prism end lens arrangement that 
projects radiant flux uniformly onto the 
96 detectors in bands 1 to 5 and 7. In 
addition, a concave a mirrar focuses the 
output from a temperature-controlled 
blackbody source onto the four detec- 
tors in band 6. To minimize inertia, orShunarRig 

the fiber optica and prism-lens arraige- (■«< Shuiwr T«no» 

ment are on a beryllium support that is 
suspended fhcm three flex pivots and 
driven at resonant frequency by a 
brushless DC motor. The flag is essei- 
tially a torsional spring mass system 
that does not use bearings and gears 
and therefore needs no lubrication. 


or Stainer Rig 
laid Staeiar Tempt 



(and Hub Tempi 


The illumination system for the in- 
ternal calibrator is shown in Figure 4. 

The three stationary, miniature tungsten 

filament lamps are at 120o intervals 3Lpnoe.fitare 

about the flag rotation axis. Prisms ^ Foidfeig (Space 

are used to arrange the images of the 
three filaments into a delta pattern, as 
seen in Figure 5a. Attenuating masks, 

which follow the prisms shown in Figure Figure 3. Internal calibrator for the TM. 

5b, modify the outputs of the lamps so 
that they are in the ratio 2:3:4 on 

on entering the circular flba* bundle. The lower end of the fiber bundle is uniformly irradiated by the radiant 
flux from each filament. Thua, as the lamps are turned on in different combinations, irradiance levels in the 
ratio 0:2:3:4:5:6:7:9 can be obtained. 


Figure 3. Internal calibrator for the TM. 


The circular fiber bundle, of diameter 3-3 mm, is divided into six smaller circular bundles, each of diameter 
1.35 mm, and each coupled across a short air gap to another bundle, of 1 x 1.5 mm rectangular cross section. 
Each bundle provides flux to a different TM band. The air gaps can be adjusted to provide small radiometric 
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THREE STATIONARY 
illuminators 



Figure ‘t. Uluffllnator optica of the TM internal calibrator. 


475-04 




Figure 5. TM internal calibrator 
lamp attenuators. 


corrections in any of the bundles. The top end of each fiber bundle is attached to a solid light pipe, which 
provides a uniform radiance area to be projected, by the lenses that follow, onto the detector plane, totween 
the first lens and prism in each band there is an aperture stop whose adjustment provides one way for matching 
the calibrator output to the selected saturation level for the system. The upper end of the flag also contains 
a black surface of known temperature that is viewed by all the detectors during the DC restoration step. 

Behind each source is a spectral filter and _a, silicon PIN detector to ounitor the source output and thereby 
control the current input to maintain constant "source radiance. This scheme counteracts the reduction in radi- 
ance that results from the buUd-up of tungsten deposited on the inside of the lamp envelope during the life of 
the lamp. The spectral filter transmittance was chosen to be in the range 0.65 to 0.66 \m, for which SI detec- 
tors have essentially zero temperature coefficient of responsivity (i0.01S/oC). 

Linatationa of the methods 

The SPOT/HRV calibration method provides a means for both absolute and relative calibration through the op- 
tical system. The condensatlan of outgassed materials on the external ends of the fibers and the effect of 
direct high energy solar radiation on these materials could change the response of the absolute calibration 
system. The relative calibration system does not include a feedback system to ensure stabilized source output. 
Because the tungsten filament becomes fragile after use, the calibration source will be turned on only briefly 
for essential calibration purposes before the HRV is in orbit. To minimize the risk of lamp failure, a sugges- 
tion was made to assume stability of response of the CCDs from laboratory calibration through launch and orbit 
insertion and to calibrate the source, at the earliest opportunity in orbit, against the known CCD response. 

Perhaps the most serious drawback of the absolute calibration method is that it does not simulate the condi- 
tions of image data collection. Only a 100-mm-diameter area of the 33(^mfflHjiameter entrance pupil receives 
flux from the calibration unit,^ and the solid angle subteided by the source lens in the calibration unit at the 
entrance pupil of the HRV is small. In contrast, when the HRV is viewing the earth, the entire entrance pupil 
is irradiated and the solid angle subtended by the surface of the earth at the HRV is about ir steradians. 
Although the actual field of view of the HRV is much smaller than this (1.3 x 10”3 sr), it sr is a measure of 
the total angular range of the irradiance on a horizontal plane in front of the HRV. Radiant flux ova* most of 
this IT sr solid angle Is not going to be imaged directly by the HRV, but some of it is going to be reflected and 
scattered off the front optical elements, their supports, and the barrel of the camera and onto the detectors. 
Such stray light is typically, a few percent in a clean, well designed optical system. However, the percentage 
can rise dramatically if surfaces become dusty or if a coating of scattering material 1s deposited on them. 
Under imaging conditions like this the transmission of the system may actually appear to increase although in 
reality the signal of interest has decreased a little and this has been more than compensated for by an in- 
crease in the DC stray light background level. Thus, the absolute calibration system on the HRVs is probably 
not sufficiently sensitive to the stray light problem. 

The main limitation of the TM calibration method is that it checks the calibration of only the filters, de- 
tectors, and electronics although it also includes the relay optics for bands 5, 6, and 7. Thus it does not 
provide, a measure of the change in the transmittance of the image-forming optics and the change in reflectance 
of the scan mirror and the scan line correction mirrors. 


The calibration layout illustrated in 
Figure 6 shows that the port in the in- 
tegrating sphere Is about the same size 
as the entrance port to the TM. ftaring 
calibration they are placed about 5 m 
apart. This arrangemeit does not simu- 
late the orbital imaging conditions of 
the TM, and it can introduce a system- 
atic error into the TM calibration. 

T^ type of OTor was discussed in 
reference to the SPOT/HRV. 

The internal calibration system con- 
tains long glass paths aid more than 20 
glass-vacuum Interfaces. To the extent 
that these are not affected by the con- 
densation of outgassed materials and 
high energy flux in the space environ- 
ment, the intomal calibrator should 
prove a stable referoice because of its 
use of stabilized sources. It should 
then be useful as an accurate reference 
for relative calibration purposes. 

Because of these potential shortcomir 
stable than the systems they are intended to calibrate, both CNES and NASA are exploring the use of a ground 
refffl’ence area for in-flight absolute calibration.^^ 

White Sands as a calibration target 

White Sands, New Mexico, is being used as a calibration site for the Landsat-A/5 sensors because this area 
provides a large uniform target, the reflectance of the gypsum sands is close to Lambertian for the solar 
angles of interest, and the concentration of atmospheric aerosols under calm wind conditions is low. These 
properties simplify radiative transfer computations. To compute the radiance at the sensor, a set of ground- 
based measurements is combined with a proven radiative transfer model. 

The sensor radianee, L^, is given by 

L, a E exp(-T'ext sec 9s) -y -r Lp 

(spectral quantities assumed). Here E is the solar spectral irradiance, ^'sxt is the total atmospheric extinc- 
tion, 9s is the sensor downlook angle, (as measured from nadir), p is the ground reflectance, and Lp is the path 
radiance. A Langley plot of In E versus sece^ (where Sz is the solar zenith angle) can be drawn from data 
taken throughout the morning of the Landsat-4/5 overflight. The slope of this plot is a meaure of r'exti snd 
the exo-atmospheric spectral solar irradiance (deioted €q) Is found ^ extrapolating to the value of E, where 
sec8z = 0. Lp is computed from data inputs such as Eo and p and from assumptions about the composition of the 
atmosphere at the time of the Landsat overflight. 

The instrument that is used at White Sands employs a precision alt-azimuth tracking stand, with stepper 
motors to drive the two axes, so that it can be pointed in almost any direction or be held in alignment with 
the sun. A microprocessor-based computer system is used to control the motors as well as the data acquisition 
and processing system. Figure 7 is a diagram of the optical layout of the instrument. Figure 8 is an artist's 
drawing of the instrument and its mounting. The ground-based measurements to be made include ground reflec- 
tance, solar irradiance versus time of day, temperature, and relative humidity. 

Previous efforts to estimate radianee levels at down-looking orbital sensors have been, at best, at the 
uncertainty level, although some comparative studies have shown that much larger uncertainties can exist. 
For example, the ratio of outputs of the 3192 to the 3191- experiments on Sky lab was 0.79 t 0.14 when the sys- 
tems were simultaneously viewing the same ground scene. With improved instrumentation and calibration accu- 
racies (to ill), it now seems feasible to use ground-based reference targets for sensor calibration. A sensi- 
tivity analysis is required, however, to predict just what this new accuracy might be. This involves a study of 
the parameters affecting the radiative transfer codes, and their effect on the predicted radiance. Such work 
has been initiated, and preliminary results are presented here. Currently the radiative transfer codes agree 
among themselves to within 1 S for a given set of input parameters. Problems arise when specific values of 
chese inputs have to be defined so as to represent a given atmosphere. Typically, approximations and seasonal 
averages are used for parameters relating to aerosols (number- density, refractive index, radius distribution, 
and vertical distribution), ozone and water vapor concentrations, and molecular gases. Often, however, these 
parameters vary considerably from their seasonal average values. This is especially 'true of aerosol and ozone 
and water vapor concentration. Fortunately, large variations in some of these parameters (10%) often lead to 
tolerable errors (<1%) in palculations of radiance at the sensor. Table 1 lists those, variables that affect 
_the transfer calculations. Comments are included to summarize their expected variability and our preliminary 
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Figure 6. The radiometric test 
arrangement for the TM. 


and because In-fUght sun calibration systems seem potentially less 
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Figure 3. Artist's concept of instru- 
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Table 1. Preliminary Examination of Input Pa'ameters to Radiative Transfer Codes 


Aerosols ; 

Total Number Density aid Vertical Distribution. Both variable, with aerosol optical thickness 
varying by as much as a factor of 2. 

Particle Size Distribution. Normalized Junge aid Oeinaenjian distributions compared. Number of 
particles at a given radius found to vary by an order of magnitude. Change in output was <0.1 S. 

Refractive Index. Extinction may vary by an order of magnitude for the- various values used to 
approximate desert a«^30ls, extinction being more affected by the real part of the refractive 
Index. Our calculations! show that a change in index from ns 1.5 to n s 1.5-0.011 changes output 
by <1*. 

Other Parameters ; 

Ozone. Can expect up to 10% daily variation from seasonal mean;i^ this changes output by 1%. 

Water Vapor. With humidity range of 30% to 90%, and temperature of 4oc to 210C, water concentra- 
tion varies from 0.15 to'3 gm cnr2 karK 

Ground Reflectance, Exo-Atmospheric Solar Spectral Irradlance. For negligible output error, need 
to know these parameters to within 0.1%. With recent detector calibraticn techniques, this Is no 
longer unreasonable. 

Field of View. A study has been done using ai infinitesimal target having the reflectance of '<Rnte 
Sands (0 s 0.7) and background of p s 0.35 and 0.7. The resultant output varied by 2%. 

Pressur--. A 10% change in pressure changes output by only 0.2%. Negligible error is anticipated, 
as .issure can be measured on a continuous basis to high accuracy. 

Time of Day. Need to know within 1 second of oeasuremant for Langley plot.^®,^® 

Polarization. Studies need to be extended to determine the influaice of scene- and atmosphere- 
induced polarization on the calibration of an orbital sensor. Polarization introduced by the 
sensor itself needs to be known from pre-flight measurements. 


study into the effect' each has on overall accuracy. It is seen that other parameters such as the ground re- 
flectance, time of day, and solar spectral irradlance should be known to 0.1%. This is not an unreasonable 
requirement, given the recent advances in calibration techniques. Our simulation was made by using an average 
desert aerosol number density, with complex index of n = 1.5-0.011. Published data^ w«"e used for the ozone 
optical depth and absorption coefficients. Up to now, water vapor variability, which affects the fourth MSS 
band, has not been analyzed. 

Results of the first attempt to make such a calibration for Landsat <*, together with further discussion of 
the uncertainties involved, can be found in Castle et al.^^ 
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Conclusions 


WLthouC appropriate pre-flight and on-hoard caLlhratlon, it now appears possible that present earth 
resources satellite systems can be radiometrically calibrated In absolute ta'ms with an uncertainty of less 
than ±31 by reference to ground targets. This conclusion is based on recent advances in the establishment of 
seif-calibrated detectors as accurate absolute radiometric standards together with the development of radiative 
transfer codes that account for multiple scattodng and Include all relevant atmospheric parameters. The 
degree of uncertainty in the calibration procedure depetds heavily on our ability to accurately characterize the 
atmosphere above the target during the sensor overfUght. 14e are investigating this matter both theoretically 
and experimentally during the image data quality analysis phase of the Landsat-4/5 program. 
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Abstract 


This paper discusses field radiometer methods for measuring (1) the reflectance factor of a surface, (2) the 
ratio of atmospherically scattered to direct Irradlance (s/d) at the ground, and (3) the atmospheric extinction 
coefficient. 

Calculations show that, under hazy or cloudy conditions, reflectance factor measurements of an unknown sur- 
face made In the field with reference to a white panel, both surfaces having nonlambertlan characteristics, can 
differ by up. to 6Z from laboratory measurements of the unknown surface. This applies to surfaces having re- 
flectance factors greater than 0.05. The error can be reduced to 0.4 Z If the direct solar component alone Is 
used for the determination. 

Measurements of surfaces with reflectance factors from 0.09 to 0.4 showed errors of lOZ and 2Z respec- 
tively when the total radiance of the target was ratloed to chat of Che reference panel. These errors can be 
reduced to 4Z and less than IZ respectively when the direct solar components are ratloed. 

The mld-lnfrared (mld-IR) bands of a commonly used field radiometer showed a high out-of-field response that 
gave rise to measurement errors on Che order of 20 Z. 

The effect of the reflectance of other surfaces in the neighborhood of the target Is demonstrated by deter- 
mining the ratio of shaded to direct Irradlances. Agricultural scenes can show changes of about 5Z In the red 
and 20Z In the near Q>. 

A commonly available field radiometer. In conjunction with a reference panel, can be used reliably to deter- 
mine the atmospheric extinction coefficients In broad wavelength Intervals. 

Introduction 

Measurement of Che reflectance factor of ground features has become a routine procedure by a number of in- 
vestigators. Mulclspeccral radiometers have been designed especially for such measurements.^ Reflectance 
factor data obtained from these Instruments are used as 'ground truth" for aircraft and satelUte experiments 
and are used to verify models chat have been developed to describe the radiation environment of surfaces. 

Such radiometers can also be used for atmospheric measurements. The separation of the atmospheric diffuse 
component from the global component Is usually made by using two pyranomecers,^ one exposed to the total 
global radiation and the ocher to the diffusa radiation. The attenuation of solar radiation by atmospheric 
scattering, usually measured with solar radiometers, can also be evaluated with field radiometers by measuring 
Che radiance of a target of known reflectance factor chat la alcamately sunlit and shaded at a number of 
solar zenith angles. The atmospheric extinction coefficient can be obtained by the Langley plot method. ^ 

It has been pointed out^*^ chat a number of factors should be considered when making field measurements, 
for example, the variation in solar Irradlance with a change In solar zenith angle, the effect of cloud move- 
ments on Che scattered Irradlance during measurements, and the effect on surface Irradlance of che proximity of 
Che observer. The errors Introduced by assuming uniformly distributed skylight and lamberclan reflectance 
characterlsclcs of che carget and reference panel have also been described.^*' 

This paper develops the equations to: U) describe che errors In reflectance factor measurements Introduced 
by these assumptions together with che assumption chat che radiometer has perfect out-of-fleld rejection char- 
acteristics; (2) determine che Influence of Che reflectance of che surrounding area on che measurement of che 
ratio of diffuse Co direct Irradlance, caking Into account che Imperfect ouC-of-fleld rejection characteristics 
of che radiometer; (3) determine che correction coefficient ' needed In che measurement of che ratio of diffuse 
CO direct Irradlance co account for che high radiance In che vicinity of che sun, which Is shielded from che 
panel during che measurement of che diffuse component; and (4) determine atmospheric extinction coefficient 
using a mulclband radiometer and horizontal reference panel. Example results are calculated from these equa- 
tions and compared co che results of field measurements. 


Development of Equatlona 




The Irradlance at a target surface located In a field envlroament la composed of three parts: (1) the direct 
Irradlance from the sun, (2) the component scattered by the atmosphere before Interaction with the ground, and 
(3) the component that comprises flux reflected from the ground and scattered by the atmosphere onto the sur- 
face to be measured. Figure 1 shows the geometry of a radiometer viewing a reference reflectance panel In a 
field environment. 


The scattered Irradlance from the hemispherical sky Is 


®3ky 


- 2 « 


*/2 

IskyC®) Sind coad dd, 

0 


where I.siiy(d} Is the component of solar flux scattered by the at- 
^ mosphere and Incident on the target without previous Interaction 

* with the ground. Then, the scattered Irradlance on the surface can 

be written as 

Edif • Egfey + (Edlrcosdc + Ejicy}P(ae,8a). (1) 

where Ecu.^'cosdj Is the solar Irradlance at the surface at a solar 
zenith angle dj, and PCRc.B^) la a coefficient representing the Irra- 
dlance contribution from the ground-reflected flux chat has been 
V subsequently scaccered downward. F(Bf,Sg) Is a function of the sur- 

round reflectance and the atmospheric conditions at the dme of 
measurement. and S 3 are the reflectance factors of the target 
and the surround, respectively. By definition, a reflectance factor 
Is the ratio of the radiant flux from a surface Co chat from an 
Identically irradiated perfectly lambertlan surface. 



Fig. 1. Measurement geometry with a 
field specCroradlomecer. 


We sec Q • (E^^coadg E 3 icy)P(Rt,Sa}, which Is the irradlance contribution due to Che surround. Hoeing chat 
this can be taken as a uniform Irradlance on the target, and considering Che angular reflection characteristics 
of the target, we can write 


2dlfat(9) 


[Eskyl 


1 + 5e(»)l + 


QjOt(d) 


QR^(8) + 2i 


»/2 

L,l(y(S) S((d) Sind cosd dd, 

0 


( 2 ) 


where the term 1 -*■ fc(9) Is a reflectance factor correction coefficient, which accounts for the nonlamberclan 
characcerlsclcs of the target and the nonuniform hemispheric scaccered Irradlance on the surface. If the tar- 
get is perfectly lambertlan, then Sf(d) • and Cc(9) • 0 . 


Referring to Fig. 1, we see chat Che radiometer voltage signal can be written as 


'^tot.c(9> - [Sdircosd, + Egityll + 5c(9)l + Qj^tC®) r(») &ia 


|E<iirCOS8i + Eskyll + 5s(9)l + q]Es( 9) KdO, 


(3) 


9 where CgCd) is a correction coefficient with the same meaning as except that it applies to the surround, 

r(o) is Che spatial responslvlcy of the radiometer, K Is a constant chat includes Che diameter of Che aperture, 
Che filter factor, and Che transmittance of Che radiometer, and dO is Che solid angle increment, defined as 


dO 


2»rdr 

(H/coso )2 


2w tana do. 


To account for out-of-field-of-view concrlbutlon to the radiometer signal, we define 


<1 


2»K 


a, 

r(o) Cano do 

0 


and 


K2 


2»K 


r(o) tana do 


(4) 
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where and IC 2 are voltages per unit of reflected radia n ce from the target and surround, respectively. They 
can be determined If the spatial responslvlty r(a) of the radiometer and the subtended angle sq are known. 
Substituting Eqs. (4) Into Eq. (3), we have 

'^tot.t^®) “ ^E^lirCoaSi + ^kyl^- * 5e(95! + S-i 

+ [Bdlrcose* + EskyU + «a(9)l + q]h,(9) K 2 . (5) 

In Eq. (3), the first term Is the radiometer voltage signal due to the flux reflected from the target. The 
second term Is that due to the flux reflected from the surround; It Is detected by the radiometer because of 
the radiometer's Imperfect out-of-fleld-of-vlew rejection eharacterlstie. 


When the surface measured la shaded from the sun with a shield, the radiometer voltage reading Is 


''shad.t(9> - [PlEskyll + Ce(9)l + q]at(9) Ei 

+ [Edircose* + FlBskyU + Ss(9)l + Q]8a<«) K 2 . 

( 5 ) 

where Fi Is the fraction of scattered skylight falling on the 
target when the shield is In place and F 2 ^ corresponding quan- 
tity for the surround (Fig. 2). ?2 and Fi are larger than 0.9 whan 

the shaded solid angle la sm^l. We discuss the Fj coeffldent In a 
later section. In this case the factor Q la small and the effect of 
a change In Q with shading la negUglhle. 

Subtracting Bq. (6) from Eq. (3), we have 

Vdlr.t(9) - Vtot.t(9) - Vshad.t(9) 

- [SdircoeS* + (l-FijEskyU + «e<9)l]Rt<9)Ki 

.a-P2)E,kyU + «s(9)1S,(9)K2. ( 7 ) 

We use Eq. (7) In the following description of a technique for 
measuring the reflectance factor In the field. 



Fig. 2. Shading geometry. Part of 
the scattered Irradlance on both 
the target and the surround Is ob- 
structed. 


Reflectance factor measurement In the field 


To measure the reflectance factor, a radiometer is directed alternately at the target surface and at a cali- 
brated reference surface, for example a BaSO^ panel, and Che results are racioed. To eliminate the effect of 
scattered Irradlance on measurement accuracy, readings are recorded when Che reference and Che target surface 
— ^are shaded and then again when they are sunlit. Although the change In irradlance during the measurement pro- 
cedure introduces an error, the error can be compensated for by caking a second shaded reading after the 
reading and averaging Che two readings. To reduce errors, Che Vg^ad ^COC readings should be taken In 

rapid succession. Taking Che ratio of Che two secs of readings, we have 


'^tot.c(9) ~ Pshad.t(9) ^ '^dir.t(9) 

'^tot.r(9) * Vgtiad.rCS) '^dlr.r(9) 

^EjlrCOSSa + (l-Fl7Egky(l + 5e(9)ljEc(9)Ki + (1-P2)^kyll + 5s(9) |Re(9)K2 
[Edircossj + (l-Fl)Eg;tyIl + fr<9)l]Er<9)Kl + (l-F2)EgkyU + 5s(9)lH9(9)K2 ’ 


( 8 ) 


where 1^(9) has the same definition as Ie(9) and Is(9) but applies Co Che reference surface of reflectance 
factor Usually In field measurements, Ec(9) • Es(9) because the target is a part of Che surround, for 
example in measurements of part of an agricultural field. However, a reference target such as a BaS 04 panel 
has a reflectance factor considerably different from the surround. 

In both Che numerator and the denominator Che last term is small compared to Che first term. In estimating 
the ratio of the second to the first term in Che numerator we refer to Table 1, which lists Che spatial respon- 
sivlcy r(a) of a Barnes Model 12-1000 modular mulciband cadiomeCer8 as a function of field of view, and Table 
2, which Lists Che K 2 /K 1 values calculated by Eq. (4) for different fields of view. 


B 
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Table 1. Spatial respoosivlcy of a Table 2. K2/K1 values as a function 

Barnes radiometer, band 3. of field of view. 


Field of 
view (deg) 


r(o) 

Field of 
view (deg) 

K2/K1 


15 deg field 

1 deg field 

15 deg field 

1 deg field 

0 

1.000 

1.000 

1 


2.549 

4 

0.984 

0.719 

4 

— 

0.858 

10 

0.939 

0.167 

10 

1.2825 

0.053 

15 

0.208 

— 

12 

0.6327 

— 

20 

0.066 

0.001 

14 

0.1710 

— 

30 

0.002 

0.001 

18 

0.0213 

— 

40 

0.001 

— 

20 

0.0023 

0.004 




30 

— 

0.004 


If 


then 




?2 - 0.95, Egky - 0.2Edircoae*. Re(9) • M®). «t(«) - «s(9) “ 0, and K2/K1 - 0.06 (IS* field of view), 


(l~F2)B3fcy&g(9) K2 

[^dijcoaflj + (l“F0Eg|ty]Re(8)Ki 


0.06j',’ 


which is negligible. 


Thus, 


Vdlr.t(9) 

^ldirnn*9z * (l"Fl)E9ity{l 5t(9)l] 

&c(9) 

’dir.r<9) 

^Edircos9j + (l-POEakyd + Cr(9)J 

Rr(9) 


(9) 


If Che measured surface has approximately lamberclaa characceriacies , Chen (((S) and Ir(®) uearly egual, 
so Che reflectance factor of Che measured surface la 

at(9) - R(9) * Elf(9). (10) 

However, when Sa(9) is much larger chan Rc(®)> as for example between rows of plants in the near Q, when 
Eg^y Is very large, or when Che ouc-of-field response of the radiometer is high, the second term in the numera- 
tor in Eq. (8) cannot be neglected. The mld-IR bands of the Sames radiometer provide an example of how a high 
out-of-fleld response can give rise to large measurement errors. Previous work^ on measuring Che radiance of 
the exit port of an integrating sphere from different distances indicated that the ouc-of-fleld response was 
high. Our results, listed in the last column of Table I, agree with this observation. If we assume severe 
measurement conditions, for example, B3|(y/E^;COs92 • O.i, 3«(9) • &Rc(®)> 9nd a l-degres field of view, repre- 
senting a high out-of-fleid response (K2/Ki “ 2.549, see Table 2), a measurement error of 20J occurs even if 
Che shaded readings are taken into account. 

Analysis of the reflectance factor error 

Ue now discuss Che reflectance factor correction coefficient I(;(9) due to the nonlambertlan characteristics 
of the target surface and reference panel and due to the angular nonuniformity of the diffuse sky irradiance. 

Prom Che notation in Eq. (2), 


I + Ct(9) 


»/2 

Lsky(9) Et(9) aln9 cos9 d9 

0 


8t(9) 


»/2 

Lsky(9) ain9 cos 9 d9 

0 


( 11 ) 


Ic(9) can be estimated if Che hemispherical distribution of the skylight and the reflectance factor of the 
target surface as a function of solar zenith angle are known. 

To calculate some typical values for It(9), we refer to two atmospheric models, a model for a diffuse 
W panel, and several models for targets. The atmospheric models are those of Dave. 10 The first (Fig. 3) is the 

light Haze L model for X • 0.415 urn, a scattering optical thickness tj, - 0,1, and Irradiance angle 9-40 deg. 
The second (Fig. 4) la the Cloud Cl model for X • 0.825 urn, a scattering optical thickness T], - 0.1, and Irradi- 
ance angle 9-40 deg. For purposes of this calculation, we made some approximations to Che distribution 
curves and also assumed chat the curves do not change shape but only translate along the abscissa as the solar 


/• • - 
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zenith angle changes. Six dlffecenc normalized reflectance factor models were used to simulate several possi- 
ble situations (Fig. S). Model 1 is of a typical painted BaS04 diffuser and is used here for reference, ^dels 
2 through 6 represent five different possible surface reflection characteristics. 





Fig. 3. Downward scattered ra- 
diance distribution as a func- 
tion of' solar zenith angle for 
Dave model Haze L. 


Fig. 4. Downward scattered radi- 
ance distrlbucian as a function 
of solar zenith angle for Dave 
model Cloud Cl. 


Fig. 5. Six reflectance factor 
models used for calculation of 
the reflectance factor error. 


We define and AHtot 4S reflectance factor errors due to hemlapheric scattered Irradiance and non- 

lambertlan characteristics of targets with shaded readings subtracted (sHtUr) subtracted (4&tot)> 

Beferring to Eq. (9), Ideally Bfs) should be equal to Bt(^)/^r(B), which we will refer to as Big(s). How the 
error will be 


4BiUr<9) 


^Edlr*°**z * (l”?l)Esliy{ 1 * Ct(8)l 

I^EdirCoaSj + (l~Fi)Est(y{ 1 + CrC8)] 


1 »id(9) ; 


(U) 


referring to Eq. (3), 


Bld(9)* (13) 


From Eqs. (12) and (13), we see that the magnitude of the errors sEcot snd A^dlr 4s dependent on the ratio 
of scattered Co direct irradiance, the nonlambertlan characteristics of the measured surface, the reflectance 
factor of Che target, and the percentage of shaded irradiance, (l-F^). The should be much smaller than 
aB^qc because the effect of scattered irradiance has been significantly reduced. 


4^ot(9) 



|EdircosSz + E 

■sltyn 

+ Et(8)l + q](Ki + Kj) 


[e<iicco8 8i + Estyd + 5r(e)l 4- qJ 

Ki + 


Bia(8)K2 


- II 


The reflectance factor errors illustrated in Fig. 6 were determined by setting Egj^y equal to O.lE^ji^cosSj* 
The results show chat if Che ratio is used to determine Che reflectance factor, Che error ABco(/Big may be 





irradiance ancle (DECREES) IRRADIANCE ANCLE (DECREES) 


20 40 60 SO 

IRRADIANCE ANCLE (DECREES) 


20 40 60 so 

IRRADIANCE ancle (DECREES) 


Fig. 6. Expected percentage error in reflectance factor measurement due to scattered irradiance and nonlamber- 
tlan characteristics of targets. Curves (2) through (6) correspond Co surface reflectance models of Fig. S. 
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2Z CO 6Z under these condicloos. The largest measurement errors can be expected for short -wavelength (high 
atmospheric absorption), low-reflectance targets, larger solar zenith angles, and hazy or cloudy conditions. 
The average error for iRfoc/^d is about 2.3 Z for a - 0.05 and 1.5 Z for a - 0.75. However, If the signal ratio 
^dlr.t ^dlr.r is used in the determination of the reflectance factor, the error Is less chan 

0.4 Z according to the calculation, assuming that - 0.2E<iirco88x caking (1-Fi) values from Table 4 (for a 

shield area 76 cm x 76 cm at a distance of 3 m). 

In summary, the reflectance factor measurement can be analyzed mors accurately by using the ratios of the 
two terms instead of the two Vfgg(a) terms because (1) the atmospherically scattered Irradlance and the 

nonlambertlan characteristic of the measured surfaces have a negligible effect on the measurement, and (2) the 
effect of Che response of the Instrument outside its nominal field of view on the measurement Is significantly 
reduced by the subtraction of from Vggg. 

Field measurement of reflectance factor 


Ue measured reflectance factors on the roof of the Optical Sciences Center, Onlverslty of Arizona, on Novem- 
ber 1, 1983, a cloudless day with high vlslbUlty. We used several panels of known reflectance factors as tar- 
gets, in order to compare them with laboratory measurements. The panels had reflectance properties similar to 
chat of model 2 in Fig. S. The reflectance factor measured in the laboratory at 45 deg Irradlance angle and 0 
deg viewing angle, 8(4S/0)]^2l„ was used as the reference to compare the effect of reflectance factor on the 
measurement error. The relative reflectance factor errors for the direct (c^lr) And total (tcot) measurements 
are. In percent. 


«dlr 


^r - ^lab 
^lab 


Z 100 


and ceot - * 100, 

Alab 


where and Bcot Are Che reflectance factors measured in the field with the diffuse component subtracted 
(^r) And not subtracted (Bcot)< 

Figure 7 shows the measurement errors relative to the laboratory measurement results. Figure 7a shows chat 
ctot increases with the irrauUance angle. The curve is similar to Fig. 6 for reflectance model 2, but the 
measured error is larger than in Fig. 6. However, Is smaller chan ttot And does not increase with in- 
creasing irradlance angle. As the reflectance factor increases (Figs. 7h through d), the measurement errors 
C(ot And cdir decrease. The ctge iA il to 3Z on average, which is in agreement with the calculated results 
shown in Fig. 6. However, is less chan IZ on average. The results emphasize the need to subtract diffuse 
irradlance readings from total irradlance readings, particularly for low reflectance targets. 
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IRRADIANCE ANCLE (DECREES) 
(A) 


IRRADIANCE ancle (DECREES) 
(b) 


R|ab(A5/0) - 0 J38 
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10 


1 

5 
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0 



20 40 ‘^60 20 ‘ 40 60 
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Fig. 7. Percentage error of the reflectance factor in a measurement on November 1, 1983. 


Scattered and direct irradlance measurements 

The ratio s/d provides a useful way to quantify atmospheric conditions in remote sensing measurements. It 
may also provide a way to check the assumed atmospheric characteristics used in radiative transfer codes when 
Che latter are used to make spectral signature corrections. 

The measurement of s/d with two pyranomecers covers Che entire wavelength range of the detector's response. 
However, in many remote sensing applications, mulclspectral data are required. The spectral measurements of 
s/d should be made within a short time of each other for accurate inCercomparlson and analysis. This ratio can 
be conveniently measured with a mulclspectral radiometer and a uniform reference panel. 
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# Taking che ratio of Eqa. (6) and (7), considering that 5^(8) - 0, and ignoring Cs(9), we obtain 

'^3had.r(^) ^ (^l^aky + (EdirCosB* + F2Esky + <yBs(8)K2 

(Ejircosaj + (l-Pi)Eaky lRr(9)Kl + a-F2)Bskya9(9)K2 ' 


One way to eliminate the second term In the numerator (in the 
denominator che second term is very small and ran be neglected) is 
CO replace che reference panel with a material of extremely low 
reflectance factor, but of Che same size as the high-reflec- 
tance reference. Then can be calculated by Eq. (7). By 

connecting the two points and and extrapolat- 

ing CO R - 0, as shown in Fig. 8, we can determine che value of 
where 


"I “ (EdirCOS®t ?22sky + Q)Rs(®) ^2* 

Vhen a panel with a zero reflectance factor of che same size as 
che reference panel is viewed by a radiometer, che Vg|^^(S) valtie is 
exactly Che contribution due to che ouc-of-fleld-of-view response. 

Thus Eq. (14) can be rewritten as 


^shad.rC*) "'ll 
'fdlr.r<9) 


^l^ky * Q 


®dir«oa*z + (l”Pl)2sky 


(15) 


If n I is not subtracted in Eq. (IS), the percentage error A in- 
troduced in che s/d measurement la, from Bq. (14): 



REFIECTANCE FACTOR 

Fig. 8. A method to eliminate che 
surround effect on the s/d meas- 
urement. 


EjirCoaSi + F2Bsky <3 . ®s(8)K2 
t'lEsky + Q ttr<»)Ki 

I + F2(s/d) ^ Rs( 8)K2 
fl(3/d) ' 


(15) 


Table 3. Average reflectance 
factors for wheat and bare 
soil. 


where s/d “ (Egicy Q}/EdirC°8 92, which is che ratio of che diffuse 
irradiance, including skylight and Che surround reflection contribu- 
tion Q, CO che direct irradiance. 


The percentage error is dependent on che reflectance factor of 
che surround, che measurement geometry, and che ratio s/d. When 
che measurement geometry stays unchanged, che mesauremenc error 
increases with s/d and with che reflectance factor of che surround. 


Band 

No. 

Wavelength 

(um) 

Wheat 

SoU 

1 

450-520 

0.019 

0.058 

2 

520-600 

0.028 

0.083 

3 

630-690 

0.017 

0.120 

4 

760-900 

0.495 

0.181 


To estimate che effect of che surround on che s/d measurement, 
we carried out an experiment on March 20, 1984, at che Maricopa 
Agricultural Center, Arizona, in which two Barnes mulciband radiome- 
ters were used Co make identical measurements in locations with 
different surrounds. One was placed at che center of a 100-acre 
field of bare soil and che ocher at che center of a 100-acre field 
of wheat. Table 3 lists che average reflectance factor values for 
che wheat and Che bare soil. 

According Co Eq. (16) and che data in Table 3, we can predict 
Chat s/d should be increasingly greater when measured in che bare 
soil field chan when measured in che wheat field for TM bands 1, 2, 
and 3 because of che higher reflectance factor of bare soil in that 
wavelength range. However, for TM band 4, the measured s/d should 
be greater in che wheat field than in che bare soil field owing to 
che higher reflectance factor for wheat in che IR spectral region. 
This prediction . was confirmed by field measurements, as shown by 
Fig. 9. For example, at larger zenith angles, Che value of s/d is 
about SZ greater for bare soil chan for wheat in band 3 but 20Z 
less in band 4, 

It is not difficult to deduce from Eqs. (14) and (16) chat 



solar zenith ancle (DECREES) 


I 


Awheat 
1 * ^soil 


'^shad.r 

'^dir.r 

'^shad.r 

'^dlr.r 


measured in wheat field 


measured in soil field 


Fig. 9. Surround effect on che s/d 
measurement in an experiment on 
March 20, 1984. 
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Figure 10 shows both the curves for (l+4uheat)/O'*’480ll) calcu- 
lated for different K2/K1 values, as given by Eq. (16) In band 4, 
and the ratio of Vjjjad.r/'^dir.r measured In the wheat field to 
^shad.r/^dlr.r measured In Che bare soil field. By comparing the 
calculated and measured data for different solar zenich angles, we 
find chat the measured ratio lies between K2/K1 values of O.OS and 
about 0.07, which is higher chan the Kj/Ci value of 0.04 that was 
determined from the measurement geometry. The discrepancy la prob- 
ably due to Che difficulty of accurately determining E2/K1. 

Equation (16) suggests that Che s/d measurement should be car- 
ded out where the surround has a low reflectance factor and/or the 
reference panel Is large enough to make the ouc-of-fleld concdbu- 
don CO Che radiometer signal negligible. 


Setting M equal to (Vaha/i r(9) 
have 


nl)/Vdir.r(®) 2q. (15). wo 


and 


s/d 


^l(^ky 0) 

^rCoaSj + (l-Fl)(Eaky Q) 

H 

’ ^1 - (l-Fi)k • 


Pl(8/d) 

1 + (I-Fi)(s/d) 


(17) 


(18) 


The magnitude of F i depends on the aCmoaphedc condidona and 
Che measurement geometry (Fig. 11). If we know the seaccerad radi- 
ance dlscdbucion f(a,4) and the geometde measurement parameter a’, 
we can determine Fi by the equadon 


l-F, 


a'/slnss 
-a'/slnsj 1 


9*+a' 


f(s,4) sine cose de de 


8*-e' 


(19) 


2» 


0 . 


f/2 


f(e,e) sine cose de de 



Fig. 10. Bado of (l+Awhoat)/ 

(I-eSsoll) es a funcdon 
of solar zenith angle. 



By measudng acmosphedc parameters (acmosphedc exClncdon co- 
efficient, air pressure, temperature, etc.), we can calculate the 
scattered radiance dlscdbudon function f(a, 4 ) from a radladve 
transfer model. 

The assumpdon chat Che scattered radiance Is uniformly dlstdb- 
uced leads to underesdmadng the (l-Fi) value because the diffuse 
radiance la much higher In the region of the solar aureole chan In 
Che rest of the sky, especially on hazy or cloudy days (see Fig. 3). 
The small area of the shield Chen shades the part of the scattered 
sky with Che strongest radiance peak. 

To calculate (l-F^), the atmospheric models were used as before. 
The (1-Fi) values are listed In Table i. These values are for a. 
shield size of 76 cm x 76 cm and for a distance between Che mid- 
point of Che reference panel and Che shield of 3 m. Table 4 also 
lists a (1-F^) value of 0.010 for an assumed uniform diffuse radi- 
ance dlscdbudon. 

The corrected s/d value may now be determined from Eq. (18) 
using Che known F^ and the measurement data. 


Fig. 11. Calculation geometry 
for (l-Fj). 


Table 4. Scattered skylight shaded 
coefficient (1-Fi) value. 


Irradlance 


angle (deg) 

Haze L 

Cloud Cl 

20 

0.158 

0.172 

30 

0.134 

0.137 

40 

0.116 

0.108 

40 

0.116 

0.108 

50 

0.100 

0.085 

60 

0.085 

0.065 

70 

0.067 

0.045 

Uniform 

dlscdbudon 

assumpdon 


0.010 


Measurement of acmosphedc extinction coefficient 

Solar radiometers are used to measure the acmosphedc eztlncdon coefficient. They have a narrow field of 
view (usually 1 co 2 deg), which Is pointed at the sun. Readings of direct solar Irradlance at Che ground are 
made at known dmes during Che day. The solar zenith angle corresponding co each dme of day is calculated and 
Che acmosphedc extlncdon coefficient Is Chen determined by Che Langley plot mechod. 


The measuremenC of solar Irradlance at Che ground using a radiometer and a reflectance panel need not In- 
volve changing che viewing angle. The measuremenC geometry was Illustrated In Fig. 1. By caking a sedes of 
global and diffuse Irradlance measuremencs in rapid succession and at known dmes from a solar zenith angle of 
75 deg CO Ics smallesc values chac day, we can again determine hy che Langley plot mechod. 




From Eq. (7) we have 


^ 99-1 


'^dlr.r<®*) * ^EiUrcoafl* + (l-Fi)Egity( 1 + tr(®r)lJSr(9*)Kl . 

Remembering chat Che cerm (1-Fi) is small (usually less than 0.1), and assuming ehac Sgicy ^ approztmacely 
proportional to E^^cosSj with a constant* coefficient Kg and chat Er(®z} * 0, we obtain 


BdlrSl 


'^dlr.r<9z) ^ I 

Sr<®z)cosez 1 + (l-Fi)B:g • 


( 20 ) 


Substituting Ej^,. a E4iro*^(~'^ezc*ac«z) into Bq. (20), where Ediro la the exoacmospherlc solar Irradlance, we 
have 


Ediro^lEl + (l-Fl)Kel * ezp(-Texcsacaj) 


and finally 


■^ext 
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'^dlr.r(9z) 

Rr<®z)coa6* 
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^dtr.r(®z) 

Rr(®z)cos6z 


( 21 ) 


( 22 ) 


If In Vdj.r.r(az)/Rr(®z)coa ®z la sec as the ozdlnats and sacSz aa the abscissa, the slope of the regression 
straight line Is exactly the Tg^c value. The nuaerlcal calculation technique used for the determination of 
Text 'vas chat proposed by Herman ec al.^^ The reflectance factor value of the reference panel should be meas- 
ured In the laboratory before and/or after the measurement. 


In comparing the two methods of tgxc measuresMnt, two points should be made. The first Is the effect of 
scattered Irradlance on the measurements. Since the solar radiometer has a small field of view, only Che 
atmosphere In the field of view will directly affect the measurements. In the Langley plot, the points that 
diverge from the straight line owing to occastonal atmospheric changes within Che field of view of Che solar 
radiometer can be discarded, thus providing a higher correlacion coeffldent for Tgxf However, when the radi- 
ometer measures the reflected radiance, the fluctuations in signal due to irradlance fluctuations in the hemls- 
^ pherlcal diffuse component oiay weaken the correlation coefficient of the regression line. Thus the measurement 

method using a reference panel may not be as accurate as chat using a solar radiometer. 


The second point is chat atmospheric gaseous absorption can be 
avoided by the appropriate choice of the center wavelengths of the 
narrow-spectral -band Inteference filters used In the solar radiome- 
ter. In this case, Tg^t Is only the scattered extinction coeffi- 
cient. However, not all field spectroradlometers are designed this 
way. For remote sensing use, the spectral bands of the field radi- 
ometer are usually chosen to match chose of space mulclspeccral 
sensors, which In some cases Include atmospheric absorption bands, 
for example, Che 0.8 to 1.1 um band of the Landsac mulclspeccral 
scanner systems. So Tgxc, determined by the reference panel 
method, will In some cases Include both Che scattered and absorp- 
tion extinction coefficients due to atmospheric gases, especially in 
the IR where water vapor, carbon dioxide, oxygen, ozone, etc., are 
strongly absorbing. 


To check the accuracy of the reference panel method for decer- 
^ext, ''B up a solar radiometer and a Sames modular mul- 
clband radiometer In the same location in Tucson, Arizona, on Sep- 
tember 21, 1983. A 3aS04 panel, Che reflectance factor of which 
had been measured In Che laboratory, was used as a reference. 
Figure 12 shows the results. Note that the two secs of tgxt data 
are consistent through the visible spectrum but that the reflec- 
tance panel method produced a higher rg^^^ value In Che near IR, 
possibly because of water absorption. 



WAVELENCIH (nm) 


Fig. 12. Comparison of Tgxt meas- 
urement with a solar radiometer and 
a Barnes radiometer. 


Concluding Ramarks 

Equations have been derived for use In determining Che reflectance factor, the ratio of scattered to diffuse 
ground Irradlance, and the atmospheric extinction coefficient. The equations account for the nonlamberclan 
characteristics of the target, nonunlform sky radiance distributions, reflectance of the surrounding neighbor- 
hood, and ouc-of-fleld radiance contributions. 


*The error Introduced by this assumption 
changed less chan 10 Z from sunrise to midday. 


Is not significant. In our summer cloudless measurement, 
which leads to an error of less chan IZ. 


Ks 



Nonunliona sky Irradlance and nonlarabertlan reflecclon characcerlsclcs give rise to a reflectance factor 
measurement error of 2Z to 6Z depending on the angular variation of the reflectance factor and atmospheric 
conditions. This result Is In general agreement with ocher authors. The method proposed in this paper for 
the measurement of reflectance factors makes this error due to nonuniform sky irradiance and out-of-field 
effects negligible. 

An error on the order of 20 Z can be encountered even If the shaded readings are taken Into account in the 
determination of reflectance factor, because of the high ouc-of-fleld response of a commonly available field 
radiometer when operated in its mld-IK bands. 

A method and a correction procedure were described for the accurate determination of the ratio of scattered 
to direct (s/d) ground Irradiance. The results showed that a multlspectral radiometer and a reference panel 
can give atmospheric extinction coeffldents in agreement with solar radiometer data In the visible spectrum. 
In the IR the results differ, possibly due to water vapor or other absorption in the spectral bands of the ra- 
diometer. This procedure should allow Tq^c l-o visible to be determined by many researchers who have 
access to a multlspectral radiometer and a reference reflectance panel but not to a solar radiometer. 
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Factor in Remote 


Radiometric Considerations in Remote 
Sensing 


PHILIP N. SLATER 

Invited Paper 


The need for accurate radiometric data for the verification and 
use of scene radiation models is emphasized. The radiometric 
problems associated with reflectance and atmospheric correction 
field measurements and sensor calibration are reviewed, estimates 
are made of the attainable accuracy in each case under favorable 
conditions. The loss in radiometric accuracy by resampling proce- 
dures in digital imagery processing is discussed. 

I. Introduction 

Coincidental with the development of comprehensive 
models of the interaction of electromagnetic radiation with 
the earth's surface and atmosphere, there is a growing 
interest in the application of radiation models to global 
earth-atmosphere-ocean studies such as COSPAR's Inter- 
national Satellite Land Surface Climatology Project and 
NASA's Global Habitability program. 

Most scene radiation models desctibed in recent reviews 
by Bunnik [1] and Smith [2] predict an upward spectral 
radiance, just above the feature of interest, on the basis of 
certain input values and then use an atmospheric model to 
predict the spectral radiance at the aircraft or satellite sensor. 
For example, in the case of a crop canopy, the input values 
usually include the leaf spearal reflectance and canopy 
geometry (leaf orientation and leaf area index, LAI), the soil 
spectral reflectance, and the geometry of irradiation and 
viewing. For many applications only the inverted form of 
this model is relevant; Given the radiance, irradiance, and 
viewing geometry, what is the constitution (LAI, vigor, etc.) 
of the canopy? Reference [3] is the latest of several publica- 
tions by Goel and coworkers in which the result of such a 
procedure are described, but little other work on scene 
radiation inversion models has been reported. 

The use of satellite-acquired image data with scene radia- 
tion models is inextricably dependent on atmospheric cor- 
rection and sensor absolute radiometric calibration. Unlike 

Manuscript received lanuary 15, 1985: revised February 21, >965. 
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computer-aided scene classification, which depends on a 
statistical analysis of the digital counts in a scene image, 
usually without correction for the intervening atmosphere, 
modeling is concerned with determining the radiance of 
the scene. For this purpose we need to know, first, the 
output digital counts from the sensor when it is imaging the 
feature of interest; second, the absolute radiometric calibra- 
tion of the sensor in order to convert the digital counts to 
radiance at the entrance pupil of the sensor; and third, the 
radiance modification introduced by the intervening atmo^ 
sphere (the atmospheric correction) in order to relate the 
entrance pupil radiance to the radiance of the ground 
feature [4], For verification of models and for purposes of 
atmospheric correction to an uncertainty of a few percent, 
we need to be able to measure the bidirectional reflectance 
of uniform ground areas to an uncertainty of about 1 
percent. We also need to be aware of the radiometric errors 
introduced if the digital image has to be resampled for 
spectral band and/or map registration purposes. 

In this paper, we review the radiometric problems associ- 
ated with field reflectance and atmospheric correction mea- 
surements and with sensor absolute calibration, and esti- 
mate the accuracies currently attainable. We then discuss 
the loss in radiometric accuracy introduced by resampling 
procedures in digital image processing. 

II. Scene Reflectance measurement 

In computer-aided scene classification we use statistical 
differences among pixel digital counts to distinguish earth 
surface features. In this context we are usually interested 
only in relative reflectance values (an exception is multi- 
temporal classification); indeed, the spectral radiant flux 
reflected by, and characteristic of, the features of interest 
may be only a small fraction of the total flux at the sensor, 
depending on the ground reflectance, solar zenith angle, 
viewing geometry, wavelength interval, and atmospheric 
conditions. In using models, on the other hand, we are 
often interested in knowing, or determining, the bidirec- 
tional radiance or reflectance distribution function of a 
surface feature, and we either correct for atmospheric ef- 
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fects or consider only the reflected radiance distribution 
immediately above the surface. 

Before we discuss some of the problems in measuring 
reflectance, it is appropriate to review the definitions and 
interrelationships among the various ways of describing the 
reflecting properties of a surface that are relevant to remote 
sensing. 


A. Reflectance Measurement Geometries 

Of the nine different geometries that can be used to 
describe the reflection of radiant flux from a surface [5], 
fortunately only four are used to any great extent in remote 
sensing. These are: 

1) the directional-hemispherical reflectance factor 

2) the bidirectional reflectance factor 

3) the directional-hemispherical reflectance p(9,,<)>,; 2ir) 

4) the bihemispherical' reflectance p(2ir;2ir) 

where 0 and <> are the angle to the surface normal and the 
azimuthal angle, respectively, and the subscripts i and r 
refer to incident and reflected angles. It is often Important 
to relate these first four geometries to a fifth description: 

5) the bidirectional reflectance distribution function 
(BRDF) f( 

Reflectance is defined as the ratio of the radiant flux 
reflected from a surface to the radiant Mux incident on the 
surface. In the most general terms, it is the ratio of the 
radiant flux reflected from a surface into a hemisphere, to 
the flux incident hemispherically onto the surface; it is 
entry 4) in the above list; Reflectance factor is the ratio of 
the radiant flux reflected from a surface to the radiant flux 
reflected from a lambertian surface of unit absolute reflec- 
tance when both measurements are made under identical 
irradiance and viewing conditions. 

The measurement of the directional-hemispherical re- 
flectance factor may be considered the fundamental reflec- 
tance measurement because it provides the calibration of 
reference surfaces used in remote sensing field and labora- 
tory work. Using refinements to the Van den Akker auxiliary 
sphere method, Venable et at. [6] at the National Bureau of 
Standards (NBS) were able to measure the directional- 
hemispherical reflectance factor of high-reflectance surfaces 
with an estimated uncertainty of ±0.15 percent. Unfor- 
tunately, the refined Van den Akker method involves com- 
plicated measurement and data-reduction procedures. Mea- 
surements of R(d;,<^-,2rr), or the equivalent R{2ir, 9, ,((>,), are 
routinely made in radiometric and colorimetric laboratories 
by use of a dual-beam spectroradiometer or spectropho- 
tometer [7], (8, pp. 174-180). These instruments provide a 
beam of narrow spectral bandwidth that is incident alter- 
nately on the surface of interest and on a reference surface 
whose calibration 's traceable to the NBS. The alternative 
signals from a detector, placed in another port in the 
sphere, are then ratioed. This ratio, when multiplied by the 
absolute reflectance of the reference surface, measured 
using the same irradiance geometry, gives <^,; 2ir). 

In remote sensing field work, directional-hemispherical 
and bihemispherical reflectance measurements have been 
made for reflectance studies and for studies of atmospheric 
and surface energy budgets (9]-(12). However, the 2w solid 


angle collection geometry approximately simulates the 
averaged response over the total field of view only for 
wide-field acquisition systems, such as the Advanced Very 
High Resolution Radiometer. The Landsat Multispectral 
Scanner System (MSS) and the Thematic Mapper (TM) are 
better simulated by a bidirectional or hemispherical-direc- 
tional measurement. Some correctly refer to "conical" in 
place of "directional" when referring to measurements 
made with a hand-held radiometer (see, for example. Kimes 
et al. [13D- In most cases, the solid angle subtended by the 
apePure of the radiometer at the surface of interest is very 
small, on the order of 10" "* sr. In these cases, and certainly 
for aircraft and spacecraft sensors, "directional" is a suffi- 
ciently accurate geometrical description. 

The instrument that is usually used to measure the hemf- 
spherical-directional reflectance factor is a nadir-viewing, 
narrow-field (1“ to 15®) multiband radiometer. Held typi- 
cally 1 to 5 m above the ground, it is moved across the 
feature of interest, and many output voltage readings are 
taken. These are averaged and ratioed to the output voltage 
reading obtained when the instrument is viewing a near- 
lambertian, horizontal reference panel. When data for the 
reflectance factor of the reference panel have been cor- 
rected for the panel's non-lambertian characteristics, the 
reflectance factor for the unknown surface is determined 
from 

/?,(2ir;0,0)--^R,(2rr;0.0) (l) 


or, if we make the approximation of neglecting irradiance 
due to the sky. 


where is the solar zenith angle and the radiometer's 
optical axis is parallel to the surface normal; and V, are 
the voltages from the radiometer (with the dark voltages 
subtracted) when the instrument is viewing the target and 
reference, respectively; and R,(tf,,i>,;0,0) is the reflectance 
factor of the reference panel with respect to a lambertian 
surface of unit reflectance. 

fof some surfaces can be measured in the 
laboratory as well as in the field. In this sense, it is more 
general than R ^{2^: because the distribution of dif- 

fuse irradiance for an outdoor measurement is difficult to 
simulate accurately in the laboratory. We will return to this 
point in the discussion of (9). Fortunately, for a near-lam- 
bertian reference panel, Rp{2ir; 9,.<^,) and Rp( 9;. 9,, <f>,) are 

•very close over the range 40® < 9, < 60® especially at 50®. 

The hemispherical-directional reflectance can be de- 
termined in a similar manner by using an upward-looking 
pyranometer, which views 2ir steradians, in place of the 
reference panel. In this case, the multiband radiometer and 
the pyranometer have to be accurately calibrated in an 
absolute sense. The hemispherical-directional reflectance 
of the target is given by 


p{2if,9,.^) 


Lr{9,.<^,) 

t,(2ir) 


( 3 ) 


where t, and L, are the radiance value readouts for the 
multiband radiometer and the pyranometer, respectively. 

The bidirectional refleaance factor <(),; 9,, is 
related to the directional-hemispherical reflectance 
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p(0;,<t>i; 2v) of a surface by 

p(^/-'^»2ir) • 2 9,.i^,)i\nd,coidrdd,. (4) 

■'0 

Hsai and Weidner [14] have shown that the reflectance of a 
reference panel can be determined to an uncertainty of 0.5 
percent by a numerical integration of (4) if R(9,,i^: is 

measured in 10® steps form 10® to 60®. The further a surface 
departs from being lambertian, the greater the number of 
measurements that must be made. 

A complete description of the reflectance properties of a 
surface is provided by the spectral bidirectional reflectance 
distribution funaion (BRDF) f{9,,<^,;9„<^;\). which relates 
the directional spectral radiance reflected from a surface to 
the directional spectral irradiance on the surface. BRDF is 
defined by 




dl;(9,.<pi; 9,,<^) 

dE{9„i,,) 


(5) 


and has the units of reciprocal steradians; dL, is the radi- 
ance reflected from the surface in the direction and 
d£ is the irradiance incident on the surface in the direction 
9;,<p;. BRDF is related to the bidirectional refleaance factor 
by 

“ rrf(9;,^;9,,<^). ( 6 ) 

Also it can be shown that 


p(9i,<^,;2v) - R(tf,;<h;2») 

(7) 

p(2v;2v) ■ R(2v;2ir). 

(8) 


The rest of this section deals with irradiance and surface 
effects on field reflectance measurements, shortcomings of 
instrumentation and measurements, the estimated accuracy 
of field reflectance measurements, and the need for the 
comprehensive description of measurement conditions. 


Table 1 Percent Change of the Reflected Radiance Due 
Either to a Change of the Distribution of the Irradiation by 1® 
of Solar Zenith Angle or by a 10 Percent Change of the 
Optical Depth of the Atmosphere. Averaged Over All 
Directions of Reflection and Over All Distributions of the 
Irradiation [16) 



Average Change of the Reflected Radiance 

Surface Type 

per degree 
change 
of the solar 
zenith angle 

per 10-percent 
change of the 
optical depth 

Savannah 

+ 1.0 

+ 1.6 

Bog 

+ 0.9 

±0.7 

Pasture land 

+ 1.7 

±1.0 

Coniferous forest 

±2.3 

±1.5 

Average over 
the four surfaces 

+ 1-3 

±1-2 


by averaging target measurements taken immediately be- 
fore and after the panel measurement. 

Irradiance fluctuations, reported by Slater [8, p. 315], of 1 
to 2 percent at frequencies on the order of 1 kHz, possibly 
due to fast-moving subvisuai cirrus clouds, would affect 
measurement accuracy only if the radiometer had a fast 
time response. Most do not. This phemonenon would nev- 
ertheless have the effect of increasing the noise in scanner 
imagery. Fortunately its occurence, at least in the author's 
experience, seems to be rare. 

Non~ Lambertian Surface Effects: The combined effect, 
on the accuracy of reflectance factor measurements, of 
anisotropic sky irradiance and a non-lambertian target and 
reference panel has been reported by Kirchner ef al. [22] 
and Robinson and Biehl (23). Their results generally agree 
with the analytical results of Che ef al. [24], who modeled 
several different non-lambertian surfaces and used Dave's 
[25] haze and cloud models to yield the results listed in 
Table 2. 


B. Irradiance and Surface Effects 

The measurement of the reflected flux from both the 
target and the reference panel is influenced by anisotropy 
in the hemispherical spatial distribution of irradiance due to 
atmospheric turbidity and clouds and changes in the solar 
zenith angle. Kriebel [15]-[17] examined the influence of 
changes in atmospheric turbidity and solar zenith angle on 
the radiant flux reflected by vegetated surfaces. He mea- 
sured the angular distribution of the spectral reflectance of 
four vegetated surfaces having significantly different canopy 
structures and used a radiative transfer program to simulate 
different atmospheric conditions. His results are listed in 
Table 1 for a wavelength of 0.52 pm. The change in reflec- 
tance with change in solar zenith angle is almost wave- 
length independent whereas the dependence of reflectance 
on optical depth decreases with increase in wavelength. 

Temporal Irradiance Changes: Duggin [11], [18] and Mil- 
ton [19], [20] have disputed the importance of simulta- 
neity of target and reflectance panel measurements. Un- 
doubtedly, simultaneity is mandatory when clouds can 
modify the irradiance level between sequential target and 
panel readings [21]. At large solar zenith angles the ground 
irradiance changes rapidly with time. Simultaneity in this 
case is desirable but not necessary because the irradiance 
change is smoothly continuous and can be accounted for 


Table 2 Maximum Percentage Errors in Hemispherical-Di- 
rectional Rellectance Factor Measurements for Five Non- 
Lambertian Surfaces Under Haze and Cloud Conditions [24| 


Solar 

Zenith 

Dave Haze t 

Dave Cloud O 

Model. \ 

- 415 nm 

Model. K 

— 825 nm 

Angle 

R - 0.05 

R-0.75 

R - 0.05 

R - 0.75 

20 

2 to -1 

3 to -1 

3 toO 

4 toO 

70 

5 to -3 

4 to -2 

6 to - 2 

5 to - 1 


Solar Zenith Angle Effects: The effect of solar zenith 
angle on reflectance or reflectance factor measurements has 
been investigated by several workers. Coulson and Reynolds 
[12] measured bihemispherical reflectance as a function of 
six vegetation canopies for six discrete wavelengths in the 
visible and near-infrared. They found that most surfaces 
showed a maximum bihemispherical reflectance between 
70® and 80® that gradually decreased as the solar zenith 
angle decreased. Egbert and Ulaby [26] measured the spec- 
tral hemispherical-directional reflectance factor for many 
solar zenith angles and nadir and azimuth viewing angles as 
part of a procedure to determine the best conditions for 
multiband photography. Kimes et al. [13] found that the 
hemispherical-directional reflectance for lodgepole pine 
and two grass canopies decreased as the solar zenith angle 
increased, under clear conditions. In a simulation of the 
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conditions, they found that the bihemispherical reflectance 
behaved in the opposite way, in agreement with Coulson 
and Reynolds. The difference was explained on the basis of 
the anisotropic BRDF of the simulated canopy. Ouggin [27] 
and Pinter et a!. [28] found that the effect of solar zenith 
angle of the bidireaional reflectance factor of several wheat 
varieties at the same stage of growth depended on the 
variety. Jackson ef al. [29] showed, for a single wheat 
variety, that the effect of solar zenith angle on hemispheri- 
cal-directional reflectance factor depended on row spacing 
and orientation, plant height, and canopy cover. Kirchner et 
al. [30] have modeled the conditions of anisotropic BRDF 
and irradiance distribution to predict direaional radiance 
for a wide range of conditions. Schnetzler [31], using data 
from Kirchner et al.. was able to predict radiance at the 
entrance pupil of a pointable orbital imager under a wide 
variety of conditions. 

Excellent sources for further information on ground mea- 
surement techniques and results are the proceedings of the 
international colloquia entitled Spectral Signatures of Ob- 
jects in Remote Sensing, edited by Cuyot and Verbrugghe 
[32], [33]. Space does not permit discussion of all the material 
in these proceedings; the reader is referred in particular to 
articles by Becker, Cuyot, Jackson, Ott et al., Vanderbilt et 
al., and Verhoef et al. 

C. Instrumentation and Measurement Shortcomings 

Field reflectance measurements are used most extensively 
in studies of changes in conditions of surface features. In 
this case, only relative reflectance values are usually re- 
quired; instrumentation deficiencies, anisotropic irradiance 
conditions, etc., generally introduce negligible errors in the 
results. Field reflectance measurements are sometimes com- 
pared with laboratory reflectance data or used, as discussed 
later, as part of the input to a method for determining the 
in-flight absolute radiometric calibration of sensor systems. 
In these two cases, and in the verification and use of scene 
radiation models, accurate absolute reflectance data are 
usually required. 

Robinson and coworkers [7], [23], [34], Collins [35], and 
Deering and Leone [36] in this country, Ahern ef al. [37] in 
Canada, and Bunnik ef al. [38] in The Netherlands have 
developed new instruments and advanced the techniques 
of accurate, convenient, and reliable field reflectance mea- 
surements. Their publications should be referred to for 
descriptions of instrumentation design and of operation 
and measurement methodology. The following discussion 


emphasizes some of the problems related to instrumenta- 
tion and measurement shortcomings. 

Avoidance of Effects Due to Sky Anisotropy: To com- 
pare field reflectance factor data with laboratory refleaance 
factor data or to minimize anisotropic sky irradiance effects, 
a simple approach is to take radiometer voltage readings of 
the target and reference surfaces when both are globally 
irradiated and then when both are shaded from direct solar 
irradiance. If we refer to the target and reference global 
voltage values as ^rc and and the target and reference 

shaded values as '^rs and then the bidirectional reflec- 

tance factor can be obtained from 

R(g;,»,;0,0) - X (9) 

^RC *RS 

where R„(tf;,<^;0,0) is the reflectance factor of the refer- 
ence panel with respect to a lambertian surface of unit 
reflectance. Equation (9) is then a more accurate form of (2). 

Che et al. [24] used several panels whose reflectance 
factors were measured in the laboratory and then in the 
field. The field measurements were made by using global 
irradiance (sun and hemispherical sky) in conjunaion with 
(1), and then by shading the target and reference panel and 
using (9). Fig. 1 shows the measurement errors relative to 
the laboratory measurements. The results show that, for a 
reflectance factor of about 0.09, the error can be as high as 
10 percent but can be reduced to 4 percent by using the 
shading method. The error shows an increase with solar 
zenith angle and a decrease with increasing reflectance. 
The residual error, when (9) is used, is dependent on the 
geometry of the measurement and the prevailing atmo- 
spheric conditions. For example, when the atmospheric 
aerosol loading increases, the solar aureole increases In size 
and radiance at the expense of the direct solar irradiance on 
the ground. The geometry of field measurements under 
these conditions will differ from a bidirectional laboratory 
measurement, and the difference between the results is 
likely to be greater than under low aerosol loading condi- 
tions. 

These experimental results of Che et al. are not in agree- 
ment with the predictions of Bauer et al. [39], who esti- 
mated the measurement uncertainties under conditions of 
moderate and of very noticeable haze. Admittedly, the haze 
conditions differed and the assumptions made by Bauer et 
al. regarding typical BRDF values and the uncertainties 
associated with field measurements may not have applied. 
Nevertheless, it is not obvious that these differences can 




R|^4S/01 *0.174 


10 

■ 

mm 

1 

BH 

Essm 

■ 



H 



23 


60 


20 


60 


miMOIANCE ancle (OECREES) 

(a) 



iRItAOMNCE ancle (OECREES) 

(b) 


23 <0 

irradiance ancle (DECREES) 


(C) 


20 <0 60 
IRRAOWnCE ancle (DECREES) 

(d) 


Fig. 1. Percentage error in reflectance factor measurements (from Che er a/. (24|). Note 
that the bidirectional reflectance factor measured in the laboratory at 45° irradiance angle 
and 0° viewing angle was used as the reference. 
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explain the disagreement. A study should be undertaken to 
reconcile the differences. 

Proximity of Nearby Objects: Ouggin [18] has pointed 
out the need to prevent secondary reflectances from falling 
onto the measurement areas and to ensure that the up- 
ward-looking pyranometer, for measuring global flux, is 
perfectly horizontal to minimize the integration of flux from 
reflecting surfaces below the horizon. This last point is. of 
course, also important in the use of a reference panel. 

Kimes et at. [40] have investigated the effect of a nearby 
object on reflectance factor data as a function of the solid 
angle subtended by the object at the surface to be mea- 
sured and as a function of the reflectance of the object and 
solar zenith angle. Both the object and surface were as- 
sumed lambertian with the object on the opposite side of 
the surface from the sun. Dave's (251 clear-sky atmospheric 
model was used for the calculations. The results, giving the 
percentage error in the spectral radiance of the surface due 
to the proximity of the object, are shown in Fig. 2. The 


Difference Between Actual and Nominal Fields of View: 
An insidious problem, with at least one multiband radiome- 
ter on the market, is that its field of view, when set to be 1" 
for TM band 5, extends well beyond the nominal field. Che 
et al. (24) measured the field of view and found, for exam- 
ple, that normalized to unity for 0® the response was as 
high as 0.7 at 2“ off-axis. They estimated that, under hazy 
conditions (a diffuse to direct irradiance ratio of 0.4) and 
with a surround four times the reflectance of the target, the 
error in the reflectance factor due to the high out-of-field 
response would be as high as 20 percent. 

Surround Effect on Irradiance: Another insidious prob- 
lem, related to the above, is the effect of the refleaance of 
the area surrounding the surface being measured. Che et al. 
(24) observed this effect when determining the ratio of the 
diffuse to direct irradiance on panels placed at the centers 
of two 100-acre fields, one of bare soil and the other of 
wheat. Fig. 3 shows how some of the upward radiance from 
the surround can be scattered by the atmosphere to be 




Fig. 3. Measurement geometry with a field spectroradi' 
ometer 124). 
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Fig. 2. Simulated percent error of radiance measurement 
due to the presence of a nearby object (opposite the sun's 
direction) as a function of solar zenith angle, wavelengths of 
0.675 and 0.796 jim. black or white objects, and the solid 
angle intercepted by the object (40). 


errors increase with increase in solar zenith angle and solid 
angle. The results show that a person in white clothing 
kneeling O.S m from the surface being measured causes an 
error of 18 percent for a wavelength of about 0.8 (im and a 
solar zenith angle of 75®. But if the person wears black 
clothing, the error is less than 2 percent for all solar zenith 
angles. The error is strongly dependent on the solid angle. 
So. when making field reflectance measurements, the rec- 
ommendation is to wear dark clothing and stay away from 
the surfaces while they are being measured. 


Table 3 Average Hemispherical-Directional 
Refleaance Faaors (or Wheat and 
Bare Soil (24) 


Band 

No. 

Wavelength 

(nm) 

Wheat 

Soil 

1 

450-520 

0.019 

0.053 

2 

520-600 

0.028 

0.083 

3 

630-690 

0.017 

0.120 

4 

760-900 

0.495 

0.181 


incident on the panel being measured. From Table 3, which 
lists the hemispherical-direaional reflectance factors for 
the two fields for the first four bands of the TM, we can 
predict that the ratio of diffuse-to-direct irradiance will 
increase more for the bare soil than for the wheat field with 
increasing wavelength for the first three TM bands. How- 
ever, for band 4 the ratio should be greater in the wheat 
field than in the bare soil field owing to the higher reflec- 
tance factor for wheat in the IR spectral region. This predic- 
tion is confirmed by the results of the field measurements 
as shown In Fig. 4. Note that these results also include the 
out-of-field effect mentioned earlier. 

Calibration of Reflectance Factor of Reference Panel: 
Many investigators calibrate their field reference panels 
with reference to BaSOj or Halon powde; pressed accord- 
ing to specifications suggested by the NBS. Some technique 
descriptions in the literature imply that the authors have 
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fig. 4. Effect of the surround in the measurement of the 
ratio of atmospherically scattered flux to direct solar flux (24|. 

assumed these pressed powder surfaces to be lambertian 
and that the directional-hemispherical reflectance has been 
used in calculating the reflectance factor of the reference 
surface. A substantial error is incurred if these assumptions 
are made, as can be seen by reference to Table 4. 

Table 4 Values for the Bidirectional Reflectance Factor 
<>,;0,0) for Halon. A is the Oifference Between the 
Quantity Immediately Above and Value of the 
Directional-Hemispherical Reflectance p(0.0;2ir) for Halon 
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Temperature Effects on Radiometer Response: The 
change in detector response with temperature, particularly 
in the 1- to 5-jtm range, is often a problem, jackson and 
Robinson [41 J have shown that, although the response of 
four silicon detectors in a multiband radiometer was nearly 
independent of temperature, the output from three PbS 
detectors, in the 1- to 2.5->im range, decreased linearly with 
increasing detector temperature at a rate of 4 to 5 percent 
per kelvin. For a 1 K change in the temperature of a silicon 
detector, the change in its response was less than 0.25 
percent, whereas the change in the response of a PbS 
detector was 2.5 percent, jackson and Robinson derived a 
correction factor that can be used if the detector tempera- 
ture is known at the time of the measurement. 

D. Estimated Accuracy of Field Reflectance Measurements 

Bearing in mind the measurement problems described 
above, the question is: how accurately can we measure the 
reflectance factor of a surface? The uncertainty in absolute 
reflectance measurement at the NBS is 0.15 percent. Trans- 
ferred to a secondary standard or by reference to pressed 
BaSO» or Halon powder, the uncertainty is probably 0.5 
percent. The calibration of a reference panel with respect to 
the secondary standard should be within an uncertainty of 1 
percent , particularly for near-normal reflectances. Probably 
1.5 percent is the lowest uncertainty we can expect to 
achieve in field reflectance factor measurements. This high 
accuracy may be achievable only in the measurement of 
extended flat surfaces with reflectances greater than 0.3 and 


for small solar zenith angles and a clear and stable atmo- 
sphere. The radiometer design to achieve this accuracy 
should have an extremely well defined field of view and an 
accurately known departure from linearity. In an effort to 
improve the precision of field reflectance measurements, 
the technique described earlier of ratioing direct solar 
irradiance measurements should be studied further. This 
method provides the bidirectional reflectance factor of the 
surface of interest. However, it should be remembered that 
although the ground irradiance in the case of remotely 
sensed data collected by aircraft or satellite sensors is mainly 
directional, there is a hemispherical component that in- 
creases with increasing solar zenith angle, ground reflec- 
tance, and aerosol loading and with decreasing wavelength. 
Because of the many factors that influence hemispheri- 
cal irradiance, it is difficult to relate measurements of 
R(2ir, ♦,) to those of R{9„ <fr). 

£ Description of Measurement Conditions 

The amount of spectral hemispherical-directional or bi- 
directional reflectance factor data characterizing many dif- 
ferent earth-surface features is growing rapidly. Unfor- 
tunately, some of these data, particularly some of those 
collected 10 years ago or earlier, cannot be used reliably 
because the measurement conditions were not well speci- 
fied. Accurate and complete documentation of the instru- 
mentation, measurement technique, and measurement con- 
ditions is mandatory for the data to be transportable. Those 
working at the Laboratory for Applications of Remote Sens- 
ing (LARS) at Purdue University suggest the following infor- 
mation should accompany refleaance data for crops: air 
temperature, relative humidity, wind speed and direction, 
field of view of the spectroradiometer, view angle, latitude, 
longitude, cloud cover, day of year, time of day (in Coordi- 
nated Universal Time), crop species and maturity, row width, 
leaves per plant, row direction, surface conditions, soil 
moisture content, leaf area index, and plants per row. In 
addition, for the reasons mentioned earlier, the irradiance 
conditions should be described as completely as possible. 
Ideally this could be done by making 1) spectral measure- 
ments of the diffuse to direct ratio of Irradiance at the 
ground, 2) an all-sky scan to determine the spectral radi- 
ance anisotropy of the sky, and 3) a measurement of the 
spectral optical depth. Because a rough estimate of the first 
two quantities can be made from a knowledge of the 
spectral optical depth and with the use of an atmospheric 
radiative transfer program, it is possible that a measurement 
of the spectral optical depth may suffice. 

Finally, the reliability of field reflectance factor data would 
gain credibility if the method and frequency of calibration 
of the reference panel were known. Unfortunately, field 
conditions can cause a rapid reduction in the reflectance of 
spray-painted BaSO« panels. In this respect the develop- 
ment by Butner ef al. [42) of washable painted Halon panels 
is noteworthy. 

III. atmospheric Correction 

There was a flurry of activity in atmospheric correction in 
the first year of the Landsat program and during the Skylab 
program, much of which has been reviewed by Slater [8, pp. 
288-317). The activity was caused by concern regarding the 
magnitude of atmospheric path radiance, that is, atmo- 
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spherically scattered light that is added to the radiance of 
the scene. A large component of the path radiance contains 
no information concerning the Kene. An example of its 
magnitude is that, under clear atmospheric conditions, at a 
wavelength of 0.55 pm and for a ground reflectance of 0.1 , 
the atmospheric path radiance at the entrance pupil of a 
satellite sensor is as large as that due to the radiance from 
the ground. 

During the past 10 years, except for a few investigations 
related to the use of models, little work has been done on 
atmospheric correction until recently. This can be attributed 
to the emphasis in the 1970s on the statistical analysis of 
image data, which did not necessitate correction for the 
atmosphere, and to the lack of accuracy of the results of 
those early investigations. The latter was due mainly to 
three reasons: 1) The atmospheric models were inadequate. 
For example, they did not account for the multiple scatter- 
ing and/or the adjacency effect (see later). 2) Because of 
the difficulty of completely characterizing the atmosphere, 
the investigator was forced to make many assumptions, 
which introduced large uncertainties into the result. 3) 
There was no convenient way to check the accuracy of the 
results or use them because of the large uncertainty in the 
in-flight absolute radiometric calibration of the aircraft and 
spacecraft sensors involved. 

A. Correction Methods 

Some different attempts to correct for or allow for the 
atmosphere have been described by Ahern ef a/. (43), 
Aranuvachapun and LeBlond [44], Dana [45], Ooak et a/. 
[46], Dozier and Frew [47], Holyer [48], Kowalik ef a/. [49], 
Kriebel [SO], Lyon ef a/. [51], Monday [52], Ofterman and 
Fraser [53], Price [54], Richardson [55], and Watson and 
Hummer-Miller [56]. Most of these papers are recent and 
are related to the use of models for analyzing satellite- 
acquired image data. Among other methods currently in 
use are the following; 

Gordon ef a/. [57] have described a method for the 
atmospheric correction of Coastal Zone Color Scanner data. 
Their method uses a Monte Carlo atmospheric radiative 
transfer model and an algorithm that includes a ratio of the 
aerosol optical depth at wavelengths of 520 and S50 nm to 
that at 670 nm, where the reflectance of the ocean is 
assumed zero. As Aranuvachapun [58] points out. the accu- 
racy of the algorithm relies mainly on the accuracy of this 
ratio, which is currently not measured by satellite remote 
sensing. The uncertainty of the method in determining 
pigment concentration is stated to be 30 to 40 percent over 
the concentration range 0.08 to 1.5 mg • m'-*. In three 
direct comparisons between ship-measured and satellite- 
determined values of the water radiance. Gordon ef a/, 
claim that the atmospheric correction algorithm had an 
average error of 10 to 15 percent. The method does not 
require any surface measurements at the time of the satel- 
lite overpass. 

Castle ef a/. [59] and Kastner and Slater [60] have de- 
scribed a method making use of ground-based measure- 
ments at White Sands, NM, to determine the in-flight 
absolute radiometric calibration of the Landsat Thematic 
Mapper. The method uses the same type of atmospheric 
data and radiative transfer code that are needed to provide 
the atmospheric correction of satellite-sensor acquired data. 


The determination of the hemispherical-direaional reflec- 
tance factor of an extended, homogeneous surface feature- 
also depends on a knowledge of the exo-atmospheric solar 
irradiance [61], the radiance at. and therefore the absolute 
radiometric calibration of, the sensor, and the angle of 
incidence of the direct solar flux at the surface. The lim- 
itations of the method are that, although it has a low 
uncertainty (about ±3 percent), its accuracy may be com- 
promised by the adjacency effect and the difficulty in 
accounting for the non-lambertian characteristics of an in- 
homogeneous surround. Analyses of the adjacency effect 
are referred to in the next section. 

Table 5 is an extract of Kastner's [62] preliminary results, 
which are relevant to atmospheric correction, of a sensitiv- 
ity analysis of the method for the absolute radiometric 
calibration of the Landsat TM over White Sands, NM. 

Table S Results’ of Sensitivity Analysis for an Atmospheric 
Correction Method that Uses Ground-Based Atmospheric 
Measurements (62) (Note that L is the total radiance at the 
sensor.) 

Particle size distribution. Change in the exponential value, k, in the 
lunge particle size distribution from 2.S to 4 corresponds to a AZ 
of 5 to 10 percent, depending on wavelength, in the 2- to 
O.AS-^m range. The inversion of spectral optical thickness, to 
provide a known > value, can reduce this range by an order 
of magnitude. 

Mie extinction optical thickness. Under good measurement condi- 
tions, reduced Langley plot data can provide with an 
uncertainty such that L can be determined to within 1 percent. 
Refractive index. Calculations show that changing refractive index 
from 1 .5 to 1 .5 - 0.01/ changes t by < I percent. The uncer- 
tainty in L associated with assuming this value is estimated to be 
2 to 5 percent, based on measured index variations at different 
locations worldwide. 

Water vapor. With measured across a band with a solar 

radiometer, the uncertainty in i is <1 percent. 

Surround reflectance. Changing the surface reflectance from 0.35 to 
0.7, which surrounds an infinitesimal target of reflectance 0.7, 
changes I by 2 percent. A Monte Carlo code is needed to 
determine the adjacency effect introduced by any given surface 
spatial distribution of radiances. 

Ground reflectance. For a high-reflectance ( > 0.5) target such as 
White Sands, a given uncertainty in ground reflectance gives rise 
to an equal uncertainty in L. A departure from lambertian 
characteristics of 5 percent gives rise to a change in L of 1 
percent. 

Polarization. Studies (63), [64| need to be extended to determine the 
influence of scene and atmospheric polarization effects. The 
polarization introduced by the sensor itself needs to be known 
from preflight measurements. 

Aerosol vertical distribution; ozone: pressure: time. The aerosol ver- 
tical distribution can be assumed, and the other quantities can be 
determined with sufficient accuracy that their effect on L is 
negligible. 


B. The Adjacency Effect 

The adjacency effect, first analyzed by Pearce [65] using a 
Monte Carlo method and later by Dave [66], describes the 
influence of atmospheric crosstalk in modifying the radi- 
ances of adjacent fields of different radiances. Pearce 
showed that the effect can extend over large distances. For 
example, if the TM were to image two semi-infinite planes 
of reflectances 0.5 and 0.1 at a wavelength of 0.55 )im under 
normal atmospheric conditions, the radiance of the lower 
reflectance area, 1 km from the. edge, would appear to be 
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10 percent more than its asymptotic value. Kaufman (67] has 
recently described the atmospheric effect on the spatial 
resolution of surface imagery. 

Methods for compensating for the adjacency effect have 
recently been described by Tanre ef a/. [68] and Kaufman 
and Fraser [69]. Three experiments have been conducted in 
attempts to verify the adjacency effect. Mekler et al. [70] 
made a laboratory simulation of the atmosphere by sus- 
pending latex spheres in water. The measured effect was 
found to be 20 percent larger than what had been predicted 
by Pearce [65]. Kaufman ef al. [71] flew an aircraft in hazy 
conditions (aerosol optical depth approximately 1.0 at 510 
nm) and demonstrated the existence of the effect at a 
wavelength of 773 nm where the reflectance difference 
between water and vegetation was high. Dyche [72] showed 
in a very clear atmosphere (total optical depth 0.3 at 440 
nm) that the effect may exist but at a level that is difficult to 
detect under such good conditions. 

Methods for the on-board determination of atmospheric 
correction; factors, which make use of multiple ground 
views frorri a pointable sensor, have been suggested by 
Diner and Martonchik [73] and Slater and Martinek [74]. 
These suggestions are preliminary and require further devel- 
opment and testing. 

IV. Sensor System Calibration 

A. Classification of Procedures 

Some of the most commonly used procedures for 
calibrating remote sensing systems are referred to in Fig. 5. 
The two major divisions in the figure are between relative 
and absolute calibration and between the static macro- 
image response and the dynamic micro-image response of 
the system. This section deals with the absolute calibration 
of remote sensing systems, and in this context only the 
static macro-image response will be considered, as is usu- 
ally the case. In Section V of this paper we will return 


briefly to the subject of the effect of sensor micro-image 
response on radiometric fidelity. 

The procedures for the absolute calibration of a remote 
sensing system fall into' the three categories shown in the 
bottom right of Fig. 5. These are described in more detail 
below. 

System Calibration, Preflight Only: In the first category, 
the absolute calibration of the system is made only before 
launch; in flight, the calibration is checked by irradiating 
the focal plane with a radiometrically calibrated source and 
optical system. The drawback to this procedure is that any 
change in the transmission of the image-forming optics of 
the sensor system, caused by the condensation of out- 
gassed contaminants, for example, will be undetected. Fur- 
thermore, the on-board calibration system is assumed to be 
stable through launch and to be unaffected by the vacuum, 
high-energy-particle irradiation, and zero-g environment at 
orbital altitudes — an assumption that can lead to error. 

The TM and the MSS on Landsat-4/5 are examples of 
remote sensing systems calibrated only in this manner. 
Barker ef al. [75] have estimated the preflight absolute 
radiometric calibration of the TM to be about ±10 percent! 
Norwood and Lansing [76] state that it is no better than 
±6.8 percent. Dinguirard and Maisonneuve [77] have esti- 
mated the absolute radiometric calibration of the 
SPOT/HRVs to be ±10 percent. The calibration methods 
employed for TM and SPOT/HRV and their shortcomings 
have been described in detail by Slater [78]. 

Laboratory calibrations, in the solar reflected spectral 
region, are made typically by reference to an integrating 
sphere that in turn is calibrated against a secondary source 
traceable to an irradiance scale at a national standards 
laboratory. In the case of the TM, the irradiance scale at the 
NBS was transferred successively to a primary lamp, a sec- 
ondary lamp, an integrating sphere, the TM focal plane, and 
the TM internal calibrator. The precision with which this 
transfer was done is illustrated in Fig. 6. Starting with the 
primary lamp, Barker and coworkers [75] estimated the 



Fig. S. Classification of radiometric calibration procedures (60). 
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Fig. 7. Arrangement for TM calibration (75]. 


uncertainties associated with the respective transfers to be 
1.3, T.O, 1.0, 2.0, and 3.0 percent. The root sum of the 
squares (RSS) of these uncertainties is ±4.1 percent; we 
emphasize, however, that this is an estimation of the uncer- 
tainty in the precision, not of the absolute accuracy of the 
calibration of the TM internal calibrator. Barker and co- 
workers report that two different measurement techniques 
used to calibrate the sphere gave results that differed inex- 
plicably, usually by 2 to 3 percent but at some wavelengths 
by as much as 7 to 8 percent. These differences coupled 
with the uncertainty in the knowledge of the N8S irradi- 
ance scale (which averages 1 .2 percent over the wavelength 
range 450 to 2500 nm) led Barker and coworkers to estimate 
an uncertainty in the absolute calibration of the TM internal 
calibrator at no less than +6 percent and conservatively at 
±10 percent. 

The layout of the arrangement used for the calibration of 
the TM is shown in Fig. 7. A 1.22-m-diameter sphere was 
required in order to obtain a sufficiently uniform response 
over the 0.41-m-diameter entrance pupil of the TM. A 


similar calibration has been done for the SPOT/HRVs using 
an integrating sphere calibrated by W. A. Hovis of NOAA. 

In-Flight Calibration Using the Son or an On-Board Source: 
In the second category, the sun or an on-board calibrated 
source is used to irradiate the focal plane through the 
image-forming optics. The drawbacks to this approach are 
the unceaainty in the knowledge of a) the irradiance of the 
sun above the atmosphere and b) the output of the 
calibrated source system, for the reasons mentioned earlier. 
Furthermore, in examples of the use of this procedure 
(MSSs 1, 2, and 3 and SPOT), the calibration beam passes 
through only a small portion of the aperture of the system, 
thus not simulating the actual operation of the system. 
When imaging the ground, the system entrance aperture is 
irradiated over its entire area by flux incident over a roughly 
3-sr solid angle. In the imaging mode, much more stray 
light is present in the system and incident on the focal 
plane. If this additional flux level is unknown, it introduces 
a substantial uncertainty into the absolute calibration of the 
system. 
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In 1972 the sun calibrator system on the Landsat I MSS, 
which was intended to provide absolute calibration, 
exhibited a remarkable change in its response. After 21 
orbits, the 0.5- to 0.6-;im band calibration had decayed to 7 
percent of its preflight value [79], The other bands showed 
pronounced but smaller changes. Although Horan et at. [79] 
hypothesized that this was due largely to an unusual amount 
of contamination occurring in preflight tests, the sun 
calibrator system has not been used on any of the Landsat 
MSSs since [80]. 

In-Flight Calibration with Respect to a Ground Surface: 
In the third category of absolute calibration, reference is 
made in flight to a ground area of known spectral radiance. 
If the atmospheric conditions are measured at the time the 
sensor system is imaging the known area, these data can be 
used with an atmospheric radiative transfer program to 
predict the spectral radiance at the entrance pupil of the 
sensor. The main uncertainty in this approach is that of 
determining the refractive index and radial size distribution 
of atmospheric aerosols well enough. The approach is also 
limited to scenes having large, uniform areas of high radi- 
ance. For example, although many bodies of water are of 
sufficient size and uniformity, they are not appropriate for 
calibration purposes because their radiance is too low to 
provide a calibration of sufficient accuracy or to cover 
much of the dynamic range of the sensor. Fortunately, some 
suitable areas do exist, particularly in the arid regions of the 
world, for example at White Sands, NM. References [59], 
[60], and [62] describe the method in detail, and the results 
of a sensitivity analysis of the method are listed in Table 5. 

S. Detector- Based Radiometric Calibration 

The recent work at the NBS on self-calibrated photodi- 
odes is described only briefly here. For more details about 
this important advance in radiometric calibration the inter- 
ested reader is referred to [81 ]-[84], 

The photodiode may be calibrated by either of two 
differently applied biasing procedures, depending on the 
wavelength region of interest. At short wavelengths, a nega- 
tive bias can be applied to remove the recombination 
centers at the Si-SiOj interface at the front of the detector. 
To do this, a temporary contact can be made with the front 
surface using an electrode immersed in a conducting liquid, 
or the surface can be exposed to a corona discharge. Re- 
cently, an inversion layer photodiode has been developed 
that obviates the need for a negative front bias [85]. For 
long wavelengths, a back bias is applied to extend the 
depletion region to a depth beyond which incident flux 
penetrates. The experimental procedure is to irradiate the 
detector with a constant monochromatic flux level and to 
increase the bias voltage until further increase no longer 
gives rise to an increase in output signal. For both the short- 
and the long-wavelength ranges, the internal quantum ef- 
ficiency saturates at a value extremely close to unity, as 
shown in Fig. 8. Thus the maximum increase in signal 
output obtained as a result of biasing can be used to 
determine the internal quantum efficiency of the detector 
without biasing, as it will be used in practice. 

The external quantum efficiency of the photodiode dif- 
fers from the internal quantum efficiency owing to reflec- 
tion losses. These can be reduced to insignificance by 
making use of four photodiodes according to the geometry 
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Fig. 8. Typical photodiode internal quantum efficiency 
without biasing (lower curve) and with biasing (upper curve) 
[83|. 
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Fig. 9. A four-diode arrangement to minimize specular re- 
flection losses (88). 


sketched in Fig. 9. The output signals from the four diodes 
are summed to provide the signal corresponding to a black 
detector of overall quantum efficiency that can be assumed 
to be unity. The second and third diodes specularly reflect 
the specular reflection from diode 1 to diode 4, which 
reflects the light back to diodes 3, 2, and 1. The incident 
flux thereby undergoes seven reflections before leaving the 
system, and if each refleaance is 10 percent, the final 
specularly reflected flux is down to 10“^ of the initial 
incident flux level. The diffuse reflectance loss for a clean 
detector is less than 10“^. The QED four-photodiode 
package manufactured by United Detector Technology has 
been specially designed as a unity quantum efficiency de- 
tector. 

The discussion at the end of this section refers to the use 
of such a defector for the spectroradiometric calibration of 
a multispectral imaging sensor. However, as shown in Fig. 8, 
the unbiased quantum efficiency is wavelength dependent. 
Because of possible changes in the passband position of 
spectral filters during space flights of long duration, it may 
be advisable to use the unity quantum efficiency detector 
to calibrate a spectrally flat, electrically calibrated pyroelec- 
tric radiometer at the 0.1- to 0.2-percent level, and for that 
to be used for the in-orbit calibration, 

Hughes [86] shows that photodiode self-calibration using 
filtered tungsten radiation, in place of laser radiation, can 
provide an absolute spectral response scale from 0.4 to 0.8 
(im with an uncertainty of less than +1 percent. The filters 
used in this work were roughly of 10-nm half-width, which 
is about a factor of 10 less than those of the MSS and TM. 
The possibility of extending this filter approach or a mono- 
chromator approach to provide a means to calibrate field 
and flight spectroradiometer instrumentation with an un- 
certainty in absolute accuracy on the order of +1 percent 
appears then to be close to reality [87], [88]. 
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An important requirement for broad-band absolute radio- 
metric calibration is the use of a spectral response normal- 
ization procedure that is as nearly independent of the 
source function as possible. (The source function in the 
cases considered in this review is the spectral distribution 
of the sun, the calibration lamp, or the radiance of the 
scene used as reference.) Fortunately, Palmer and Tomasko 
[69] have developed a moments normalization method that 
gives exact results if the source function is quadratic and 
results that are superior to those of other methods for 
nonquadratic continuous source functions. Their method 
has the additional advantage over other normalization pro- 
cedures of determining both the height of the rectangular 
function equivalent to the given spectral response and also 
the wavelength limits. Palmer [90] has applied this moments 
normalization procedure to specify the effective band- 
widths of the Landsat 4 and S TMs and MSSs, and Castle ef 
al. [59] have used it in their data reduction for TM calibra- 
tion. 

Calibration in the Laboratory and in Orbit: Similar con- 
cepts have been proposed [60] for sensor calibration in the 
laboratory and in orbit. The main difference between them 
is that an artificial source is used in the laboratory and the 
sun is used in orbit — simply a matter of convenience and 
stability in the laboratory case and of convenience and 
reliability in the latter case. In the laboratory, redundancy is 
not at a premium and our requirements for a source arc 
simply power, spectral flatness, and stability. We do not 
need a standard source although an array of standard N6S 
FEL tungsten halogen lamps could be used. A xenon arc 
stabilized for minimum arc wander and with a highly stable 
power supply and a feedback loop would suffice. 

The source would be used to irradiate a near-lambertian, 
near-unit-reflectance, white-surfaced panel perhaps 1 m by 
O.S m in size in front of the system. (An integrating sphere 
could be used, but it would have to be very large and 
therefore very inefficient, and uniformity checks can some- 
times themselves introduce nonuniformities.) A calibrated 
radiometer whose area-solid angle product and spectral 
response are well defined and accurately measured would 
be used to determine the radiance of the panel in each 
wavelength band of the sensor. The sensor would image 
the panel out of focus, but since it is an extended object, its 
image would have exactly the same irradiance in or out of 
focus. The arrangement is sketched in Fig. 10. 

To avoid problems due to the nonuniform irradiation of 
the panel, the sensor should be rotated to sequentially 
irradiate the focal plane with the image of the same small 
area that is sampled by the radiometer. The reason for a 
large panel is to simulate the viewing conditions from space 



Fig. 10. Factory procedure for absolute spectroradiometric 
calibration [60|. 


in which, depending on the baffle design, significant out- 
of-the-field-of-view stray light could, be incident on the 
image plane to modify the calibration. For this reason it 
would be worthwhile to conduct at least one calibration 
using a white panel several meters in diameter. To check 
for linearity of response, several different irradiance levels 
on the panel should be used. The irradiance level can be 
changed conveniently by inserting heat-resistant neutral 
density filters in -front of the stabilized xenon arc source. 

If the sensor uses spectral filters to define spectral bands, 
the wavelength profiles of the spectral filters should be 
measured by a spectrophotometer using the same f-num- 
ber beam as that for the sensor. If the filters are integral 
with the array, they should measured using a double 
monochromator, again with the sensor f-number, before 
installation in the focal plane. If the sensor is an imaging 
spectrometer, the wavelength calibration can be made by 
reference to spectral line sources or sharp cutting spectral 
filters. For both types of sensor, care should be taken to 
cover the whole wavelength sensitivity range of the detec- 
tors, the off-band suppression being particularly important 
for detectors with the wide spectral response of silicon. 

The procedure proposed here for in-flight calibration is 
similar to the panel method sketched in Fig. 10, but it uses 
the sun as the source. The irradiance over the panel will be 
uniform and known spectrally to better than 1 percent. 
(Several solar measurement programs are currently being 
conducted with this accuracy as a goal.) The calibration 
would be carried out in the few minutes while the spacecraft 
is sun-illuminated but before it images the sun-illuminated 
earth. The absolute radiometer containing a self-calibrated 
detector or an electrically calibrated pyroelectric radiometer 
would now be needed only to check for any deterioration 
in the refleaance of the panel owing to exposure to the 
space environment and short exposures to unattenuated UV 
and other high-energy radiation from the sun and from 
space. In this last respect, the panel would usually be 
stowed in a well-shielded compartment and exposed only 
during actual calibration checks. Also, when deployed, it 
would not interfere with normal operation of the system. If 
a pointing and/or stereo mirror were used, the panel would 
be viewed in an extreme position of the mirror(s). In this 
respect it is fail-safe. Keene [91] has determined that the 
uncertainty of such a calibration procedure could be as low 
as 1.4 percent. 

V. RADiOMeTRic Errors in image Data Collection and 
Processing 

The dynamic micro-image response of remote sensing 
imaging systems and the radiometric effects of digital image 
processing are of importance in any pixel-by-pixel analysis 
procedure. As was indicated in Fig. S, the dynamic image 
response includes 1) the system modulation transfer func- 
tion (MTF), which includes those MTFs due to the imaging 
optics, the effective size of the detectors, and image smear; 
2) the signal-to-noise ratio of the imagery; 3) the electrical 
frequency response of the system; and 4) aliasing. 

A. Radiometric Errors Introduced in Image Data Collection 

Several general investigations have been made into this 
topic, notably a discussion by Norwood [92] on the in- 
fluence of MTF on the errors in radiance measurements of 


SlAItR- RADIOMETRIC CONSIDERATIONS IN REMOTE SENSING 


1007 


different sized ground features, a recent paper by Landgrebe 
and Malaret [93] on the effects of system noise and the 
atmosphere on the probability of classification error, and 
the long-awaited definitive analysis by Park and Schowen- 
gerdt [94] on the MTF of sampled imaging systems that 
accounts for each step in the process of image formation, 
detection, sampling, and reconstruction. Each of these three 
references uses the Landsat MSS or TM as an example. In 
addition, for more detailed system-specific investigations, 
the reader is referred to a special issue of the IEEE Transac- 
tions ON Geoscience and Remote Sensing edited by 
Salomonson [95], in particular the papers by Anuta et al.. 
Bernstein et al., Markham, and Wrigley ef al. These refer to 
analyses of the MSS or TM on Landsat 4 as part of the 
Landsat image data quality assessment program. A summary 
of these papers has been compiled by Schueler and 
Salomonson [%]. The spectral response of sensors has been 
examined by Markham and Barker [97] and Slater [98], who 
show that spectrally dependent striping can occur for sys- 
tems like the MSS due to small departures in the spectral 
responses of channels within a band. 

S. Radiometric Errors Due fo Resampling 

Analysts who use digital imagery from aircraft or 
satellite-borne multispectral imaging systems generally 
purchase it in the processed form, that is, after the image 
has been radiometrically corrected for striping and resam- 
pled to provide systematic distortion correction and reg- 
istration between spectral bands, or to a specified map 
projection. For the past 10 years, a cubic convolution method 
due to Rifman and McKinnon [99] has been used for the 
resampling of Landsat imagery. Park and Schowengerdt 
[100] showed that a parametric implementation of cubic 
convolution is generally superior and can be optimized to 
the frequency content of the image. In later papers, 
Schowengerdt et al. [101], [102] compared the radiometric 
errors introduced for different values of the parameter a, 
the slope of the cubic function at its first crossing. Fig. 11 
shows the mean square radiometric error due to sampling 
and reconstruction as a function of a for digitized aerial 
photographs and Landsat MSS data. The improvement in 
choosing a — —0.5 instead of the usual a — -1 is between 
1 and 3 percent, and in all cases the error associated with 
a — —0.5 or -1 is less than that for bilinear reconstruction. 
Note that the errors are larger for the higher frequency 
content digitized imagery than for the Landsat MSS imag- 
ery, and that in general any resampling error will be worst 
in the vicinity of image edges. 

VI. Conclusion 

An important aspect of remote sensing is the quantitative 
determination of earth surface characteristics by spectro- 
radiometric analyses of imagery collected by aircraft and 
spacecraft sensors. Such analyses are often supported by 
similar ground-based measurements. We have seen from 
this review that uncertainties and errors occur, or noise is 
introduced, in each step in the process of radiative transfer: 
from the sun to the surface, from the surface to a field 
radiometer or an aircraft or spacecraft imaging system, 
through the sensor to the detector, and to a recorded 
radiance level; and finally, errors are introduced in the data 
processing step. 
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Fig. 11. Mean square radiometric error due to sampling and 
reconstruction as a function of a for digitized aerial photo- 
graphs and Landsat MSS data [102|. 

The accuracy of the verification and exploitation of scene 
radiation models depends critically on minimizing these 
various error sources and on completely characterizing the 
surface irradiance conditions, the surface and atmospheric 
conditions, the macro- and micro-image calibration of the 
sensor, and the signal processing techniques employed. 
More specifically, improvements are needed in 1) simulat- 
ing the geometry of space imagery with ground measure- 
ments, 2) further investigating the bidirectional reflectance 
factor method for characterizing surface reflectance to avoid 
sky anisotropy effects, 3) improving ground-based methods 
for characterizing the atmosphere and thereby correcting 
remotely sensed data for its effects. 4) devising a method 
for making atmospheric corrections from orbit, and 5) defin- 
ing the requirements for, and improving the accuracy of, 
the absolute calibration of ground-based, aircraft, and space 
sensors. 

Finally, the magnitude of the error introduced by even 
optimal resampling techniques, such as parametric cubic 
convolution, is alarmingly high and should be of concern to 
all interested in the radiometric fidelity of remote sensing 
data. There is little point in striving for highly accurate 
sensor calibration and atmospheric correction if the radio- 
metric quality of the processed image data is poor. Those 
requiring high radiometric accuracy should use unresam- 
pled data and. where possible, use pixels that are well 
removed from high-contrast edges in order to avoid atmo- 
spheric and sensor-related degradations that are difficult to 
account for because of their dependence on the spatial 
distribution of the surface. 
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ABSTRACT 

Source-based and detector-baaed methods for 
the absolute radiometric calibratlaa of a broad b an d 
field radiometer ace briefly described. Using such 
a radiometer, calibrated by both methods the 
calibration of the Integrating sphere used in the 
pre-flight calibration of the Themadc Mapper was 
re-determined. The results are presented. 

The in-flight calibration of space remote 
sensing instruments is briefly discussed. A method 
is described which uses chs results of ground- 
based reflectance and atmospheric measurements as 
input to a radiative transfer code to predict Che 
radiance at the inscrumenc. A calibrated, 
helicopter-mounted, radiometer is used Co 
determine the radiance levels at intermediate 
altitudes to check the code predlcdoas. Results 
of such measurements for the calibration of the 
Themadc Mapper on Landsat 5 are descdbed 
together with an analysis chat shows the value of 
such measurements. 

Introduction 

To determine ground spectral radiances from 
Che digital counts provided by a mulclspectral 
system in a spacecraft or aircraft or to compare 
data from Che same area collected at several 
different times, we need to sake correcclons for 
Che atmosphere and we need to know the 

relationship between chs radiance at the entrance 
pupil of Che system within each mulclapeccral band 
and Che output number of digital counts. The 
latter is called absolute radiometric callbradoa. 
Its use in quandcatlve remote sensing has been 
discussed by Slater (1, 2). 

In this paper, we descdbe two methods for 
preflighc absolute radiomecrlc callbradoa of a 
mulclspectral system and two methods for inflight 
calibradon. The four methods are fundamentally 
different. One of the preflighc methods is 
referred to as source-based and the other is 
detector-based. One of the inflight methods 
depends on measurements of the spectral 
hemispherical directional reflectance factor of the 
ground, and Che ocher makes use of measurements of 
Che spectral radiance of the ground. In both 

inflight methods the Influence of the intervening 
atmosphere on the measurements has to be taken 
into account. We will provide examples to compare 
Che different methods. 


Source-Based and Detector-Based Pre-Flight 
Calibration 

A. Source-Baaed 

The steps in the source— based calibration of 
sensors such as the Thematic Mapper and the SPOT 
HSV are described briefly below with reference to 
?tguvs 1. The procedure Involves Che transfer of 
calibrations from the source to an Integrating 
sphere and then to the imaging sensor or any 
broadband radiometer. 



Figure 1. The arrangement for calibrating the 
integrating sphere and TM at Santa Barbara Research 
Center. 


1. A 1000-W tungsten halogen light source 
connected to a stabilized power supply is used to 
irradiate a diffuse reflector. The calibration of 
Chs source output and the bidirectional spectral 
reflectance factor of the diffuser are traceable 
to a national standards laboratory. The distance 
from the source to the diffuser (halon in Figure I) 
is accurately measured, so the spectral irradiance 
on Che diffuser is accurately knotrn. 





2. The Integrating sphere, which la calibrated in 
the first transfer, la arranged such that a 
oonochrooator , sampling In 10~nm steps, can 
altematel/ sample and compare the sample to the 
radiance of the diffuser. This la accomplished by 
rotating the mirror shown In Figure 1 through 90°. 
The lenses In front of the monochromator ' are 
adjusted so that a large fraction of the area of 
the sphere exit port and the halon surface are 
sampled by the entrance slit of the monochromator. 
The aperture la adjusted to limit the cone of 
radiation passing through Che monochromator. Thus 
Che comparison is between spectral radiant flux 
levels and the spectral radiance of the sphere 
I*Xa. 


I-Xs 


aU;0°.»5°) 

t 


Bx 



( 1 ) 


where p(X;0°,4S°) is the spectral bidirectional 
reflectance of the halon; S\ la the spectral 
Irradlance of Che standard lamp; and S(X )3 and 
S(.\)i are the signals from the monochromator when 
It Is viewing Che sphere and the standard lamp, 
respectively. This procedure la repeated for 

several sphere output radiance levels achieved by 
Cuming on different numbers of lamps Inside the 
sphere. The highest power level is 1.3 hH. The 
spectral coverage la typically from 0.4 to 2.6 urn 
using SI and PbS detectors Co measure Che oucpuc 
of Che monochromator. Any departures from 
linearity in Che response of the detector system 
oust be accounted for. Thus the spectral radiance 
ac Che exit port of the integrating sphere la known 
at several levels. 


3. The responslvley of the system to be calibrated 
is determined by aligning it so that Its field of 
view is contained within the area of Che exit port 
of Che sphere. The radiance at the exit port of 
Che sphere must not change due Co the proximity of 
Che sysCem to be calibrated. 

4. The ouCpuc digital counts or voltages Fj| 
corresponding co each band, n, are recorded for the 
various radiance levels produced ac Che port of 
Che sphere. 

3. The normalized spectral responses r(X>Q of the 
system are determined by performing a wavelength 
scan using a monochromator chat irradiates the 
entrance pupil of Che system; r(X) is unlcless. 


6. The effective radiance measured by the 
radiometer la given by t ;^Lxt(X)QdX, where the 
integration Is performed over all wavelengths co 
which the system responds. The relationship 
between this quantity and for each band 

compares Che total radiance at the entrance pupil 
of Che sysCem co che oucpuc volcage or digital 
counts. This constitutes che absolute radiometric 
calibration of Che system. 

3. Detector-Based 

Method 1. The steps in one method of detector- 

based calibration of a broadband radiometer are 
described as follows: 


detector of accurately known spectral 
responslvley, described by Zalewskl and Ouda (3). 

2. A wavelength scan of che monochromator then 
provides a calibration of Its output radiant flux, 
♦(X). 

3. The radiometer to be calibrated is inserted 
into che beam In place of the QED, where again the 
radiometer detector area must totally encompass 
the beam. The output volcage 7(X) is recorded as 
a function of wavelength. 

4. Applying Palmer's (4) moments normalization 
method co che ratio 7(X)/*<1) provides a measure of 
the effective bandwidth and responslvley Bg of che 
radiometer. 


3. A peak normallzaclon method Is used to 
determine che effective field of view, which is 
used CO calculate the solid angle, 3, subtended by 
che detector ac Che entrance pupil. A simple 
measurement of Che diameter of che entrance pupil 
provides che area A. 

6. The absolute calibration of che cadlomacer Is 
given by the ' relation between lx* spectral 

radiance ac che entrance pupil, and V(X), the 
voltage out of che radiometer according co 


I-X “ 


RgAB 


( 2 ) 


Method 2. An alternative detector-based procedure 
is as followst 


1. The end of a cube with two small circular 
apertures in It Is attached to Che front of a QED 
detector as shown in Figure 2. The geometry of 
Che cube defines Che A3 product mentioned earlier. 
The filter la removed from the radiometer to be 
calibrated and placed over che cube. 
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Figure 2. The four-diode QEO-100 with throughput 
adaptor. 

2. This QED spectroradlomecer is used co measure 
che radiance of a panel Irradiated with a stable, 
hlgh-power whlte-Ught source. The voltage, V, 
from che detector is recorded. 

3. The spectral transmittance of the filter is 

measured with a spectrophotometer and che moments 
normalization method is applied to che 
product t(X) R(X), where R(X) is che known spectral 
response of che QED. This gives che normalized 

responslvley of che filtered QED detector Rg. 

4. Applying Equation (2) yields che spectral 
radiance of che panel. ’ 


1. A very stable, hlgh-power tungsten halogen 
lamp is used co irradiate che entrance slit of a 
monochromator. The output beam from Che 
monochromator is adjusted Co be smaller chan che 
area of a United Detector Technology QEP-100 


5. The filter is replaced in the radiometer Co be 
calibrated, and the instrument is directed at the 
same angle co Che panel as was che QED 
spectroradlomecer. 
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TABLE 1. PERCENTAGE DIFFERENCES IN CALIBRATIONS* ueTYIEEN: 



OSC DETECTOR- 

OSC SOURCE- 

SBRC MONOCHROMATOR (SM) 


BASEO (OOB) AND 

BASED (OSB) AND 

AND SBRC SOURCE-BASED 


SBRC SOURCE- 

SBRC SOURCE- 

(SSB), COMPARED AGAINST 


BASED (SSB) 

BASED (SSB) 

INTEGRATING SPHERE 

WAVELENGTH SS8-OOB 

in nm TOT"*’®® 

SSB-OSB...« 

— SSB“*’®® 

SM-SSB 

— S35|-*’®® 


•• 

•• 

*• 

440.3 

-10.8 

-9.4 

•►4.1 

525.4 

-2.6 

-0.7 

+2.9 

605.4 

-2.9 

-O.S 

+2.S 

862.1 

-2.0 

•H>.8 

+2.1 

780.7 


-►1.4 

+0.6 

861.7 

•►8.7 

■►2.5 

-1.0 


* All tlM comparisons ara for tfw sama radiomatar, axcapt for tha last column whara tha 
radlometar and monochromator wars cailbratad against tha sama source and halon diffusar 
and than both maasurad tha Intagrating sphara output. Tha valuas listed raprasant tha 
dlffaraneas batwaan tha raadings from tha taro Instrumants whan nnaasuring tha sphara output. 
** Low signal to noisa ratio conditions. 

*** Oatactor haatar not working. 


6. The relatioaahlp becween the voltage fcoo the 
radiometer and tha spectral radiance of the panel 
provides the absolute calibration of tha 
radiometer. 

C. Comparison o£ Results 

Table I compares source-based and detector-based 
calibrations using a Elltec spectropolarlmeter 
described by Castle (3) with a 3* field of view. 
Some of Che calibrations were conducted at the 
Optical Sciences Center (OSC) of the University of 
Arizona; others were conducted at Santa Barbara 
Research Center (SBRC), a division of Hughes 
Aircraft Company by Vlcman (6). The Integrating 
sphere and associated equipment (see Figure I) used 
In the pref light absolute radiometric calibration 
of Thematic Mappers 4 and 3 are at SBRC. 

The comparisons In Table 1 show chat good 
agreement exists for most of the visible and near 
Infrared wavelengths used. At shorter wavelengths 
Che uncertainty In che measurements increases 
because of chelr low signal to noise ratio. This 
is che reason che first detector-based method 
described can be used co calibrate only 
spectroradlomecers chat use large detectors. Tha 
irradlance across che exit silt of che 
monochromator Is low, so co obtain a large signal 
a long silt is needed. 

Inflight calibration by reference to the Earth's 
surface. 

A. Reflectance-Based Methods 

The Inflight calibration methods considered here 
utilize simultaneous measurements of che earth's 
surface and atmosphere and then predict radiance 
levels at the sensor. The relation of chose 
radiance levels co che digital counts In che Image 
of che measured ground area provides che absolute 
calibration. The advantage of such measurements 
is that they are made under the exact conditions 
che system Is _ operating under and therefore all 
sources of transmission loss and stray light ace 
accounted for. 


Gordon et al. (7) have described a method for the 
ataospherle correction of Coastal Zona Color 
Scanner (CZCS) data Chat can. In conjunction with 
ocean reflectance measurements, provide a 
calibration of che sensor. Their method uses a 
Monte Carlo ataospherle radiative transfer model 
and an algorithm chat Includes a ratio of che 
aerosol optical depth at wavelengths of 320 am and 
330 OB to chat at 670 am, where the ocean 
reflectance Is assumed zero. As Aranuvachapon (8) 
paints out, the accuracy of che algorithm relies 
mainly on che accuracy of this ratio, which Is not 
currently measured by sacelUce remote sensing. 
The uncertainty of che method In determining 
pigment concentration Is stated co be 30 co 40 
percent over the concentration range 0.08 co 1.3 
mg m'’. In three direct comparisons becween ship- 
measured and sacelllce-decermlned values of water 
radiance, Gordon and associates claim chat che 
atmospheric correction algorithm had an average 
error of 10 to 13 percent. 

Koepke (9), using a radiative transfer program chat 
accounts for multiple scattering In conjunction 
with ground reflectance data and atmospheric data, 
performed an Inflight calibration of Meceosat I for 
use In radiation budget studies. Atmospheric water 
vapor was calculated using data from the nearest 
radiosonde, and che climatological values of ozone 
and oxygen were assumed. Ground level barometric 
pressure defined che Rayleigh optical depth, and 
aerosol optical depths were derived from visibility 
data. Four test sices were used: rough ocean, che 
savanna In tiamlbia, pascureland In northern 
Germany, and freshly fallen snow. Sldlrecclonal 

reflectance data were collected for these sites 
over a period of several months. Histograms of 
these areas typically filled less chan three 
columns, and often one column contained 60 per 
cent of che pixels. We note, however, that 
Meceosat quantizes data co 6 bits, so for low 
reflectances a three column spread represents a 
large variation In reflectance. nevertheless by 
using four targets and pointing che sensor Into 
deep space, which fixed che offset point on the 
councs-versus-radlance curve, Koepke estimated che 
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uncertaiacy to be on the order of 6 percent. Tble 
resourceful use of colacerel data provided a 
valuable calibration of Heteosat 1. Koeplce's 
estimate of an uncertainty of 6 per cent probably 
represents a lower limit on the uncertainty of the 
method because of Its dependence on reflectance 
data that were not collected simultaneously with 
the Heteosat Imagery and the use of v i si bil ity 
estimates that provide only approximate values for 
optical depths. 

During the past three years, we have used this 
technique several times at White Sands, Hew Mexico, 
but we have taken simultaneous ground and 
atmospheric measureaents as the Thematic Mapper 
acquired Images of the' ground area. This work has 
been described by Kastner and Slater (10) and 
Castle et al. (11). Originally the measurements of 
the ground were only of Its spectral hemispherical 
directional reflectance. Recently ground level and 
helicopter measurements of the spectral radiance 
of the surface have been made. Results of such 
measurements are shown alongside the reflectance— 
baaed results for October 23, 1984 In Figure 3. 
The conditions at White Sands and the reflectance- 
baaed results are listed In Table 2. In Figure 3, 
the Fre and I.C. bars represent, respectively, the 
spectral radiance values of the ground determined, 
from the unprocessed digital counts of the Image 
using the preflight calibration (see the 
comparisons In Table 1) and the Internal calibrator 
updates provided by Barker (12). The Code bar 
refers to the results of the reflectance-based 
method, which makes use of a radiative transfer 
code developed by Herman and Browning (13). The 
Hell bar represents the radiance at 3,000 m as 
measured by a spectroradlometer. The H— Hell bar 
represents modified helicopter results. In this 
case the ground reflectance used In the code was 
modified to give a radiance value at 3,000 a chat 
matched chat measured by the helicopter, and the 
results In the bar chart are the corresponding 
radiances at the cop of the atmosphere. The 
Rayleigh bars represent the radiances at the cop 
of the atmosphere for a Rayleigh atmosphere. Note 
chat Che Rayleigh results are within 2 per cent of 
the result for the no atmosphere case for TM band 
1 and within I percent for bands 2, 1, and 4. 

The results in Figure 3 show anomalously high 
radiance levels In bands 2 and 4 based on preflight 
and I.C. calibration data. That the two 
Independent reflectance and radiance-based methods 
and Che Rayleigh case give results chat are 
closely bunched Is strong evidence chat the 
preflight and I. C. results are mlsrepresencatlve 
of Che calibration of TM on chat dace. 
Furthermore, results on ocher daces (July 3, 1984, 
May 24, 1985, and August 28, 1985) fall to show 
this anomaly although there was no change In our 
measurement methodology, reference panel 

calibration, and data reduction methods. We have 
no reason to question Che preflight or Che Internal 
calibrator gain and offset values themselves; we 
believe Che anomalies are more likely Che result 
of a source of radiance In the atmosphere 
surrounding the spacecraft, perhaps due to 
luminescence or Ionization phenomena. High energy 
solar UV or X-ray Irradiation could be one cause 
and Che passage of Che spacecraft through a gas 
cloud created by an earlier orbit adjustment has 
been suggested as another. This last hypothesis is 
based upon the observed deceleration of Che 
spacecraft as It passes through the region In space 
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TM calibration at White Sands 
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Figure 3. Comparisicn of results for TM cattbratian at 
White Sands, see text. 

where the orbit adjustment jets had been 

previously activated. These decelerations occur 
for two Co ctoee weeks after an orbit adjustment. 

The fact chat such anomalies were not detected by 
Che internal calibrator Is an Important reason for 
using ground-based callbraclon techniques whenever 
possible. The following discussion shows chat one 
such technique Is easy to Implement given a ground 
reflectance of about 0.5 in each band. 

B. Radiance-Based Methods 

Figure 3 shows chat Che radiances as measured at 
3000 m above White Sands are about . Che same as 
chose measured by the Thematic Mapper. This 
result prompted a more detailed examination of Che 
variation of radiance with altitude as a function 
of wavelength, reflectance, and visibility. 

Starting with the atmospheric conditions of October 
28, 1984, at White Sands, and with a solar zenith 
angle of 35 deg, radiance values were calculated 
for altitudes between ground level (1196 m at 
White Sands) and 21,000 m (the operational altitude 
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Figure 4. (^»rcefttage differences betvnoen the radiances at the top of the atmosphere and at 
inteimatfiate aititudes for the White Sands condtions of Octoher 28, 1984, tee Tabie 2. 


above sea level of the HA8A 0-2 or E&-2 aircraft) 
and for above the acaoaphere. Graphs ware then 
plotted of the percentage difference (the 

differences between the radiance above the 

ataoaphere and at an Intarmedlace altitude divided 
by the radiance above the acaoaphere, the whole 
quantity expressed as a percentage) against 
altitude above sea level. These are shown In 
Figure 4, for the first four solar reflective 
TH bands. As anclclpaced, the results for 
reflectances less chon 0.1 show large percentage 
differences, but what Is Intriguing is the sign 
change In the percentage difference for 

reflectances between 0.2S and 1.0 for bands 1, 2, 
and 3 and between 0.1 and 0.23 for band 4. These 
results show chat, for a reflectance of 0.3 and an 
altitude of 3000 a ASL, Che percentage differences 
are -^2.3, -1.0, -1.0, and -2.0 for bands 1, 2, 3, 
and 4 respectively. This Indicates chat for the 
White Sands conditions of 28 October, 1984, 

helicopter radiance values at an altitude of 
3000 n are close to chose for space and require 
only a small correcdon, which can be made 
accurately. 

With this promising result, Che Investigation was 
extended co determine how general this condition 
might be. Visibilities of 10 and 23 Icm were 
considered for a reflectance of 0.3 at all 
wavelengths. The Input values are listed In Table 
3, Che refractive Index was 1.34 - 0.011 In all 
cases and Junge v values of 2.9 and 2.6 were used 
for Che vlslblUcles of 10 and 23 Icm respectively. 
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The output plots of percentage difference, as 
defined earlier, against altitude above sea level 
are shown as Figures 3 and 6 together with Che 
Rayleigh atmosphere curve for a wavelength of 0.49 
pm. The positive difference for 1 ■ 0.49 pm la 
due CO Che fact that the large path radiance term 
contributed more co the total radiance chan is 
lost by the accenuatloa of the ground reflected 
term. The ocher curves show increasing negative 
differences for altitudes approaching ground level, 
owing CO Che heavy aerosol loading In Che lower 
atmosphere. The bowing of Che X ■ 0.37 pm and 
0.66 pm curves at Intermediate altitudes la due co 
reduced pech radiance from Rayleigh scattering and 
Che presence of ozone absorption at chose 
altitudes. Note chat for an altitude of 3000 m 
and a vislblllcy of 23 Icm Che percentage 
differences are +3.0, -1.4, -0.9, -0.8, -2.3, and 
-1.6 for wavelengths of 0.49, 0.37, 0.66, 0.84, 
1.67, and 2.22 pm respectively. Furthermore the 
percentage differences are less chan twice these 
values for a vislblllcy of 10 Icm, as can be seen 
from Che remaricable slmUarlcy becween Figures 3 
and 6. 



ALTmjOe ABOVE 3EA UEVEUKxIOOOft. or 304.8m) 

Figure S. Percentage differences between the radiances 
at the top of the atmospere and at intenmediate 
altitudes for paO.3 and a visibility of 10 km. 
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Figure 6. Percentage differences betvtieen the radiances 
at the top of the aemospere and at intermediate 
attitudes for and a vaibdity of 23 km. 


ConcluBlon 

The callbracloa of tanging and nooianglng broadband 
radloaeeera by dececcor-based and aoucce-baaed 
aechoda needa to be continued to check and refine 
Che accuracy of Che farmer. Although 

cheoreeically more accurate, In practice the 
meaauremenca aaaoeiaeed with Che detector-baaed 
method are aoaetiaea difficult becauae of low 
signal Co noiae racioa and ill-defined radiometer 
geomeeriea. Further work needs co be done to 
determine the practical liaicaclons of dececcor- 
based calibrations accuracy for broadband 
radiometers. 

The Inflight calibrations of several satellite 
mul clspeccral sensors have shown substantial 
changes with time. The CZCS has shown what is 
believed co be a monotonic decrease In response 
for all wavelengths, with the decrease greater for 
shorter wavelengths (14). The same is true for TM 
for the solar reflective bands with the following 
two exceptions: The mid-IR bands show oscilliCory 

changes (13), and our results for October 28, 1984, 
show anomalously high responses in bands 2 and 4. 

The inflight calibration of satellite sensors by 
reference to ground reflectance measurements can 
be conveniently checked by the use of low altitude 
(3000 m) helicopter radiance measurements. Without 
correction the latter are within 3 percent of the 
satellite values for the solar reflective bands of 
TM, and with correction they are limited by the 
uncertainty in the absolute calibration of the 
helicopter radiometer. They do require a ground 
reflectance of about 0.3 but they are insensitive 
to atmospheric conditions for a wide range of 
visibilities above about 13 km. Compared co 
satellite calibration by reference co high altitude 
aircraft measurements, Che helicopter method 
enjoys the advantage of lower cost, longer dwell 
time over Che area of interest, and Che capability 
of averaging over a 30 m diameter pixel with a 1 
degree field of view instrument from 3000 m 
altitude. 

Again, because of Che many difficulties chat beset 
Che inflight absolute radiometric calibration of 
sacellice systems, we are not proposing chat one 
method be used co Che exclusion of others. 
Whenever possible, effort should be made Co 


coordinate different but simultaneous measurements 
over, a given site with instruments chat have been 
intercompared in the laboratory. Only by such a 
concerted effort will we be able to identify 
systematic errors and reduce them, and be 
confident of the uncertainties chat we ascribe co 
our calibrations. 
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ABSTRACT 

Variations teported in the in-flight absoluts 
radiometric calibration of the Coastal Zone Color 
Scanner (C2CS) and the Theaaclc Mapper (TM) will 
be reviewed. At short wevelangchs, chase sensors 
show a gradual reduction in response while In the 
mld-IR Che TM shows osdllacory varladona. One 
see of asasureoenes made at White Sands, Maw 
Mexico shows anoaalous results in TM bands Z and 
4. 

Tha results of a reflaeeanea-based and a 
radiance-based callbraelaa oeehod at White Sends, 
Mew Mexico will be described. An analysis of tbs 
radiance-based oeehod is presented showing the 
value of such oeasureaenes from helicopter 
altitudes. 


Incroductioo 

It is now well established chat the absolute 
radiometric calibration of sacelUce remote sensing 
systems changes with time. The need to \oiow the 
absolute calibration and methods to obtain it have 
been described elsewhere (ref. 1,2). 

In this paper, we review Che vartadons in the 
in-flight absolute calibration of Che Coastal Zone 
Color Scanner (CZCS) and Che Themaclc Mapper' (TM). 
We outline two methods used for In-fllghc absolute 
calibration by reference co ground sices, one which 
uClUaes ground spectral reflectance measurements 
and Che other which utilizes ground spectral 
radiance measurements. An analysis of Che 
radiance-based method is presented showing Che 
value of such measurements from helicopter 
altitudes. 

Variations in Radiometric Calibration 

Ue briefly describe here che variations noted 
in che calibration of che CZCS and TM. With 
respect co che CZCS, Ic Is worth noting chat Gordon 
(ref. 3), from an error budget analysis, showed 
chat a 3 percent error In sensor calibration and a 
10 percent uncertainty in che exoacmospheric solar 


irradlance can comhina co give rise co errors 
exceeding ZOO percent in che retrieval of water 
reflaecance values. Although che 10 percent 
uncertainty should perhaps be 1-2 percent, Prohlleh 
(ref. 4), and water raflaccaneea are low making 
than relatively hard co measure accurately, che 
point remains chat large errors in reflectance 
measuremancs can occur due co sensor calibration 
errors. Gordon showed chat, by waking 
slmulcaneous measuremancs of aerosol optical 
chlckaasses and water radiance with che satelUte 
osasurement of cecal radiance, che arror in 
dacermlnlng water refleetanca can be reduced by an 
order of magnicuda. 

In che summer of 1979, ViolUer (ref, 3) 
compared CZCS derived water reflectance values 
with measuremancs mads nesr the surface. He found 
chat the satellite value was less chan chat of che 
surface by 3.5 percent at 443 nm and greater chan 
chat at the surface by 6.3 and 12 percent at 320 
and 350 am respectively. VlolUer comments chat 
Cordon reported corresponding values of -2, >2 and 
7 percent at che Tlfceench CZCS MET Meeting but 
used a different set of calibration constants and 
exoacmospheric irradlance values. 

Gordon ec al. (ref. 6) have described a method 
for che atmospheric correction of (CZCS) data chat 
can. In conjunction with ocean . reflectance 
measurements, provide a calibration of che sensor. 
Their method uses a Monte Carlo atmospheric 
radiative transfer model and an algorlchm chat 
includes a ratio of the aerosol optical depth at 
wavelengths of 320 nm and 330 am co chac at 
670 am, where che ocean reflectance is assumed 
zero. As Aranuvachapun (ref. 7) points out, che 
accuracy of che algorithm relies mainly on che 
accuracy of this ratio, which is not currently 
measured by sacelUce remote sensing. The 
uncertainty of che method in determining pigment 
concentration is stated co be 30 co 40 percent 
over che concencraclon range 0.08 co 1.5 mg m~*. 
In three direct comparlsone between ship-measured 
and sacelUce— determined values of water radiance, 
Gordon and associates claim chac che atmospheric 
correccian algorlchm had an average error of 10 co 
15 percent. 

3y computing che radiance at che sacelUce 
from iuiown water radiance values, Gordon ec. al 
(ref. S) proceeded co force agreement between che 
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Figure 1. The decay factor with orbit number for the first band 
of the CZCS, see text (Cardan ec al., ref. 8). 


3eaaor~reeorded radiance an<i the coopueed radiance 
by adjuaclng the seneor caiibracion. Their reeulca 
are shown In Figure I for band 1 at CZCS. The 
error bars repreaenc the error Induced by the 
uncercalncy In the pigment coaceneraeion of the 
water. 

Hovla at al. (ref. 9) uoed the second oa-boerd 
calibration lamp, used only occaaionally, to check 
Che calibration of CZCS after It had been In orbic 
for S 1/2 years. They found no meaeurable change 
in the calibration, but it la Important to tuite 
chat Che calibration lamp checka only the 
calibration of the focal plane detectors and 
associated electronics, as Is also the case In the 
Hulclapectral Scanner System and TM. In compering 
CZCS data with chose from a calibrated radiometer 
In a high altitude aircraft. Howls and aasoelates 
found a 23 percent degradation In Che blue band, 
centered at 443 am, afcer CZCS had been In orbit 
for four years and seven months. The longer 
wavelength bands exhlblced progressively smaller 
degradations as shown In Figure 2. Because of che 
calibration lamp results, the degradation is 

attributed to a reduction In che reflectances of 
che telescope optics and che scanning mirror. 



Figure 2. The percentage degradation in the response of the 
CZCS with time for the first four spectral bfuids (Hovis et alv 
ref. 9). 


oortheast Pacific subtropical gyre. The squares 
are data points generated using che radiometric 
sensitivity decay correction models of Gordon ec 
al. (ref. 3) with random noise, and che dashed 
lines are chair models. Again the tendencies are 
as noted by che authors deed earlier, although It 
Is Interesting to note the linear fit used by 
Hueller compared to the quadratic fit used by 
Gordon ec al. (ref. 3). 
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Figure 3. Corrections for che radiometric sensitivitv of rhn 
CZCS with time for the first three spectral bands (Mueller, ref. 
10 ) 


Frequently recorded values for Che 
calibration gains In che solar reflective bands of 
TM on Landsac 4 have been accumulated by Lae and 
Barker (ref. 11). Figure 4 compares che pre-launch 
calibrations with che values from che Internal 
calibrator (IC). Again che results, as for che 
CZCS lamp-based calibration mendoned earlier, only 
apply CO changes In che focal-plane fllcers and 
dsceccors and che associated electronics. The 
graphs In Figure 4 are divided into three sections. 
The left secclon Includes che pre-launch 
calibration conducted under thermal vacuum 
condlcloos. The middle section shows Che results 
of In-orblc IC calibrations during che so-called 
Scrounge Image processing era and che third secclon 
shows the results during che Thematic mapper Image 
Processing System (TIPS) era. The gains, G, are 
given from che equation F • (G x L) 0 where ? is 
the observed IC pulse, L Is the IC spectral 
radiance and 0 Is che offset. 


Mueller (ref. 10) has also reported on the 
change in che calibration gains of Che first three 
channels of che Mlmbus-7 CZCS during its first four 
years of operation. Figure 3 shows che gains for 
the three channels centered at 443, 320 and 3S0 am 
plotted against year and orbit number. The solid 
lines are linear regression models and represent 
che average fits over 20 replications with 
Independent samples of Gaussian random noise. The 
circles represent decay coefficients calculated for 
pixels In che central water masses of che 


The results In Figure 4 show decreases In che 
apparent gains for band 1 amounting to 7 percenc 
and 13 percent after 300 and 600 days from launch, 
respectively. Sands 2 and 3 show a similar 

pattern buc with a smaller loss In gain. Band 4 
deports from che trend exhibited by bands 1,2 and 
3 In showing a gain Increase for che last 
caiibracion. Sands 3 and 7 show osclllacory 
changes In gain which are out of phase to start 
with buc Chen appear to be In phase. The 
amplitude of che oscillations Is about 7 percenc. 
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Figure 4. The apparent gain change with time for the solar reflective bands of the TM on Landsat 4, The 
gain values are in digital numbers per mW/on'* sr urn (Lee and Barker, ref. 11). 

It cnusc be repeated chat these results are only superlaipoaed on them should be the changes that 

Indicative of changes in the filter-detector- occur in the reflectances of the optical 

electronics part of the system and that components of the telescope and scan mirror. 












In-fUght eallbraelon by ’ reference to eha F.arrh'a 
surface. 


In addition to the work described earlier on 
Che In-flight callbracioa of the CZCS by reference 
to know water reflectances, two ocher reflectance- 
based calibrations are described here, chose of 
Heceosac and TM. These are followed by a 
description of a radiance- baaed calibration oachod. 

A. Sefleccance-Sased Methods 

Eoepka (12) using a radlaelva transfer program 
chat accounts for multiple scaccerlng In 

conjunction with ground reflectance data and 
atmospheric data, performed an In-flighc 

callbracloo of Meceosac 1 for use In radiation 
budget scudlea. Acmoapherlc water vapor was 
calculated using data from the nearest radiosonde, 
and Che climatological values of osone oxygen 
were assumed. Ground level barometric pressure 
defined the Bayleigh optical depth, and aerosol 
optical depths ware derived from vislbillcy data. 
Four test sites were used: rough ocean, the 

savanna In Hamlbta, paacuraland In northern 

Germany, and freshly fallen snow. Sidiracciooal 

reflectance data ware collected for these sices 
over a period of several months. Histograms of 
these areas typically filled lass chan three 
columns, and often one column contained 60 par 
cent of Che pixels. Ua note, however, chat 
Meceosac quantlaes data to 6 bits, so for low 
reflectances a three column spread represents a 
large variation In reflectance. Meverchelass by 
using four targets and pointing the sensor into 
deep space, which fixed the offset point on the 
counts-versus-radianca curve, Hoepka ascimated the 
uncertainty to be on the order of 6 percent. This 
resourceful use of colateral data provided a 
valuable callbracloo of Meceosac 1. Koepka’s 
eadmacs of an uncertainty of 6 par cent probably 
represents a lower limit on the uncertainty of the 
method because of its dapaodenee on reflectance 
data chat were ooc collected simultaneously with 
Che Meceosac imagery and the use of visibility 
estimates chat provide only approximate values for 
optical depths. 

During Che past three years, we have used 
this technique several times at Uhice Sands, Mew 
Mexico, but we have taken simultaneous ground and 
atmospheric measurements as Che Thematic Mapper 
acquired Images of Che ground area. This work has 
been described by Kascner and Slater (13) and 
Castle ec al. (14). Originally the maasucemencs of 
the ground were only of its spectral hemispherical 
directional reflectance. Hecencly ground level and 
helicopter measurements of the spectral radiance 
of Che surface have been made. Hesulcs of such 
measurements are shown alongside the reflectance- 
based results for October 23, 1984 In Figure S. 

The conditions at White Sands and the reflectance- 
based . results are listed In Table 1. In Figure 3, 
Che Pre and IC bars represent, respectively, the 
spectral radiance values of the ground determined 
from Che unprocessed digital counts of the Image 
using Che pre-flight calibration and Che internal 
calibrator updates provided by Barker (11). The 
Code bar refers to the results of the reflectance- 
based method, which makes use of a radiative 
transfer code developed by Herman and Browning 
(IS). The Hell bar represents the radiance at 

3,000 m as measured by a spectroradlometer. The 
M-Hell bar represents modified helicopter results. 
In this case Che ground reflectance used in the 
code was modified to give a radiance value at 
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TABLE 1 : White Sands CaUbratian of TM on October 28, 1964. 
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Figure S. Comparision of results for TM calibration at White 
Sands, see text. 

3,000 m chat matched chat measured by Che 

hellcopcer, and Che results in che bar chart are 
Che corresponding radiances at che cop of che 
atmosphere. The Bayleigh bars represent che 
radiances at Che cop of che atmosphere for a 
Bayleigh atmosphere. Mote chat Che Bayleigh 

results are within 2 per cent of che result for 
che no atmosphere case for TM band I and within 1 
percent for bands 2, 3, and 4. 


The results in Figure S show anomalously high 
radiance levels in bands 2 and 4 based on pre- 
flight and IC calibration data. That che two 
independent reflectance and radiance-based methods 
and che Rayleigh case give results chat are 
closely bunched is strong evidence chat che pre- 
flight and IC results are miarepresentatlve of che 
calibration of TM on chat date. Furthermore, 
results on other dates (July 8, 1984, May 24, 1983, 
and August 28, 1983) fall to show chls anomaly 
although there was no change in our measuremenc 
methodology, reference panel calibration, and data 
reduction methods. We have no reason co question 
che pre-flight or che internal calibrator gain and 
offset values themselves; we believe che anomalies 
are more likely the result of a source of radiance 
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In Che acoosphara surrounding eha spaeeeralt, 
perhaps due co luolnaseaaea or lonixacloa 

phaaomana. dlgh energy solar U7 or X-*ray 

Irradiaclon could bo oao causa and eha paaaaga of 
eha spaeeerafe ehrough a gas cloud craaead by an 
earlier orbie adjuaeoene has bean suggasead as 
anoehar. This last hypoehaals Is based upon eha 
observed daealeraclon of eha spaeecrafc as it 
passes ehrough eha raglan ta space where eha orbie 
adjusenane jees had been previously acelvacsd. 
Thase deealeradons occur for evo eo ehrea weeks 
after an orbie adjusCaane< 

Tha fact Chat such anooallas wars not 
detected by eha Ineemal callbraeor Is an ioportane 
reason for using ground-based callbracion 
cechnlquss whenever possible. The following 
discussion shows ehae one such ceehnlqua Is easy to 
ImplenenC given a ground rafleceance of about 0.5 
In each band. 


and 3 and between 0.1 and 0.25 for band 4. These 
results show chat, for a reflectance of 0.5 and an 
aldtude of 3000 m ASL, eha percentage differences 
are '^2.5, -1.0, -1.0, and -2.0 for bands I, 2, 3, 
and 4 respeedvaly. This Indicates chat for eha 
UU.ee Sands condlclons of 28 October, 1984, 
helicopter radiance values at an aldtude of 
3000 a are close eo chose for space and require 
only a soall eorreedon, wUch can be made 
accurately. 

With this promising result, the Invesdgaclon 
was extended co daceraine how general this 
condition lalght be. Vlslbllldes of 10 and 23 km 
were considered for a refleccancs of 0.5 at all 
wavalangche. The Input values are listed In Table 
2, Che refractive Index was 1.54 - 0.011 In all 
cases and Junge v values of 2.9 and 2.6 were used 
for eha vlslbllldes of 10 and 23 km respeedvaly. 


3. Eladlance-Baaed Hechods 

Figure 5 shows ehae the radiances as aeasured 
at 3000 a above UUCe Sands are about Che same as 
chose aeasured by the Themadc Mapper. This 
result prompted a aore detailed examlnadon of Che 
variation of radiance wleh aldtude as a function 
of wavelength, reflectance, and vlslblUcy. 


The output plots of percentage difference, as 
defined earlier, against aldtude above sea level 
are shown as Figures 7 and 8 together with Che 
Baylelgh atmosphere curve for a wavelength of 0.49 
ua. The posldve difference for 1 • 0.49 um Is 
due CO Che face that Che large path radiance term 
contributed aore to Che total radiance chan Is 
lost by Che accenuadon of Che ground reflected 


Starting with Che atmospheric conditions of 
October 28, 1984, at UUce Sands, and with a solar 
zenith angle of 35 deg, radiance values wars 
calculated for aldcudes between ground level 
(1196 a at White Sands) and 21,000 a (the 
operadonal aldtude above sea level of the NASA 
U-2 or EE- 2 aircraft) and for above Che 
atmosphere. Graphs were Chen plotted of Che 

percentage difference (the differences between the 
radiance above Che atmosphere and at an 
Intermediate altitude divided by the radiance above 
Che atmosphere, Che whole quandcy expressed as a 
percentage) against altitude above sea level. 
These are shown In Figure 6, for the first four 
solar reflective TM bands. As anticipated, Che 
results for reflectances less chan 0.1 show large 
percentage differences, but what Is intriguing Is 
Che sign change in the percentage difference for 
reflectances between 0.25 and 1.0 for bands 1, 2, 
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differences for aldcudes approaching ground level, 
owing CO the heavy aerosol loading in the lower 
atmosphere. The bowing of the \ • 0.57 urn and 
0.66 urn curves at intermediate altitudes is due to 
reduced path radiance from Rayleigh scattering and 


PERCENT fUFFERENCE PERCENT DIFFERENCE 


ViaBUTY > lOkffl 



Figure 7. Percentage differences betvxeen the radiances at the 
top of the atmosphere and at intermediate altitudes for a^.S 
and a visibility of 10 bn. 



Figure 3. Percentage differences between the radiances at the 
top of the atmosphere and at intermediste altitudes for saO.S 
and a visibility of 23 bn. 

Che presence of ozone absorpclon ac chose 
alclcudes. Note chac for an alclcude of 3000 a 
and a vlsltilllcy of 23 km Che percentage 
differences are +3.0, -1.4, -0.9, -0.3, -2.3, and 
-1.6 for uavelengchs of 0.49, 0.37, 0.66, 0.34, 
1.67, and 2.22 urn respecclvely. Furchernore Che 
percencage differences are leas chan cwlce chese 
values for a vlalbillcy of 10 ka, as can be seen 
from Che reoarkable slailaclcy becween Figures 7 
and 3. 

Conclusion 

The In-flighc calibrations of several 

satelUee muldspeccral sensors have shown 

substantial changes with cine. The CZCS has shown 
uhac Is believed co be a nonoconic decrease in 
response for all wavelengths, with Che decrease 
greater for shorter wavelengths. The sane is true 
for TM for che solar reflective bands with che 
following cwo exceptions: The nid-IE bands show 

osclllicory changes, and our results for October 
28, 1984, show anonalously high responses In bands 
2 and 4. 

That che degradation In response of che CZCS 
was noc detected by che on-board calibration 
system and our anomalous results for TM were noc 
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detected by che IC Is evidence of che desirablllcy 
of using sons calibration procedure chac involves 
ima g in g che Earth's surface and applying 
acmoapberlc corrections. 

The in-flight calibration of satellite sensors 
by reference co ground reflectance measurements 
can be conveniently checked by che use of low 
aldcude (3000 m) helicopter radiance neasuremenca. 
Without correction the latter are within t3 percent 
of che secallica values for the solar raflecdve 
kawda of TM, and with correction they are liolced 
by the uneercalncy in che absolute callbradon of 
the helicopter radiometer. They do require a 
ground reflectance of about 0.3 but they are 
insansldve co atmospheric condidons for a wide 
range of vlsibilides above about 13 km. Compared 
CO satellite calibration by reference co high 
aldcude aircraft measurements, the helicopter 
method enjoys che advantage of lower cost, longer 
dwell time over che area of interest, and che 
capability of averaging over a 30 m diameter pixel 
with a 1 degree field of view Inacrusienc from 
3000 m aldcude. 

Again, because of che many diffleuldes chac 
beset che in-flighc absolute radiomscrie 
calibration of satellite systems, we are noc 
proposing chac one method be used co che exclusion 
of others. Whenever possible, effort should be 
made to coordlnace different but simtilcaneous 
measurements over a given site with Instruments 
chac have been incereompared in the laboratory. 
Only by such a concerted effort will wo be able to 
identify syscamade errors and reduce them, and be 
confident of che uncertainties chat we ascribe co 
our calibradons. 
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Survey of Multispectral Imagiiig Systems for Earth Observations 
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Fifty-siz mulUspectral imagiiig systems are described in terms of their instantaneous fields of view, spectral bands, 
fields of view, and number of quantization levels. Iliese systems have been used during the past decade for earth 
resources studies from aircraft or spacecraft, or are currently in the proposal or design and development stage. 



Introduction 

An increasing interest woridwide in re- 
search and applications studies in remote 
sensing, together with the development of 
new types of optical detectors and sys- 
tems, has led to a recent proliferation of 
multispectral imaging systems^ for use on 
aircraft and spacecraft Most of these sys- 
tems are experimental, their purpose 
being to determine the optimum selection 
of spectral bands and the S|)atial and ra- 
diometric resolutions suitable for different 
research and applications studies. 

The purpose of this paper is to provide 
a listing of these systems as a convenient 
reference for remote sensing investigators 
and systems engineers. For space systems, 
the listing includes, to the author’s knowl- 
edge, all past and present systems and 
those planned for the future. With re- 
spect to airborne systems the listing is not 
exhaustive. For example, military systems, 
such as forward-looking infrared systems, 
are not included, nor are multispectral 
video cameras, which are becoming in- 

‘ For the purpose of this survey, the term " multispec- 
trai imaging system" indudes any optical system that 
''scans the object or image surface electromechanically or 
electronically in more than one wavelength interval It 
excludes photographic film cameras and nonimage- 
forming spectroradiometers. 

0E3sevier Science Publishing Co.. Inc., 1985 
.52 Vanderbilt Ave., New York. NY 10017 


creasingiy popular for aerial remote sens- 
ing in the visible and near infrared. Also, 
when improvements have been made to a 
particular system, only the latest model is 
described. 

In addition to the references cited in 
the listing, there are several general ref- 
erences worthy of note. Due (1982) 
describes in considerable detail many 
military and civilian systems, both past 
and present, among them the AVHKR-2, 
M^S, Landsat MSS, S-192, TM, and 
VISSR, which are listed in the appendix. 
(’The acronyms are defined in the appen- 
dix.) Doyle (1984) and the Landsat Data 
Users Notes (1982) describe the status of 
space remote sensing plans woridwide for 
scanners, film cameras, and microwave 
sensors. Slater (1980) describes in detail 
many past and present space systems. 
Holmes (1984) discusses several advanced 
sensor systems. 

Many authors have suggested spectral 
bands for earth observation purposes, 
among them the foUowing: Begni (1982) 
describes the rationale for the choice of 
the SPOT HRV bands; Cihlar et ai. ( 1984) 
describe the choice for Radarsat; and Col- 
vocoresses (1979) describes the choice 
for Mapsat. Stoner and Baiungardner 
(1980) discuss the location of bands most 

0034-1257/8.5/.S3.30 
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suitable for soil classifications. Siegrist and 
Schnetzler (1980) discuss the location of 
bands for geological applications. 
Wiersma and Landgrebe (1979) discuss 
the location of bands suitable for vegeta- 
tioix monitoring. Tucker (1978). lists the 
desirable bands for vegetation monitoring 
in order of priority. 

Vane et aL (1982) summarize the rec- 
onunendations in some of the publish- 
ed literature and in the Multispectral 
Resource Sampler (MRS) workshop, de- 
scribed in the ORI Report (1979), regard- 
ing the spatial and spectral requirements 
for future sensor systems. Cox (1982) 
summarizes discussions of the Multispec- 
tral Imaging Science Working Croup con- 
vened by NASA, which included a 
Botanical Sciences Team, a Geographic 
Sciences Team (including panels on land 
use and land cover, geomorphology, and 
cartography), a Hydrologic Sciences 
Team, a Workshop for Lithologic Map- 
ping, an Imaging Science Panel, and an 
Information Science Panel 

System Characteristics 

The appendix lists the multispectral 
imaging systems and satellites for earth 
observations as well as spectral band rec- 
ommendations of science working groups 
concerned with the definition of future 
systems. The listing includes references to 
publications that describe these systems 
or recommendations in detail or to the 
individuals who supplied the information. 

The following comments relate to the 
Appendix: 

1. Thirteen coimtries and the European 
Space Agency (ESA), have been or 
are involved in the development of 
one or more of the systems. 
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2. Of the 56 systems listed, 11 are no 
longer in use, 23 are currently in 
use, and 22 are proposed or in the 
design and development phase. 

3. Of the 56 systems, 32 are spacecraft 
systems, of which three are primari- 
ly meteorological systems, 23 are 
aircraft systems, and 1 has been used 
both as an aircraft and a spacecraft 
system. 

4. Of the 56 systems, 24 are mechani- 
cs scanners, 20 are pushbroom 
scarmers using linear detector arrays, 
5 are television systems (four of 
which use vidicons and one of which 
uses an area detector array); 5 are 
imaging spectrometers using area 
detector arrays; and 2 are spin-scan 
systems. 

Discussion 

All numerical data listed in the appen- 
dix should be taken as nominal or pro- 
jected values. Depending on the system 
and the source reference, however, they 
may be accurately measmed quantities. 
In some cases the systems are in the 
planning stage and may change substan- 
tially before fabrication. The reader should 
refer to the listed references to obtain 
further details and to form an opinion on 
the reliability of the data quoted in the 
references. 

This survey shows that many different 
kinds of multispectral imaging systems 
are now available for research and appli- 
cation purposes in remote sensing, and 
that the number, availability, and im- 
proved performance of linear and area 
arrays of detectors, using mainly CCD 
and CID readouts, has encouraged several 
groups to explore their application as 
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aircraft and spacecraft remote sensing 
systems. 

The advantages of using these detec- 
tors are that a scanning mirror is no longer 
required, the system can be pointed for 
stereo imaging or for more frequent imag- 
ing of a given area, and the integration 
time and therefore the signal-to-noise ratio 
and radiometric sensitivity are higher than 
for the equivalent mechanical scanner. 

The shortcomings of silicon CCOs and 
CEDs are that their long-waveiength cut- 
off is I.l /xm and, unless carefully selected, 
they exhibit large nonuniformities in re- 
sponse and noise properties. However, 
improvements in the yields of high per- 
formance CCDs and CIDs and the pro- 


duction of IRC CDs are subjects of 
substantial research and development ac- 
tivity. So, although few linear or area 
CCD or CED systems are currently in 
use, we can expect them to rival the 
worldwide data output of mechanical 
scaimers, both in volume and spectral 
coverage, by the end of the decade. 

In addition to those persons mentioned 
in the references as having provided 
infimnation on certain systems, I should 
like to thank K. J. Ando, W. L. Barnes, 
C. D. Mason, and M. S. Maxwell for their 
inputs and comments, and acknowledge 
financial support from MASA grant NAG 
5-196, 


Appendix: List of Multispectral Imaging Systems and Satellites for Earth Observa- 
tions, and Spectral Band Recommendations of Science Working Groups Concerned 
with the Definition of Future Systems 


AADS-1260 Daedalus airborne multispectral scaimer. The AADS series of mech- 
anical scanners is commercially available for airborne use; refer- 
enced here is the 11-band system. IFOV 2.5 mr, field of view 86”, 
8-bit quantization. Spectral bands: 0.38-0.42, 0.42-0.45, 0.45-0.50, 
0.50-0.55, 0.55-0.60, 0.60-0.65, 0.65-0.69, 0.70-0.79, 0.80-0.89, 
0.90-1.1, 10.4-12.5 jam (Daedalus Enterprises, Inc.). 

AADS-1268 This is an extensively modified version of AADS-1260 (above) for 
use as a Thematic Mapper (see TM) simulator. The IFOV can be 
selected to be 2.5 or 1.25 mr with field of view 86” or 43°, 
respectively. Spectral bands: 0.42-0.45, 0.45-0.52, 0.52-0.60, 
0.605-0.625, 0.63-0.69, 0.695-0.75, 0.76-0.90, 0.91-1.05, 1.55-1.75, 
2.08-2.35, 8.5-13 pm. Note that the thermal band is wider than that 
of the Thematic Mapper and that there are four additional bands 
(Daedalus Enterprises, Inc.). 


AEROS-A Advanced Earth Resources Observational Satellite-A, propos^ for 

1986 launch financed by Space .America. The payload consists of 
two 4-band linear array sensors, one looking 22.5° forward and 5° 
starboard, the other 22.5” aft and 5” port. The 0.63-0.69 ,um band 
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in the former and the 0.6-0.7 /xm band in the latter will have a 45-m 
IFOV and a 149-km swath; the other bands will have an 80-m IFOV 
and a 1654an swath. The second sensor will image the same ground 
swath as the first but one day later, providing stereo data with a 
base-to-height ratio of 0.83 and a height resolution of 104 m. The 
orbit is that of Landsats 1-3 (see Landsat). 8-bit quantization, 
18-Mb/s bandwidth using X-band. Spectral bands: 0.45-0.52, 
0.52-0.6, 0.63-0.69, 0.78-0.9 fim for the forward looking sensor and 
0.433-0.53, 0.54-0.56, 0.6-0. 7, 0.7-0.8 /xm for the aft looking sensor 
(Lowndes, 1983). 

Airborne Imagmg Spectrometer. This system was recently built and 
tested at the Jet Propulsion Laboratory for measurements in the 
1.2-2.4 fim range. IFOV 2 mr, field of view 3.6°, 8-bit quantization. 
A 32 x 32 detector area array is stepped four times during data 
coUection to provide 128 bands in the in-track direction, each 
having 10-tun bandwidth (Vane et aL, 1983). 

Australian Multispectral Scaimer System. (The system is actually 
known as MSS, but AMSS is used here to avoid confusion with 
Landsat MSS.) A mechanical scanner for aircraft use under develop- 
ment in Australia, primarily for ocean color studies. IFOV 4 mr, 
field of view 40°, 8- or 10-bit quantization. Spectral bands: 0.345- 
0.395, 0.43-0.45, 0.45-0.47, 0.47-4).48, 0.49-0.51, 0.51-0.53, 0.53- 
0.55, 0.68-0.69, 0.975-1.025 (Cartwright, 1983). 

Advanced Solid-state Array Spectroradiometer. Modification of an 
existing aircraft system by NASA Johnson Space Center to include a 
512 X 32 CID array in the focal plane of a grating spectrometer that 
will provide 32 contiguous spectral bands 14-nm wide covering the 
range 0.40 to 0.85 fim or 0.50 to 0.95 /xm. The IFOV will be 0.34 or 
0.854 nu, the field of view 10° or 25°, and there will be 12-bit 
quantization with variable gain and offset before quantization. The 
system is designed to map chlorophyll, determine chlorophyll pro- 
files, and study ocean surface effects. The spectral bands for this 
. system are to be determined (Stewart and Buntzen, 1982). 

Advanced Very High Resolution Radiometer. A 5-band mechanical 
scaimer system on a NOAA meteorological satellite. IFOV 1.1 km, 
field of view ±56°, inclination 98.9°, orbit altitude 862 km, repeat 
coverage interval 12 h. To map cloud cover, surface, and cloud top 
temperatures, land-water boimdaries, and snow and ice coverage. 
Spectral bands: 0.58-0.68, 0.725-1.0, 3.55-3.93, 10.2-11.3, 11.5-12.5 
fim (Due, 1982; Schneider and McGinnis, 1982). 
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AVIRIS 


Bhaskara-2 


BST 


CAESAR 


CBML/CSmO 


CZCS 


Airborne Visible-InfraRed Imaging Spectrometer. A mechanical 
scanner imaging spectrometer capable of operating in as many as 
224 spectral bands in the 0.4 to 2.4 /xm region. Under construction 
at the Jet Propulsion Laboratory as an airborne research instrument 
for development of new data analysis methodologies for use with 
imaging spectrometer data. IFOV 1.0 mr, field of view 30°, 10-bit 
quantization, data rate 15 Mb/s (Vane et aL, 1984). 

An Indian system launched on 20 November 1981 that consisted of 
a 2-band TV camera system with 1-km IFOV. Eccentric orijit 
varying between 514 and 557-km altitude at 50° inclination. Spec- 
tral bands: 0.54-0.65 and 0.75-0.85 /xm (Doyle, 1984). 

Botanical Science Team, a group of 21 participants convened by 
NASA in 1982 to recommend research activities and parameters for 
future sensors. Recommended spectral bands: 0.54-0.56, 0.65-0.67, 
0.84-0.86, 1.52-1.74, 2.08-2.32, 3.5-3.9, 10.5-12.5 ;xm (Cox, 1982). 

CCD Airborne Experimental Scanner for Applications in Remote 
sensing. A modular CCD pushbroom system in the development 
phase at the National Aerospace Laboratory NHL and at the In- 
stitute of Applied Physics TPD, both in the Netherlands. IFOV 0.3 
mr, field of view 25°, 12-bit quantization. Capability of 54° forward 
look for canopy and water studies. Spectral bands: 0.405-0.415, 
0.51-0.53, 0.535-0.565, 0.555-0.575, 0.62-0.64, 0.655-0.685, 
0.675-0.695, 0.77-0.80, 0.825-0.875 /xm (Bunnik, 1984b). 

A modular mechanical 15-channel scanner for aircraft use being 
developed by Carr Boyd Minerals Limited and the Commonwealth 
Scientific and Industriai Research Organization, Australia. IFOV 5 
mr, field of view 90°. Can be used in an active mode using a laser. 
Prototype band selections: 0.45-0.55, 0.55-0.6, 0.65-0.7, 0.83-0.87, 
0.93-0.97, 1.55-1.65, 1.98-2.08, 2.13-2.23, 2.3-2.4, 7.75-8.75, 
8.75-9.35, 9.25-10.25, 10.75-11.25, 11.75-12.25 /xm (Lyon, 1983). 

Coastal 2^ne Color Scanner. A 6-band mechanical scanner orbited 
on Nimbus 7 on October 1978 in a 955-km circular, sun-synchronous 
orbit. IFOV 800-m, swath 1600-km, 8-bit quantization. Modifica- 
tions to the system are being studied with the anticipation of 
orbiting a modified 8-band unit on the TIROS-N series (NOAA-H 
and/or I). The Tiros system may have a 980-m IFOV. For determin- 
ing global oceanic chlorophyll concentration and diffuse attenuation 
coefficients and sediment transport in coastal zones. Spectral bands 
for Tiros and Nimbus systems: 0.43-0.45, 0.48-0.50, 0.51-0.53, 
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0.55-0.57, 0.58-0.60, 0.66-0.68, 0.75-0.79, 0.84-0.89 /im (Barnes, 
1983). 

FILE Feature Identification and Location Experiment. A 2-band CCD 

system flown on the second Shuttle. The system consisted of two 
filtered CCD area array cameras, each array consisting of 100 X 100 
detectors. The IFOV was 1 x 0.75 km and the field of view 17° by 
23°. The orbit was circular at an altitude of 255 km with a period of 
1 hr 29 m 40 s. Data were collected on 12-14 November 1981. 
Spectral bands: 0.64-0.66 and 0.84-0.^ jam (Silverton et aL, 1982). 

FLD Fraunhofer Line Discriminator. A variable single-band, aircraft, 

mechanical scanner for detecting solar-stimulated luminescence. Op- 
erates in Fraunhofer absorption lines, typically at 656.3 nm, Na 
at 589.0 am. Mg at 518.4 nm, at 486.1 nm, Ca at 422.7 nm. 
Uses Fabry-Perot filter with < 0.1-nm halfwidth. IFOV 1°, field of 
view 37° (Hemphill, 1981; Slater, 1980; Watson and Theisen, 1981). 

FLI Fluorescence Line Imager. A progranunable spectroscopic imaging 

system with five two-dimensional CCD array sensors operating in a 
pushbroom mode being designed by Moniteq Ltd. of Canada for 
delivery in January 1985 to the Canadian Department of Fisheries 
and Oceans and to be flown in an aircraft of the Canada Centre for 
Remote Sensing. The FLI can operate over the range 0.42 to 0.81 
fim with 1.9-nm bandwidths and a choice from 256 spectral chan- 
nels. Eight spectral bands can be used in two configurations, one for 
general chlorophyll mapping (A) and one for studying fluorescence 
line properties (B). IFOV 1 mr. Each sensor covers about 15° to 
give a total field of view of 73°, 12-bit quantization. Typical spectral 
bands: 0.433-0.453 (A&B), 0.510-0.530 (A), 0.540-0.560 (A&B), 
0.633-0.644 (A), 0.66-0.671 (A&B), 0.671-0.682 (A&B), 0.682-0.685 
(B), 0.685-0.690 (B), 0.708-0.714 (A&B), 0.746-0.754 (A), 0.772- 
0.790 (B) /un (Edel, 1982; Gower, 1982). 


Fragment-2 8-band mechanical scanner orbited as part of the Meteor spacecraft 
payload (see Meteor), by the USSR. IFOV 80 m, swath width 85 
km. Spectral bands: 0.40-0.80, 0.50-0.60, 0.60-0.70, 0.70-0.80, 
0.70-1.1, 1.2-1.3, 1.5-1.8, 2.1-2.4 fim (Avanesov et al., 1981). 


HCM Heat Capacity Mapper. A 2-band mechanical scanner used on a U-2 

in support of the HCMR (see below). IFOV 2.8 mr, field of view 
90°. Spectral bands: 0.51-0.89 and 10.44-12.52 jam (Price, 1984). 
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Heat Capacity Mapping Radiometer. A 2-band mechanical scanner 
launched 26 April 1978 into nearly sim-synchronous 620-km circular 
orbit at 97.6® inclination, repeat coverage interval 16 days. IFOV 
0.5 km in visible, 0.6 km in m, swath width 716 km. Equatorial 
crossing times 0230 and 1330. Spectral bands: 0.55-1.1 and 10.5-12.5 
fim (Price, 1978; Slater, 1980). 

Hi^ Resolution Visible. Two HRV systems, which are CCD 
pushbroom scanners, are under construction for the French Centre 
National d’Etudes Spatiales in association with Belgium and Sweden. 
The HRV systems will be launched on SPOT (see SPOT). Each 
HRV will operate in 3 spectral bands with an IFOV of 20 m or 1 
panchromatic band with an IFOV of 10 m. Each covers a 60-km 
swath with 84Dit quantization. Pointing mirror to ± 27° cross-track 
can provide stereo or repeat coverage interval at least every 5 days. 
Data rate 48 Mb/s for both HRV systems, transmitted in X-band. 
Spectral bands: 0.50-0.59, 0.61-0.68, 0.79-0.89 /im. Panchromatic 
bluid: 0.51-0.73 ;im (Midan et aL, 1982). 

Indian Resources Satellite, see LISS. 

Japanese Elarth Resources Satellite. Proposed by the Japanese Sci- 
ence and Technology Agency for a Launch in 1989-90, it is a linear 
array stereo camera with 30-m IFOV and four spectral bands 
between 0.51 and 1.1 (im. No further details are available (Doyle, 
1984). 

A NASA program and five associated satellite instrument payloads. 
The first payload launched on 23 July 1972 (then called the Earth 
Resources Technology Satellite and now Landsat-1), included a 
multispectral Return Beam Vidicon (RBV) system, a four-band 
Multispectral Scaimer System (MSS), and a Data Collection System 
(DCS). (See the Landsat Data Users Handbook, 1979.) Landsat-2, 
launched on 22 January 1975, contained instruments with the same 
nominal characteristics as Landsat-1. Landsat-3, launched on 5 
March 1978, contained a two-camera panchromatic RBV system and 
an MSS with a fifth band in the thermal infrared. Landsats 1, 2, and 
3 were inserted into a circular, near-polar orbit at a nominal altitude 
of 919 km, with an inclination of 99.11°, a period of 103.27 min, a 
repeat coverage interval of 18 days, and an equatorial crossing time 
of nominally 9:30 a.m. (the crossing times could change by an hour 
dming the system's lifetime). Landsat-4, launched in July 1982, 
contains an MSS with similar characteristics to MSS 1 and 2. It also 
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contains the 7-band Thematic Mapper but no RBV camera or DCS. 
The circular orbit is at an altitude of 705 km, inclination is 98.2°, 
repeat coverage interval is 16 days, and the equatorial crossing time 
is 0945 ±15. Landsat-5, identical to Landsat-4, was launched in 
March 1984 (Landsat Data Users Notes, 1982). 

Linear Array Pushbroom Radiometer. A pushbroom airborne system 
built and currently under test at NASA Goddard Space Flight 
Center. IFOV 1.54 mr, field of view 42.5°, 10-bit quantization. 
There are four lens-detector modules in the system. A filter wheel 
allows for the selection of 1 of 6 bandpasses per module in flight. 
The system is pointable through ± 45° both in-track and across-track. 
Spectral bands: 0.40-0.42, 0.45-0.52, 0.49-0.51, 0.52-0.60, 0.54- 
0.56, 0.63-0.69, 0.67-0.69, 0.71-0.73, 0.73-0.75, 0.76-0.90, 0.78- 
0.80, 0.84-0.86, 0.88-0.90, 0.93-0.95 inn (Irons et al., 1982). 

Linear Imaging Sell-Scanned sensors, three of which are to be flown 
on the Indian Resources Satellite (IRS-IA), planned launch date 
1985-86. LISS-1 will have four spectral bands in the visible and near 
IR, an IFOV of 73 m, and a swath of 148 km. Two LlSS-2 systems 
will cover adjacent swaths 74-km wide with IFOVs of 37 m in four 
spectral bands. IRS-IA will be placed into a 904-km sun-synchro- 
nous orbit with a 99° inclination, and repeat coverage will be 
obtained every 22 days. S-band and X-band telemetry. Plans call for 
the sensors on IRS-IB to provide IFOVs of 15 to 20 m and those on 
IRS-2 to include middle and thermal ER bands (Doyle 1984). 

Modular Multiband Scanner. An 11-band airborne mechanical 
scanner built by Bendix Corp. IFOV 2.5 mr, field of view 100° ± 10°. 
Spectral bands: 0.38-0.44, 0.44-0.49, 0.49-0.54, 0.54-0.58, 0.58- 
0.62, 0.62-0.66, 0.66-0.70, 0.70-0.74, 0.77-0.86, 0.97-1.06 nm, and 
a band of filter-selectable width at approximately 11 /xm (Due, 
1982). 

University of Michigan Scanner. A 12-band airborne mechanical 
scanner used extensively during the Com Blight Watch experiment 
in 1971. The scanner has been the property of the Environmental 
Research Institute of Michigan since 1973 and has been modified 
several times since 1971. IFOV 2 mr, field of view 90°. Included 17 
bands, only 12 of which could be used at one time. Spectral bands; 
0.33-0.35, 0.46-0.49, 0.48-0.51, 0.50-0.54, 0.52-0.57, 0.54-0.60, 
0.58-0.65, 0.61-0.71, 0.72-0.92, 1.0-1.4, 1.5-1.8, 2.0-2.6, 9.3-11.7 
jam (Lowe, 1975; Landgrebe, 1981). 
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An orbital mapping system, proposed by the U.S. Geological Survey, 
that can operate in a 3-band mviltispectral or a panchromatic 
stereoscopic CCD pushbroom mode. Orbital characteristics are those 
of landsats 1, 2, and 3 (see Landsat). Imaging proposed in an 
epipolar orbit to facilitate stereo compilation. Mixed resolution, the 
finest being an IFOV of 10 m, 8-bit quantization. Swath width 60 to 
185 km depending on imaging mode for bandwidth not to exceed 48 
Mb/s. Nadir and ±23*^ viewing to obtain 20-m contoxus. Spectral 
bands: 0.47-0.57, 0.57-0.70, 0.76-1.05 fim (Colvocoresses, 1979, 
1982; Snyder, 1982). 

Multispectral Electro-optical Imaging Spectrometer. An airborne 
pushbroom CCD system built in 1982 by Macdonald Dettwiler and 
Associates Ltd., for the Canada Centre for Remote Sensing. A 
modular 54)and system with future expansion to 8 bands including 
IR. Output IFOV (after resampling) 0.71 mr, input IFOV 0.53 mr, 
13-jam detectors. Fidd of view 50® for 24.6-mm lens, 84jit quantiza- 
tion. Real time resampling to provide band-to-band pixel registration 
to ± 1/8 pixel Real time dark cturent subtraction and relative 
radiometric calibration. Spectral bands: 0.45-0.52, 0.50-0.60, 0.52- 
0.60, 0.60-0.70, 0.63-0.69, 0.70-0.80, 0.76-0.90, 0.80-1.1 nm 
(Hrebenyk, 1982; Macdonald Dettwiler and Associates). 

Multispectral Electro Optical Stereo Scanner, a system under study 
by Deutsche Forschungs- und Versuchsanstalt ftir Luft und Raum- 
fahrt (DFVLR), West Germany, and India. Five CCD linear arrays 
are located across the focal plane of a Hasselblad camera using a 
Zeiss Biogon lens of 61-mm focal length. This arrangement will 
provide two forward, two aft, and one vertical view of the ground. 
Depending on the arrays and lens chosen, the IFOVs will be 67 or 
150 m and the swath 180 or 350 km. The system will be placed in 
orbit by the Indian Rohini launch vehicle in April 1985 at an 
altitude of 400 km and at a 44° inclination. Telemetry bandwidth 
12.3 Mb/s. Spectral bands: 0.5-0.6, 0.6-0.7, 0.7-0.8 /xm (Doyle 
1984). 

Multispectral Electronic Self-Soanning Radiometer. Two 4-band 
CCD pushbroom MESSRs are being built by Nippon Electric Co., 
Inc., as part of the payload for the Japanese Marine Observation 
Satellite I to be launched in 1986. Sun-synchronous circular orbit, 
altitude 909 km, inclination 99.1°. Time of descending node be- 
tween 10 and 11 a.m. IFOV 50 m, swath for both MESSRs 
combined is 200 km, 6-bit quantization, data rate 8 Mb/s for each 
system. Spectral bands: 0.51-0.59, 0.61-0.69, 0.72-0.80, 0.8— 1.1 /xm 
(Matsumoto, 1981; Ishizawa et al., 1980). 
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An earth observation and meteorology spacecraft orbited by the 
USSR in June 1980. Eccentric orbit between altitudes of 589 and 
678 km widi a 90** inclination. Payload contained 5 multispectral 
systems: MSU-M and MSU-S (together referred to as RTVK), MSU-E 
and MSU-SK (together referr^ to as BK-E), and Fragment-2. Refer 
to each for further details (Sidorenko, 1981). 

A European Space Agency geostationary meteorological spin-scan 
radiometer located at longitude 0°. East-to-west scan provided by 
spacecraft rotation round an axis parallel to earth’s axis, south-to- 
north scan provided by small stepping mirror. First system launched 
in November 1977; second system launched in May 1981; four 
additional launches piaimed for the 1980s. From 36,000-km altitude 
the IFOV is 2.5 km in the visible and near-IR band and 5 km in the 
IR bands. 6-bit quantization in spectral bands: 0.5-1.0 and 5.7-7. 1 
jxm, and 8-bit in the 10.5-12.5 /im band (Jouan, 1982). 

A Multispectral Linear Array system for the Space Transportation 
System under development at NASA Goddard Space Flight Center 
for shuttle flights in the late 1980s. The system has six spectral 
bands in the 0.45-1.65 ixm region: 0.46-0.48, 0.54-0.56, 0.66-0.68, 
0.87-0.89, 1J28-1.25, 1.64-1.66 jam. IFOV is 15 m in the first four 
bands, 30 m in the last two. Swath width is 30 km. The system 
provides 8-bit quantization with a data rate of 48 Mb/s. Along-track 
pointing achieved using a steerable mirror, across-track pointing (of 
±45”) achieved by roiling the shuttle (Barnes and Salomonson, 
1983). 

Modular Opto-electronic Multispectral Scaimer. A CCD pushbroom 
system built by Messerchmitt-Bolkow-Blohm for the Shuttle Pallet 
Satellite (SPAS-01) and flown in July 1983 on the Space Transporta- 
tion System-7 (STS-7). A two-band system expandable to four bands 
for subsequent flights, IFOV 20 m, field of view 140 km from 
shuttle altitude of 296 km; 7-bit quantization. Spectral bands; 
0.575-0.625 and 0.825-0.975 /xm (Hofmann et al., 1982). 

Marine Observation Satellite, see MESSR. 

Multispectral Resource Sampler. A workshop involving 70 attendees 
representing user interests met in 1979 to suggest requirements for a 
new orbital earth resources research system. Spectral bands: 0.36- 
0.40, 0.40-0.42, 0.41-0.43, 0.43-0.45, 0.44-0.46, 0.48-0.50, 0.53- 
0.55, 0.54-0.56, 0.58-0.62, 0.59-0.61, 0.66-0.68, 0.67-0.69, 0.73- 
0.78, 0.74-0.76, 0.76-0.78, 0.78-0.82, 0.84-0.90 urn (ORI Report, 
1979). 
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Multispectral Scanner System. Five mechanical scanner systems 
have been orbited as part of the Landsat program (see Landsat). 
EFOV 76 m for MSS 1 to 3 and 80 m for MSS 4 and 3, swath width 
185 km, 6-bit quantization, data rate 15 Mb/s. Spectral bands: 
0.5-0.6, 0.6-0.7, 0.7-0.8, 0.8-1.1 nm (Due, 1982; Landsat Data 
Users Handbook, 1979; Lansing and Cline, 1975; Norwood et aL, 
1972; Markham and Barker, 1983; Slater, 1979). 

A 4-band television system on the Meteor spacecraft (see Meteor). 
IFOV 1.7 X 1 km, field of view 1930 km. Spectral bands: 0.50-0.60, 
0.60-0.70, 0.70-0.80, 0.80-1.0 /im (Selivanov and Tuchin, 1981). 

A 2-band television system on the Meteor spacecraft (see Meteor). 
IFOV 240 m, field of view 1380 km. Spectral bands: 0.50-0.70 and 
0.70-1.1 [im (Selivanov and Tu chin, 1981). 


A 4-band conical image scaimer on the Meteor spacecraft (see 
Meteor). IFOV 170 m, field of view 600 km. Spectral bands: 
0.50-0.60, 0.60-0.70, 0.70-0.80, 0.80-1.0 jitm (Selivanov et aL, 
1981). 


A 3-band system using electronically scanned arrays on the Meteor 
spacecraft (see Meteor). IFOV 30 m, field of view 30 km. Spectral 
bands: 0.50-0.70, 0,70-0.80, 0.80-1.0 jam (Selivanov et al., 1981). 

Multispectral Fast Area Digitizing Scaimer, a system comprising 
four filtered cameras each using 128 x 128 area arrays as detectors, 
under development for aircraft use in New Zealand. The system 
operates in the snapshot mode with an IFOV of 0.35®, field of view 
of 35®, Obit quantization. Spectral bands: 0.5-0.6, 0.6-0.7, 0.7-0.8, 
0.8-1.1 /im (Hodgson et al., 1981). 

An 8band airborne mechanical scanner thematic mapper simulator 
system flown from NASA Johnson Space Center during 1978 to 
1981 and now from NASA Ames Research Center in an NC-130 
aircraft. IFOV 2.5 mr, field of view 100°. Spectral bands: 0.45-0.52, 
0.52-0.60, 0.63-0.69, 0.76-0.90, 1.0-1.3, 1.55-1.73, 2.08-2.35, 
10.4-12.5 Jim. 
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OCM Ocean Color Monitor. A 13-band CCD pushbroom system proposed 

by the European Sp^^e Agency. Sun-synchronous orbit altitude 777 
km, equatorial crossing time 11 a.m., inclination 99°, repeat cover- 
age interval of 3 days. £FOV 800 m (with 200 m a design goal), field 
of view 80°, 11-bit quantization. Data rate 20.8 and 2.5 Mb/s for 
DFOVs of 200 and 800 m, respectively. Spectral bands: 0.39-0.41, 
0.435-0.455, 0.51-0.53, 0.555-0.575, 0.63-0.65, 0.675-0.695, 0.77- 
0.8, 0.99-1.05, 1.5-1.7, 3.5-3.9, 8.25-8.75, 10.3-11.3, 11.5-12.5 pm. 

OCS (I-m) Ocean Color Scanner. Three 10-band mechanical scanners flown in 
different configurations by NASA Goddard Space Flight Center on 
a U-2, a Lear Jet, and the Shuttle. [The Shuttle system is referred to 
as the Ocean Color Experiment (OCE) and uses eight slightly 
different bands than the OCS.] IFOV 3.5 mr, field of view 90°, 
10-bit qtianti^tion. Spectral bands: 0.455-0.475, 0.485-0.505, 
0.515-0.535, 0.546-0.566, 0.577-0.597, 0.609-0.629, 0.642-0.662, 
0.673-0.693, 0.727-0.749, and 0.763-0.783 pm modified to 0.498- 
0.518, 0.528-0.548, 0.560-0.580, 0.588-0.608, 0.620-0.640, 0.655- 
0.675, 0.683-0.703, 0.717-0.737, and 0.822-0.842 pm (Kim et al., 
1980). 


OCTS Ocean Color and Thermal Scanner, a Japanese mechanical scanner 

for a post-1992 proposed launch of MOS-2. IFOV 1 km, swath 1400 
km. Spectral bands: 0.43-0.45, 0.5-0.53, 0.54-0.58, 0.65-0.67, 
0.75-0.79, 10.5-11.5, 11.5-12.5 pm. 

Push-broom A 4-band airborne pushbroom scanner built by CNES. System can 
be rotated by ±26° to provide lateral stereo. Each band uses a 
sin^e 1728 detector linear array with a choice of lenses of focal 
lengths 55, 32, or 18 mm. IFOV 0.24, 0.41, or 0.72 mr, field of view 
22.6, 37.9, or 62.9°, respectively. 10-bit quantization. Selectable 
spectral bands in the range 0.4-1.1 pm. CNES uses a set that 
nominally matches that of the SPOT/HRV. Mid-IR bands are 
planned. (Laporte, 1981; Begni, 1983). 

Radarsat The primary payload proposed by the Canada Centre for Remote 

Sensing is a C-band synthetic aperture radar. The secondary payload 
proposed comprises three optical sensors: a high resolution (10 m) 
and a low resolution (1000 m) subsystem, and a cloud sensor. Only 
the former is described here. From an altitude of 1000 km the IFOV 
will be 10 m and the swath 100 km. The NEAp will be 0.1% of full 
scale with 10-bit quantization. A mirror will provide cross-track 
pointing to ±43°. Initial data rate is 448 Mb/s. The data will be 
compressed before transmission. The seven spectral bands tentar 
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lively proposed are: 0.43-0.45, 0.51-0.53, 0.54-0.56, 0.64-0.69, 
0.80-0.90, 1.55-1.75, 2.08-2.35 /xm (Cihlar et al., 1984). 

Return Beam Vidicon. Three differently filtered return beam vidi- 
con cameras were boresighted on Landsat-1 and -2 to provide 
muitispectral imagery. Only a small amount of data was collected on 
Landsat-1 because of an electrical difficulty, and it was found that 
the imagery from the MSS was preferred from a radiometric stand- 
point The muitispectral system was replaced on Landsat-3 by a 
panchromatic system of half the IFOV. The mudtispectrai RBV had 
an IFOV of about 70 m and a swath of 185 km. Spectral bands: 
0.475-0.575, 0.58-0.68, 0.69-0.83 /im (Landsat Data Users 
Handbook, 1979; Slater, 1980). 

A 13-band conical scanner flown as part of the Earth Resources 
Experiment Package on Skylab. Data were collected from May 25, 
1973, to February 8, 1974. Circular orbit altitude 435 km, inclina- 
tion 50®, repeat coverage every 5 days. IFOV 87 m, swath width 
11®. Spectral bands: 0.41-0.46, 0.46-0.51, 0.52-0.56, 0.56-0.61, 
0.62-0.67, 0.68-0.76. 0.78-0.88, 0.98-1.08, 1.09-1.19, 1.20-1.30, 
1.55-1.75, 2.10-2.35, 10.2-12.5 ;im (Due, 1982; Abel and Reynolds, 
1974). 

Shuttle Imaging Spectrometer Experiment, proposed by Jet Propul- 
sion Laboratory for a late 1980s Shuttle flight. The system has 128 
contiguous spectral bands in the 0.4 to 2.5-fim region with 10-nm 
spectral resolution in the 0.4 to 1.0-/xm region and 20-nm spectral 
resolution in the 1 to 2.5-/im region. IFOV 30-m, swath width 12.1 
km from 250-km altitude. ±20° cross-track and ±45° pointing 
along-track obtained by means of steerable mirrors. The system 
provides 8-bit quantixation and an initial data rate of 103 Mb/s 
reduced with on-board editing to 50 Mb/s (Wellman et al., 1983). 

Systeme Probatoire d’Observation de la Terre. The SPOT satellite is 
a multimission platform to be first laimched by an Aiiane-2 launcher 
in 1985 and placed in a circular sun-synchronous orbit of 832-km 
altitude and 98.7° inclination. Payload will consist of two HRV 
systems (see HRV). Equatorial crossing time will be 10:30 a.m. 

Shuttle Thermal Infrared Mxiltispectral Scanner. A pushbroom ver- 
sion of TIMS proposed by JPL as a Shuttle experiment for 1988. 
IFOV 30-m, swath 12-km, 8-bit quantization. Spectral bands: 8.2-8.6, 
8.6-9.0, 9.0-9.4, 9.4-10.2, 10.2-11.2, 11.2-12.2 /xm (Kahle, 1984). 
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A 4-band airborne mechanical scanner for the middle and thermal 
IR built by the Societe Anonyme de Tel^mmunication (SAT). 
Focal lengA 145 mm, EFOV 1.72 mr, field of view 90®. Analog 
output. Spectral bands: 3.6-4.3, 7.3-S.3, 8.4-9.5, 9.5-11.3 iim (Be- 
gni, 1983). 

Short Wave IniraRed Linear Array Pushbroom Radiometer. A 
sin^e-band aircraft system with three in-flight commandable spec- 
tral filtm used at NASA Goddard Space Fli^t Center for geologi- 
cal studies. IFOV 5 mr, field of view 18.5°. Spectral filters: 1.0-1.5, 
1.61-1.89, 2.14-2.30 Mm (Mason, 1984). 

Tropical Earth Resources Satellite, a joint Indonesia-Netherlands 
program. The sensor is a three spectral band system of 20-m IFOV 
and panchromatic band of 10-m IFOV with a swath of 100 km. 
High equatorial orbit, altitude 1680 km, 0° inclination, providing 

2.5- day repeat coverage. Use of a pointing mirror can increase 
frequency to a few repeated acquisitions each day, a capability that 
is desirable for regions fiequently covered by cloud. A thermal IR 
sensor has also been suggested to provide a 100-m IFOV over a 
100-km swath. NEAp 0.5%, NEAT 0.5 K at 300 K, 8-bit quantiza- 
tion. Spectral bands: 0.49-0.59, 0.61-0.685, 0.75-0.835, 10.4-12.5 
(tm Mritb a panchromatic band from 0.49-0.685 pm (Bunnik, 1984a). 

Thermal Infrared Multispectral Scanner. A 6-band mechanical 
scarmer delivered by Daedalus as the AAOS 1285 TIMS System to 
NASA-NSTL Earth Resources Laboratory in 1982. IFOV 2.5 mr, 
field of view 76.6°, 8-bit quantization. NEAT between 0.09 and 0.32 
K at 300 K, depending on spectral band. Spectral bands: 8.2-8.6, 

8.6- 9.0, 9.0-9.4, 9.4-10.2, 10.2-11.2, 11.2-12.2 pm (Kahle and 
Goetz, 1983). 

Thematic Mapper. A 7-band mechanical scanner orbited by NASA 
with MSS-4 and -5 as part of the Landsat-4 and -5 payloads in July 
1982 and March 1984, respectively. Sun-synchronous circular orbit, 
altitude 705 km, inclination 99.1°, equatorial crossing time 9.45 
xm., repeat coverage interval every 16 days. IFOV 30 ra (except for 
band 6, which is 120 m), swath 185 km, 8bit quantization, data rate 
85 Mb/s. Spectral bands: 0.45-0.52, 0.52-0.60, 0.63-0.69, 0.76-0.90, 
1.55-1.75, 2.09-^38, 10.4-12.5 pm (Engel and Weinstein, 1982; 
Slater, 1980). 

Thematic Mapper Simulator flown on a U-2 or an ER-2 aircraft at 
20-km altitude by NASA Ames Research Center. Of the 11 spectral 
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bands available, 7 are the same as those for the Thematic Mapper. 
IFOV 1.3 mr, field of view 43°, 8-bit quantization. 

VISSR Visible-Infrared Spin-Scan Radiometer. A 2-band spin-scaimed sys- 

tem in geostationary orbit providing images every 30 min. System is 
on the Synchronous Meteorological Satellites (SMS 1 and 2), the 
Geostationary Operational Environmental Satellites (COES 1, 2, 
and 3), and Japan’s Geostationary Meteorological Satellite (CMS 1 
and 2). EFOV 0.025x0.021 miad, field of view 18° x20°. Spectral 
bands: 0.55-0.75 and 10.5-12.5 /im (Due, 1982). 

VTIR Visible and Thermal Infrared Radiometer. A Japanese system to be 

part of the Marine Observation Satellite (see MOS) payload for 
measuring sea surface temperature. Visible band, 0.5-0.7 /xm, has 
IFOV of 0.9 km; three bands in the IR, 6-7, 10.5-11.5, and 
11.5-12.5 fim, have IFOV of 2.7 km. Swath width 500 km. NEAT 
< 0.5 K at 300 K (Ishizawa, 1981). 
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Abstract A method is described whereby held reflectance radiometers can be calibrated in 
an absolute sense using equipment available at most agricultural or enviromnental research 
locations. A radiometer is ptMitioned directly above a calibrated standard reflectance panel 
that is horizontal to the Earth's surface. The sun's direct beam is separated from the total by 
measuring the total, shading the panel with a nontransparent shield held between the sun 
and the panel, measuring the diffuse component, and subtracting the diffuse from the total. 
These measurements are repeated periodically from shortly after suiuise to near solar noon. 
A graph of the logarithm of the radiometer respoitse to the direct beam versus the secant of 
the solar zenith angle (known as a Langley plot) yields the spectral-extinction optical thick- 
ness of the atmosphere as the slope and the logarithm of the exoatmospheric irradiance 
divided by the calibration factor as the intercept. C^bration factors for two radiometers were 
within 10 percent of those obtained by other methods, indicating that this technique is a 
viable method for the absolute calibration of field radiometers. 


INTRODUCTION 

S MALL, UCHTWEICHT, RADIOMETERS that mea- 
sure radiation in the reflected solar portion of 
the electromagnetic spectrum are currently being 
used to obtain spectral data for environmental and 
agricultural research projects. For the most part, these 
radiometers are us^ in conjvmction with a refn- 
ence reflectance panel made horn flat metal plates 
coated with highly reflecting substances such as 
BaSO« (Robinson and Biehl, 1979) or Halon (Schutt 
et al., 1981). With a calibrated reflectance panel, the 
reflectance Actors of targets can be calculated by 
dividing the target radiance by the panel radiance, 
if the radiances from both sur^ces are measured at 
nearly the same time. For some purposes, however, 
it is necessary to quantitatively measure the radi- 
ance reflected from targets. For example, the cal- 
culation of the net amount of radiation absorbed by 
a plant canopy can be made if the irradiance at, and 
the reflected radiance from, the canopy are known 
(Jackson et al.. 1985). To accomplish this using re- 
flectance radiometers, calibration factors for each 
chaimel of the radiometer must be known. 

Calibration factors for some radiometers are avail- 
able from their manufacturer, whereas other radi- 
ometers are delivered uncalibrated. Calibration 
procedures are best carried out in well equipped 
optical laboratories. Such facilities are not ubiqui- 
tous, and are essentially unavailable to a number of 
researchers who routinely use small field radiome- 
ters. In the absence of a precise laboratory calibra- 
tion, a field calibration procedure would be of benefit. 
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This report describes a field procedure for the cali- 
bration of reflectance radiometers using equipment 
already available at most enviromnental research 
stations. 

METHOD 

The proposed method is an adaptation of a tech- 
nique for determining the solar constant from ground 
based measurements (Shaw et al., 1973; Slater, 1980). 
For any solar zenith angle (0J, the spectral irradi- 
ance on a surface perpendicular to the direct solar 
incident flux is. 


where is the exoatmospheric irradiance (the 
spectral irradiance outside the Earth's atmosphere 
on a plane one astronomical unit from the sun and 
perpendicular to the incidence flux), and t(\) is the 
spectral-extinction optical thickness at the wave- 
length \. 

Observations of E^ are made with a radiometer 
for several zenith angles from just after surmse to 
near solar noon on a single day. It is assumed that 
t(\) remains constant during the measurement pe- 
riod, and that the radiometer output signal is linear 
with respect to input radiant energy. The graph of 
ln(E;^) versus sec(9j), known as a Langley plot, will 
be linear if the assumptions are sufficiently met; The 
slope of the line is the spectral-extinction optical 
thickness, and the intercept is the natural logarithm 
of the exoatmospheric irradiance. The solar zenith 
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Fiq. 1. Solar spectral irradlance at the top of Itte atmosphere. 


angle can be directly measured, but it is readily cal- 
culated from a time measurement and the ephem- 
eris of the sun (List, 1958, Table 169). 

The iiradiance (E^) can be measur^ by pointing 
a radiometer, having a field of view that just in- 
cludes the solar disk, directly at the sun. The ii»- 
diance is the produce of an output voltage (V) times 
a calibratioi^ factor (c). 

For an uncalibrated radiometer. Equation 1 be- 
comes 

Vc = E^e''^^’^*’ (2) 

If ln(V) is plotted versus sec9^ the slope is — r(\) 
and the intercept is InfE^g/c). If E»„ is known, the 
calibration factor (c) is 

c = E,o/e^ (3) 

where A is the intercept of the Langley plot. 

Equation 2 requires that V' be in response to the 
direct solar incident flux. A radiometer having a field 
of view (FOV) of 3” pointed directly at the sun would 
see less than 2 percent diffuse radiation (Shaw et 
al., 1973). With field radiometers such as the Exo- 
tech Model lOO-A* and the Barnes 12-1000 Modular 
Multispectral Radiometer (MMR)*, the available Rev's 
are 1° and 15®. Using a 1® FOV, alignment of the 
radiometer optics with the solar beam would re- 
quire precision solar tracking equipment. Using a 
15® FOV would allow a significant portion of diffuse 
radiation to reach the detector. 

An alternate procedure is to measure the direct 


* Trade names and company names are included for the 
benefit of the reader and imply no endorsement of the 
product or company by the U. S. Department of Agricul- 
ture or the University of Arizona. 


radiant fltu reflected from a horizontal panel of 
known reflectance. This removes the stringent field 
of view requirement and allows radiometers having 
any field of view less than, say, 20® to be calibrated. 
Fields of view larger than about 20® would require, 
either an excessivdy large standard reflectance panel 
or placement of the ra&ometer close to the panel. 
The latter condition could cause errors by the in- 
strrunent blocking a significant amount of diffuse ’ 
radiation and/or re-reflecting radiation from the 
panel. 

The procedure depends upon the assumption that 
the standard reflectance panel is a perfect lamber- 
tian reflector. This requirement can be approxi- 
mated if the reflectance of the panel is accurately 
known at all illumination angles that may be en- 
countered during field measurements. 

The evaluation of E^^^, the corrections necessary 
for the reflectance panel to approximate a perfect 
lambertian reflector, and experimental procedures 
are discussed in the following sections. 

Exoatmospheric Irraoiance 

Data for solar spectral irradiance at the top of the 
atmosphere were obtained for the wavelength 
interval from 0.33 to 1.25 p.m from Neckel and Labs 
(1981), for 1.25 to 2.95 pm from Pierce and Allen 
(1977), and for wavelengths below 0.33 pm from 
Slater (1980). The data were interpolated to yield 
values for each nanometre between 0.25 and 3.0 p.m 
(Figure 1). 

Response functions for the four Exotech bands 
and the seven reflectance bands of the Barnes 12- 
1000 MMR are shown in Figure 2. Response function 
data are usually found in the radiometer instruction 
manual. The Exotech simulates the four Landsat 
multispectral scanner system (MSS) bands and the 
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Fk 3. 2. Spectral response functions for two field radiometers. 


MMR simulates the six solar reflective Landsat 
Thematic Mapper (TM) bands. The MMR also has a 
band at 1.15 to-L30 pm. Nominal wavelength 
intervals for the two instruments are given in Table 
1 . 

The exoatmospheric irradiance within the 
wavelength interval (the symbol will now apply 
to the wu) for each band of the two instruments 
was obtained by summing the product of the 
irradiance (Figure 1) and the response functions 
(Figure 2) for each radiometer band of interest. The 
results are presented in Table 1. The data are for 
the mean Earth-sim distance. varies with the 
square of the Earth-sun distance, oiffering by about 
7 percent between January and July. The values for 
E^„ given in Table 1 can be corrected to the actual 
Earth-sim distance by dividing them by the square 
of the radius vector of the Earth (r) for each 
measurement day. Tabular values of the radius vector 


are available (List, 1958, Table 169), but they can be 
approximated using the relation (Gurney and Hall, 
1983) 

r = 1 + 0.0167sin(2ir(D - 93.5)/365] (4) 

where D is the day of year. With this correction, 
and dividing by -ir to put the calibration factor in 
terms of radiance per volt. Equation 3 becomes 

c = E^(r^) (5) 

with c having units of W m** sr"* V*'. 

When ephemeris tables are used to calculate the 
solar zenith angle, the solar time must be noted at 
the prec^ time that the measurements are made. 
An uncertainty in the time measurement of about 
10 s would cause an error of about 0.1 percent in 
the zenith angle calculation for angles less than 75° 
(Thomason et al., 1982). 

Reflectance Panel 

Reflectance ^ctor data for a panel painted with 
BaSO« was provided by LARS, Purdue University (L. 
L Biehl, personal commurucation). Reflectance factor 
data for 10 and 20^ incidence angles were interpolated 
to yield data at 15° for each band of the Exot^ and 
the MMR (Table 1). Panel reflectance factors were 
measured as a function of incidence angle from 15 
to 75° using a device that allowed the panel to be 
held at a Imown angle to the sim's rays. With the 
MMR radiometer positioned about 1 m above and 
perpendicular to the panel, a measurement of the 
dir^ solar radiation was made (using a shading 
technique described in the following section). 
Subsequently, the incidence angle was changed and 


Table 1 . Nominal wavelength intervals (wu), Exoatmospheric irradiance within the specired wu of the 
Exotech Model ioo-a and the Barnes i2-iooo MMR radiometers, and BaS04 panel reflectance factors at 9; = 

15° FOR each wu. 


Radiometer 

Exotech 


Band 

Nominal 

wavelength 

(tun) 

Exoatmospheric 
irradiance 
(W m-2) 

Panel 

reflectance 

factor 

1 

0.5 -0.6 

167.6 

0.942 

2 

0.6 - 0.7 

188.2 

0.936 

3 

0.7 - 0.8 

157.3 

0.929 

4 

0.8 -1.1 

197.5 

0.914 

1 

0.45 - 0.52 

112.4 

0.948 

2 

0.52 - 0.60 

134.4 

0.941 

3 

0.63 - 0.69 

71.1 

0.935 

4 

0.76 - 0.90 

145.0 

0.922 

5 

1.15-1.30 

69.1 

0.897 

6 

1.55 - 1.75 

49.9 

0.855 

7 

2.05 - 2.30' 

22.0 

0.757 

3 

10.5 - 12.5 

- 

- 


Bames MMR 
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Fio. 3. Relative reflectance factors of a BaS 04 panel as a 
function of inddenca angle for two bands of a Held radi- 
ometer. 


the measurement repeated. Relative reflectance 
hictors were calculate as the ratio of radiometer 
output voltage to the voltage at 15°. The choice of 
15° as the the reference an^e was dictated by the 
fact that, for smaller incidence angles, the radiometer 
shadowed the panel. A measurement sequence for 
eight incidence angles required about 15 minutes. 

Values of relative reflectance hictors for the BaSO« 
panel as measured by MMR radiometer bands 1 and 

6 are shown in Figure 3. Data for MMR bands 2 to 5 
are intermediate to those shown, and those for MMR- 

7 are slightly higher than for MMR-6. Third degree 
pol 3 moirial equations were statistically fit to the data 
for each MMR band. The lines in Figure 3 indicate 
how well the data from two measurements were 
represented by the equations. Equations for MMR 
bands 2, 3, and 4 were used as representative for 
similar bands of the Exotech. 

The reflectance factor of the panel at any zenith 
angle R(0J is the product of the reflectance factor 
at 15° and the relative reflectance factor. In practice, 
the polynomial approximations of the relative 
reflectance factor data (Figure 3) were used in the 
calculations. 

Because the panel was kept horizontal during the 
measurements, instead of perpendicular to the solar 
beam as required by Equation 1, the voltage that 
would result if the standard reflectance panel were 
a perfect lambertian reflector held perpendicular to 
the sun's rays is approximated by 

V = V^/[R(0Jcos0j (6) 

which is the value of V to be used in the Langley 
plots. The term is the radiometer voltage when 
viewing a horizontal panel. 


If radiometer gain settings are other than 1 during 
measurements, the voltage should be adjusted to 
represent a gain of 1. llie output from the PbS 
detectors (B«mds 5 to 7) of the MMR are ambient 
temperature sensitive. A procedure to adjust the 
out]mt to a reference tempmture was described by 
Jackson and Robinson (1985). The Lan^ey plot for 
bands 5 to 7 of the MMR will not be linear if the 
temperature effect is not compensated for. 

Measurement of Direct Solar Radiation 

The direct component of the irradiance can be 
separated from the total by measuring the total, 
shading the panel with a nontransparent shield held 
between the sun and the panel, and measuring the 
diffuse component, with the direct component being 
the difference between the total and the diffuse. The 
use of a shading device (which should be held at as 
great a distance from the panel as practical) blocks 
a portion of the sky, reducing the amount of diffuse 
radiation that stril^ the panel by a small amount. 
The error caused by the use of the shade and the 
error caused by the time difference between 
measurements of the total and the disuse can be 
minimized by holding the shield so that the shaded 
area is to the side of the panel while the total 
irradiance is measured, then moved sideways to 
shade the panel for the diffuse measurements, then 
moved back to the original position for another total 
irradiance measurement. The before and after total 
irradiance measurements are averaged to yield a 
value corresponding to nearly the same time as that 
when the diffuse measurements were made. By 
holding the shield the same distance from the panel 
during ail measurements, the portion of the sky 
biock^ by the object is similar, thus reducing the 
error caus^ by its use. For a detailed discussion of 
this type of measurement see Che et al. (1985). 

Measurement Procedure 

An Exotech and a Barnes MMR (serial no. 119) 
were mounted side-by-side at the end of a rotatable 
boom which held the instruments about 1.7-m above 
a 1.2 by 1.2-m BaSO, painted horizontal panel. 
Outputs from the two radiometers were recorded 
using a portable data acquisition and storage device 
that also noted the time of measurement. 
Measurements were made on six dates, 19 November 
1983, 15 December 1983, 13 April 1984, and 8, 10, 
and 20 June, 1984, beginrting shortly after sunrise 
and continuing perio^cally until about an hour 
before solar noon. A 1.3 by 1.3-m flat shield mounted 
at the end of a 3.5-m pole was used to shade the 
panel during measurements of diffuse radiation. The 
measurement sequence was total (sunlit panel), 
diffuse (shaded panel), and total. 

When the measurements were completed, the data 
were dowtrloaded to a computer. The voltage data 
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Fig. 4. Langley plots for the four Exotech bands. Oats were 
taken on 13 April 1984. 



Fig. 6. Langley plots for mmr bands S to 7. Data were taken 
on 19 November 1983. 



Fig. 5. Langley plots for mmr bands 1 to 4. Data were taken 
on 15 December 1983. 


were, when necessary, adjusted to a gain of 1 and 
the ambient temperature sensitive channels adjusted 
to a reference temperature of 25®C (Jackson and 
Robinson, 1985). For each channel on each radiometer 
at each measurement time, the two total (sunlit) data 
were averaged and the diffuse subtracted to yield a 
voltage value (V^^) due to direct radiation at the time 
of the diffuse measurement (the time at which 8^ 
was calculated). Using 9„ the relative reflectance 
' factor was calculated using the appropriate 
polynomial equations (see discussion of Figure 3), 
the panel reflectance factors at 15® were taken from 
• Table 1, and the voltage was adjusted to V using 
Equation 6. Plots of ln(V) versus secfO,) were made 


to obtain the intercept A which, when used in 
Equation 5, yielded the calibration Actors. 

RESULTS AND DISCUSSION 

Representative Langley plots for eleven bands and 
for three measurement dates are shown in Figures 
4 to 6. In all cases the data are linear, a necessary 
but not sufficient condition for the underlying as- 
sumptions to be met. Reagan et al. (1984) showed 
by model simulations that a linear relation can result 
when the spectral-extinction optical thickness varies 
temporally. This emphasizes the necessity of mak- 
ing measurements under clear, stable weather con- 
ditions. 

If the temperature correction had not been made 
to the PbS detector outputs for MMR bands 5 to 7, 
the data would have been markedly non-linear, being 
concave downward. 

The calibration factors for the two instruments are 
given in Tables 2 and 3. Table 2 includes the cali- 
bration factors for the Exotech as furnished by the 
manufacturer. The percent difference between the 
company measured and the Langley plot results 
range from -4.0 percent for band 2 to 9.6 percent 
for band 3. There were no data available with which 
the MMR calibration factors (Table 3) could be di- 
rectly compared. 

A secondary comparison was made with a second 
MMR (S# 116) at White Sands, New Mexico on 28 
October 1984. Langley plots of the data provided a 
caiibrarion factor for the first four MMR bands. These 
data, along with the means for bands 1 to 4 of MMR 
S# 119 from Table 3, are presented in Table 4. Also 
included are calibrations for bands 1 and 3 for MMR 
S# 116 from a laboratory-based procedure (Phillips, 
1985). The detector temperature was not monitored 
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Tabi£ 2. Caubration factors for the Exotech Model lOO-A determined from Langley plots and response 

FUNCTIONS GIVEN IN THE MANUFACTURERS INSTRUCTION MANUAL. AlL CAUBHATION FACTORS HAVE UNITS OF 

Wm-2 SR-' V-'. 


Date of 


Band 



measurement 

1 

2 

3 

4 

19 Nov 83 

U.5 

11.8 

13.6 

13.5 

15 Dec 83 

12.1 

11.5 

13.3 

13.4 

13 Apr 84 

12.8 

12.1 

13.6 

13.1 

08 Jun 34 

12.0 

11.6 

13.3 

13.2 

10 Jun 84 

13.2 

U.5 

14.1 

13.9 

20 Jun 34 

13.6 

U.7 

14.2 

13.7 

Mean (SD) 

12.7 (0.6) 

12.0 (0.5) 

13.7 (0.4) 

13.5 (0.3) 

Company calibration 

12.7 

12.5 

US 

U9 

Percent difference 

0 

-4.0 

9.6 

4.7 


Table 3. Caubration factors for the Barnes MMR 12-1000 S#1 19 determined from Langley plots and 
RESPONSE functions GIVEN IN THE MANUFACTURERS INSTRUCTION MANUAL. ALL CAUBRATION FACTORS HAVE UNITS OF 

W M-^ 3R-’ V-’. 

Date of 522^ 

measurement 1 2 3 4 5 6 7 


a 

a 

a 

a 

I 


Table 4. Comparison of caubration factors for 4 bands of two mmr radiometers and reld and laboratory 

CAU BRATTONS FOR 1 BANOS OF ONE RADIOMETER. THE FACTORS HAVE UNITS OF W/M'^ SR'' V*'. 

Date of 

measurement 1 


*The internal gain adjust was changed from the factory setting. 

for S# 116, precluding the calculation of calibration The spectral-extinction optical thickness (t), being 
factors for bands 5 to 7. On this radiometer, the gain the slope of the Langley plots, is readily obtained 

on band 2 had been internally adjusted from the from the same data as the calibration factors. This 

factory setting, invalidating any comparison of that parameter is a measure of the extent to which the 

band. For the following comparison we make the atmosphere scatters and absorbs the irradiance within 

• assumption that, because S# 116 and S# 119 were the spectral bands of the MMR. Values of t for each 

made at the same time, their calibration factors are of the eleven bands are given in Table 5. Differences 

similar (for bands 1, 3, and 4). between the six measurement dates are evident, with 

The calibration factors for these two instruments 8 June 1984 having the highest values. Two days 

differed by 12 percent in band 1, by 3.6 percent in later the next to lowest values for the sue days were 

band 3, and were identical for band 4. The factors recorded. Although both days were cloud-free, 10 ’ 

for both instruments are reasonably close to the lab- June apparently had less atmospheric contaminants 

# oratory-based calibration of bands 1 and 3. These than did 8 June. Che et al. (1985) presented a de- 
data support the premise that the Langley plot tailed discussion of the measurement of t using a • 

method is a viable means of obtaining calibration field radiometer and a calibrated reflectance panel, 

factors using field data. The data and the standard deviations of the cal- 


SN#119 I 

SN#116 (field) 

SN#116 (lab) V 


Band 

; 2 3 i_ 

.0 7.00 8.00 17.0 

.8 13.5* 7.72 17.0 

.3 - 7.57 


45 

41 

40 

41 
48 
44 

43 

03 


19 Nov 83 

10.3 

7.16 

7.77 

17.0 

6.78 

7.26 

15 Dec 83 

10.4 

6.91 

7.58 

16.7 

6.72 

7.00 

13 Apr 84 

11.2 

7.04 

8.09 

17.0 

6.59 

6.94 

08 Jun 84 

10.4 

6.62 

7.83 

16.7 

6.91 

7.10 

10 Jun 34 

11.4 

7.13 

8.37 

17.4 

7.26 

7.39 

20 Jun 84 

11.5 

7.16 

8.34 

17.2 

7.00 

7.23 

Mean 

11.0 

7.00 

8.00 

17.0 

6.88 

7.15 

SD 

0.5 

0.21 

0.32 

0.3 

0.24 

0.17 
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Ta81£ 5. Values of the spECTBw.-exTO»cTtON optical thickness (t) osTAiNaj ftom slopes of Langley plots fop 
^VEN SANDS OF THE BaRNES MMR S#1 19 AND TNE FOUR ExOTECH BANOS. 


Date of 
measurement 




Band 




1 

2 

3 

4 

5 

6 

7 

Barnes MMR 








19 Nov 83 

0.211 

0.169 

0.105 

0.074 

0.(B9 

0.025 

0.050 

15 Dec 83 

0.243 

0.197 

0.122 

0.095 

0.086 . 

0.049 

0.087 

l5Apt84 . 

0.286 

0.239 

0.157 

0.122 

0.094 

0.054 

0.072 

OB Jun 84 

0.353 

0.292 

0.195 

0.149 

0.116 

0.067 

0.095 

lO Jun 84 

0.225 

0.185 

0.112 

0.089 

0.074 

0.031 

0.063 

20 Jun 84 

Q.233 

0.19Q 

0.U7 

0.089 

0.079 

0.039 

0.063 

Exotech 








l9 Nov 83 

0.165 

0.123 

0.078 

0.072 




15 Dec 83 

0.187 

0.140 

0.0% 

0.095 




13 Apr 84 

0.232 

0.179 

0.127 

0.120 




05 Jun 84 

0.298 

0.229 

0.164 

0.153 




lO Jun 34 

0.181 

0.133 

0.092 

0.089 




20 Jun 34 

0:183 

0.136 

0.093 

0.090 





ibiation factors given in Tables 2 and 3 indicate that 
the method is reasonably precise. The variation be- 
tween the six measurements does not appear to be 
related to time nor to the degree of atmospheric 
* scattering (Table 5). 

The accuracy of the measurement depends upon 
the accuracy of the exoatmospheric irradiance ^ta 
and the reflectance data for the standard panel. Er- 
rors in will cause a proportional error in c. The 
most probable error in exoatmospheric irradiance 
•within the wavelength interval as given in Table I 
would be in the response functions used in the cal- 
culation of the values. These functions could be de- 
termined for a particular instrument using equipment 
available in most optics laboratories. 

If the reflectances of the standard panel are in 
error by a multiplicative amount, the calibration fac- 
tors will be in error by the negative of that amoimt. 
Thus, if the relative reflectance of the panel is known 
with sufficient accuracy, but the reference reflec- 
tance at 15° is in error, the calibration factors can be 
easily corrected when the true reflectance values are 
known. If the relative reflectances are in error, the 
Langley plots may not be linear. 

Equations 1 and 2 are expressions of Beer's law 
which was derived for monochromatic radiation. 
Thomason et al. (1982) showed that less than 0.1 
percent error would result from using a 0.01-um 
wavelength interval. The wavelength intervals used 
ranged horn 0.06 to 1.1 urn. However, a twenty- 
fold larger error (2 percent) can be tolerated in the 
field method discussed here. 

CONCLUDING REMARKS 

The Langley plot technique is a viable field method 
for absolute radiometric calibration at the 10 percent 
level. Requirements are that the reflectance of a 
standard panel be known as a function of incidence 


angle, and that the exoatmospheric irradiance be 
known for the wavelength interval of each channel 
of the radiometer. The method is rather simple to 
implement but requires several hours to obtain one 
set of calibration Actors. 
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Forum 


Comments on Evaluations of Simulated SPOT Data 


The August 1985 issue of Phatogramnutric Engi- 
neering and Remote Sensing contained a number of 
arddes evaluating products from the 1983 SPOT Sim- 
ulation Campaign. Because condusions from the 
simulation pertain also to the utility of the satellite 
data, it is appropriate to clarify differences between 
the simulation and satellite cUta which may afiect 
some applications. 

The paper by Ackleson and Klemas discusses the 
utility ot the simulation data for disciiininating water 
masses, emphasizing the low noise content of the 
simulation data. However, the simulation data have 
lower noise content than is to be expected horn the 
satellite data; e.g.. Saint and Weill (1984) present 
noise values for the three 20-m channels (SI, S2, 
and S3) and the panchromatic channel (P). 


Radiometric Noise 

SI 

S2 

S3 

P 

SPOT (Nominal) 
NeAL 

0.14 

0.15 

0.11 

0.15 

Simulation 

NeAL 

0.05 

0.05 

0.03 

0.05 

Radiometric noise is given 

in mW cm*^r‘|juir') 


"The simulatiotts are in this respect significantly 
better than real data, but this can be taken into ac- 
count by using the calibration coeffidents on the 
CCTs to compute the noise equivalent digital count 
corresponding to the noise and quantization levels 
of SPOT." This computation requires knowledge of 
the satellite calibration coeffidents, as describe by 
Price (1984). However, in the author's experience 
with two scenes, the calibration coeffidents are not 
correct in the simulation data sets. Thus, the con- 
version to satellite equivalent data is problematic. 


The difference between simulation and satellite data 
is expected to be most significant in cases similar to 
that studied by Ackleson and Klemas having low 
surface refiectivities. , 

In addition, Ackleson and Klemas state that "It is 
anticipated that banding will not be a problem within . 
the operational SPOT data because of ffie linear anay , 
configuration of the sensor." However, the simu- 
lation data again present an optimistic piettue, as 
the simulation data were obtained by a Daedalus 
scanner with a single detector per channel, with 
scaniung by means of a rotating minor, while the 
SPOT satellite sensor will have some thousands of 
detectors. It is likely that minor residual banding 
will remain after c^bration, unless special tech- 
niques are used (Bernstein et al.. 1985). For data 
which have not been geometrically corrected, this 
effect will be present as vertical striping associated 
with the pushbroom scan of the sensor along the 
satellite track. 

Several authors, e.g., DeGloria, have carried out 
photointerpretation of the imagery produced in 
conjimction with the SPOT Simulation Campaign. 
Several caveats apply to the availability of such sat- 
ellite products: 

• Although the simulation imagery included 10-metie 
mulfispectral image products produced (tom the high 
resolution aircraft data, the 10- and 20-metre satellite 
data will not be coregistered by the satellite instru- 
ments. This task must fall to the user or to the SPOT 
Corporation, at least until the launch of SPOT 3 and 
4, which will acquire 10- and 20-m data in registered 
form. In contrast, the planned Landsats 6 and 7 will , 
be able to acquire 15-m panchromatic data coregis- 
tered with the mulfispectral 30-m data. 

(continued on page 211) ^ 


